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UNCLASSIFIED  ABSTRACT 


(U)  The  objective  of  the  engineering  design  study  of  the  close  air  support 
weapon  (CASW)  was  to  provide  design  considerations  for  the  new  close  air 
support  missile  (CASM).  The  derivation  of  the  missile  was  undertaken  based 
on  the  modification  of  an  existing  missile.  This  study  incorporates  opera¬ 
tional  requirement  and  warhead  effectiveness  studies  for  various  close  air 
support  targets  leading  to  warhead  and  launch  envelope  recommendations. 

A  thorough  analysis  of  the  system  performance  and  terminal  accuracy  was 
conducted.  Missile  simulation  models  and  a  system  description,  including 
missile,  launcher,  avionics,  arid  aerospace  ground  equipment  (AGE)  are  pro¬ 
vided.  A  cost  analysis  exercise  was  conducted  for  the  design,  development, 
test  and  evaluation  (DDTfcE)  and  production  of  the  candidate  approach.  This 
report  consists  of  six  volumes:  Management  Summary,  Operational  Analysis 
and  Warhead  Effectiveness,  System  Analysis,  System  Design,  Cost  Analysis, 
and  Missile  Simulation. 
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SECTION  I 


CLOSE  AIR  SUPPORT  MISSILE  SIMULATIONS 


1.1  INTRODUCTION 

(U)  The  basic  simulation  tools  used  for  the  CAS  weapon  system  analysis 
included  (1)  a  basic  six-degree-of-freedom  simulation  digital  program, 

(2)  a  modified  version  of  the  six-degree-of-freedom  including  a  Monte  Carlo 
version,  and  (3)  and  a  simplified  adjoint  system  model  (described  in  Vol¬ 
ume  III,  System  Analysis). 


1.2  SIX-DEGREE-OF-FREEDOM  DIGITAL  SIMULATION 

(U)  The  basic  simulation  program  used  in  the  performance  analysis  eval¬ 
uation  of  the  close  air  support  missile  system  was  a  six-degree-of-freedom 
digital  computer  program  which  has  been  constructed  by  modifying  the  simu¬ 
lation  of  the  AGM-65A.  The  objective  of  the  simulation  was  to  provide  a 
complete  and  intensively  detailed  representation  of  the  entire  missile  system 
which  could  be  used  for  final  design  verification,  performance  evaluation, 
and  spot  checks  of  parameter  optimization  results  from  simpler  simulations. 

(U)  The  present  simulation  represents  a  highly  sophisticated  and  analyti¬ 
cal  model  of  the  entire  missile  system.  The  simulation  was  developed  using 
the  system  analysis  by  digital  simulation  using  analog  methods  (SADSAM) 
programming  system.  This  system  is  used  with  the  FORTRAN  IV  compiler 
language  and  provides  an  extensive  library  of  functions  and  operations  which 
lend  themselves  well  to  handling  the  bookkeeping  and  computational  problems 
of  engineering  systems.  The  computational  speed  of  a  simulation  developed 
around  this  system  is  much  greater  than  one  programmed  in  a  more  conven¬ 
tional  manner. 

(U)  Wide  use  has  been  made  of  the  6  DOF  program  in  evaluating  the  sys¬ 
tem  performance,  especially  as  an  analytical  tool  in  defining  the  miss 
weighting  function  as  affected  by  heading  error,  launch  velocity,  target 
motion,  motor  temperature  effects,  and  seeker  drift  effects. 

(U)  A  listing  of  the  basic  6  DOF  programs  used  in  the  study  together  with 
system  nomenclature  and  input  data  coefficients  is  presented  herein.  This 
document  represents  a  complete  and  comprehensive  description  of  the  6  DOF 
program,  including: 

1)  Program  listing 

2)  Mathematical  model  description 

a)  Block  diagram 

b)  Parameter  definitions 
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c)  Transfer  functions  where  applicable 

d)  Program  input  requirements 

e)  Program  flow  charts 

(U)  Figure  1  illustrates  a  simplified  block  diagram  of  the  entire  simula¬ 
tion  model  and  indicates  the  depth  and  scope  that  have  been  included  in  this 
simulation  package. 

(U)  One  of  the  specific  laser  seeker  model  capabilities  include  the  ability 
to  evaluate  the  effects  of  laser  spot  size  image  variations  which  cause  a 
variation  in  the  angle-tracking  loop  gain.  Figure  2  shows  a  typical  seeker 
gain  curve  varying  with  range  and  spot  size  that  has  been  modeled.  The  » 

compensation  networks,  as  indicated  in  Figure  1,  can  be  placed  in  the 
seeker  forward  loop  to  increase  the  stability  margin  and  to  reduce  degrada¬ 
tion  in  angle-tracking  response  resulting  from  the  effects  of  spot  size  growth. 


1.3  MONTE  CARLO  SIMULATION 

(U)  Paralleling  the  approach  used  for  AGM-65A  performance  evaluation, 
a  six-degree-of-freedom  digital  computer  simulation  incorporating  Monte 
Carlo  techniques  has  been  developed  for  the  close  air  support  missile  con¬ 
cept  formulation  study.  For  any  given  set  of  launch  conditions  against  a 
particular  target,  there  will  be  some  statistical  variation  of  the  miss  dis¬ 
tance.  This  results  from  target  designation  errors,  missile  parameter 
tolerances,  uncertainties  in  the  launch  conditions,  and  uncertainties  in 
ambient  flight  conditions  such  as  temperature  and  winds. 

(U)  A  statistical  description  of  each  parameter  is  stored  in  the  computer 
and  sampled  by  a  Monte  Carlo  process  which  randomly  selects  a  value  of 
each  parameter  within  its  own  distribution.  A  homing  encounter  is  then  run 
with  this  set  of  parameters  and  with  all  error  sources  present,  including 
tracker  noise,  seeker  drift  limit  cycle,  steering  unbalance,  and  target 
motion.  This  results  in  a  miss  distance.  The  process  is  repeated  many 
times,  each  with  a  new  set  of  parameters  selected  by  the  Monte  Carlo  tech¬ 
nique.  As  a  result,  a  distribution  of  impact  points  about  the  aimpoint  is 
obtained,  as  illustrated  in  Figure  3. 

(U)  The  Monte  Carlo  simulation  was  primarily  used  to  verify  results 
obtained  from  the  adjoint  system  simulation.  The  Monte  Carlo  simulation 
was  used  to  check  the  wind  disturbance  effects  obtained  from  the  adjoint. 

It  has  also  been  used  to  simulate  and  verify  random  noise  disturbance 
occurring  from  tracker  noise  or  designation  signal  noise. 

(U)  The  results  from  these  simulations  are  reported  in  Volume  III, 
System  Analysis. 
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Figure  1.  Close  Air  Support  Missile  Six- Degree  Of-F  reedom  Simulation,  Simplified  Block  Diagram 


RATIO  OF  SEEKER  RANGE  TO  SPOT  DIAMETER,  R  /De 
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Figure  2.  Forward  Loop  Gain  Variation  Model 


Figure  3.  Designation  and  Impact  Point  Dispersions 
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SECTION  II 


SK-DEGREE-OF-FREEDOM  SYSTEM  PERFORMANCE 

DIGITAL  MODEL 


2.  1  OBJECTIVES,  CONCEPTS,  REQUIREMENTS,  METHODS,  AND 
TECHNIQUES 

2.  1.  1  Objectives 

(U)  The  simulation  objectives  are  to  provide  a  complete  and  highly 
detailed  simulation  of  the  entire  missile  system  which  can  be  used  for  final 
design  verification  and  spot  checks  of  parameter  optimization  results  from 
simpler  simulations. 

2.1.2  Concepts 

(U)  This  simulation  represents  the  most  complete  all-analytical  model 
of  the  entire  CAS  missile  system.  A  bare  minimum  of  simplifying  assump¬ 
tions  are  made  in  subsystem  hardware  representation.  Seeker  drift  phe¬ 
nomena  are  included  in  their  entirety.  Autopilot  transfer  functions  are  not 
approximated  in  any  way.  In  short,  this  simulation  is  the  best  overall 
missile  system  dynamics  model. 

2.1.3  Requirements 

(U)  The  operation  of  this  program  requires  the. use  of  an  object  deck,  a 
data  deck,  and  a  library  of  special  functions  which  constitutes  a  special 
computational  package  called  SADSAM.  All  these  inputs  must  be  compatible 
for  use  with  the  GE635  Computer. 

2.  1.4  Methods  and  Techniques 

(U)  As  previously  mentioned,  the  simulation  was  developed  using  a  pro¬ 
gramming  system  called  SADSAM.  This  system  is  used  with  the  FORTRANIV 
compiler  language  and  provides  an  extensive  library  of  functions  and  opera¬ 
tions  which  lend  themselves  well  to  handling  the  bookkeeping  and  computa¬ 
tional  problems  of  many  engineering  systems.  Typical  of  the  available 
functions  or  operations  are  integrations,  differentiations,  linear  transfer 
functions  of  almost  any  order,  orthogonal  transformations,  Euler  angle 
computations,  and  a  variety  of  non-linear  operations  such  as  limiting  and 
backlash.  The  computational  speed  of  a  simulation  developed  around  this 
system  is  much  greater  than  one  programmed  in  a  more  conventional 
manner. 
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2.  2  DESCRIPTION  OF  WEAPON  SYSTEM  EQUIPMENT  USED 
2.  2.  1  Introduction 


(U)  The  dynamic  performance  of  the  CAS  guided  missile  has  been  simu¬ 
lated  by  a  digital  computer  model.  The  purpose  of  this  section  is  to  describe 
this  model. 

(U)  The  digital  computer  simulation  model  is  organized  in  modules  and 
written  in  FORTRAN  IV  compiler  language.  The  CAS  simulation  model 
makes  use  of  SADSAM  III,  a  programming  technique  which  is  specifically 
intended  for  dynamic  simulations  and  which  achieves  both  high  dynamic 
accuracy  and  high  speed  operation.  In  addition,  it  provides  preprogrammed 
subroutines  typical  of  those  used  in  missile  simulations,  as  an  aid  to  the 
analyst. 

(U)  The  CAS  simulation  program  is  organized  into  four  principal  modules: 

(1)  universal  seeker,  (2)  autopilot,  (3)  control  surface,  and  (4)  aerodyna¬ 
mics  modules.  In  addition,  there  are  six  other  minor  or  supporting  modules 
incorporated  into  the  simulation.  These  are  (1)  initial  conditions  computer, 

(2)  aimpoint  wander  (target  motion),  (3)  angle  restoration  bias  (guidance 
law),  (4)  blind  range  filter,  (5)  track,  and  (6)  gyro.  Each  of  the  eight  sub¬ 
routines  and  the  main  program  are  discussed  in  the  paragraphs  that  follow. 
Also,  there  is  a  description  of  how  the  model  is  used  as  well  as  a  descrip¬ 
tion  of  the  SADSAM  programming  technique. 

(U)  Paragraph  2.  2.  2  provides  a  brief  overall  description  of  the  CAS 
model. 

2.  2.  2  Background 

(U)  The  dynamic  model  of  the  CAS  missile  describes  the  motion  of  the 
missile  in  three  dimensions  and  makes  use  of  all  six  degrees  of  freedom: 
three  positions,  three  velocities,  three  attitude  angles,  and  three  angular 
rates.  The  two  vector  equations  (translation  and  rotation)  applying  Newton's 
Second  Law  to  the  rigid  missile  are  rigorously  applied,  and  the  kinetic  and 
kinematic  behavior  of  the  gimballed  seeker  is  also  described  in  great  detail. 
All  significant  contributors  to  seeker  drift  are  represented  as  well  as  all 
significant  aerodynamic  forces  and  moments. 

(U)  The  program  is  arranged  in  four  basic  modules  describing  the  seeker, 
the  autopilot,  the  control  surfaces,  and  the  aerodynamics.  The  modular 
representation  is  used  because  (1)  it  .permits  the  model  to  be  programmed 
and  checked  out  more  easily;  (2)  it  permits  changes  and  substitutions  to  be 
made  more  easily;  and  (3)  it  provides  a  good  correspondence  with  the  actual 
hardware  elements  of  the  missile,  so  that  subsystem  specialists  can  parti¬ 
cipate  in  the  performance  evaluation  process  in  a  more  direct  fashion. 

(U)  The  simulation  is  exercised  by  the  operator  inputing  values  of  all  the 
system  constants  and  setting  the  initial  conditions  of  the  system  to  the 
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desired  values.  Due  to  the  extensiveness  of  the  simulation,  numerous  inputs 
are  required  for  a  complete  initial  conditions  set.  These  are  provided  as 
the  output  of  the  initial  conditions  subroutine  in  a  form  convenient  for  use 
in  the  simulation.  The  input  required  for  this  subroutine  establishes  the 
missile  configuration  at  the  time  of  launch  and  is  fully  described  in  para¬ 
graph  2.  3.  3.  The  simulation  begins  with  all  dynamic  elements  at  their 
steady  state  conditions.  The  printout  interval  and  maximum  simulation 
time  is  also  input  for  each  run,  and  at  the  option  of  the  operator,  any  or 
all  the  system  variables  may  be  printed  in  any  sequence. 


2.  3  DESCRIPTION  OF  DATA  USED  IN  SIMULATION 
2.  3.  1  Coordinate  Systems 

(U)  Four  different  coordinate  systems  are  in  use  in  the  simulation: 

(1)  earth,  (2)  missile  body  (control  surface),  (3)  autopilot,  and  (4)  seeker. 
Earth  coordinates  are  simply  fixed  in  inertial  space  with  the  missile  located 
at  the  origin.  The  Z-axis  is  vertical  downward,  and  the  X-axis  is  aligned 
with  the  ground  projection  of  the  original  line-of-sight  vector.  The  missile 
body  axes  are  fixed  in  the  missile  with  the  X-axis  aligned  longitudinally  and 
the  Y-  and  Z-axes  aligned  with  the  control  surfaces.  Since  the  control 
surface  orientation  is  nominally  at  45- degree  angles  with  the  horizontal  and 
vertical,  these  axes  are  also  rotated  in  this  manner.  The  autopilot  axes 
are  also  fixed  in  the  missile  body  with  the  X-axis  oriented  longitudinally 
but  with  the  Y-  and  Z-axes  rotated  45  degrees  from  the  missile  body  axes. 
The  seeker  coordinates  are  fixed  to  the  seeker  head  and  are  aligned  with 
the  autopilot  axes  when  the  seeker  gimbal  angles  are  set  to  zero. 

(U)  A  fifth  set  of  coordinate  axes  is  used  in  the  aerodynamics  calcula¬ 
tions.  These  are  the  maneuver  axes  which  are  obtained  by  rotating  the 
missile  body  axes  about  the  X-axis  through  the  aerodynamic  roll  angle. 

In  addition,  a  sixth  coordinate  set  is  also  used  for  miss  distance  calcula¬ 
tions.  This  miss  distance  coordinate  set  may  be  obtained  from  the  inertial 
coordinate  set  simply  by  a  rotation  about  the  Y-axis  which  aligns  the  X-axis 
with  the  initial  line  of  sight.  When  miss  distance  is  measured  in  this 
coordinate  set,  it  is  taken  as  the  missile-to-target  distance  at  the  point 
where  the  x- component  of  range  reduces  to  zero. 

(U)  The  coordinate  sets  described  above  are  listed  for  convenience  in 
Table  I.  Figure  4  depicts  the  Euler  angle  relationships  among  the  various 
sets  by  means  of  piograms  (or  resolver  chains).  With  the  exception  of 
earth-fixed  coordinates,  these  coordinate  systems  are  also  shown  graphic¬ 
ally  in  Figure  5. 

2.  3.  2  Main  Program 

(U)  Within  the  simulation,  each  subroutine  or  module  deals,  for  the 
most  part,  with  only  a  single  set  of  coordinates.  Transformations  between 
these  coordinates  are  therefore  performed  largely  in  the  main  or  call 
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TABLE  I.  SIMULATION  COORDINATE  SYSTEMS 


Coordinate  System 

Desc  ription 

Inertial  (Earth) 

Fixed  in  earth  with  the  origin  at  initial 
missile  location.  X-axis  is  horizontal  and 
aligned  with  ground  projection  of  initial  line 
of  sight  and  positive  in  direction  toward 
target.  Z-axis  is  vertical  and  positive  down¬ 
ward,  Y-axis  horizontal  and  positive  in  the 
sense  to  complete  a  right-handed  system. 

Miss  Distance 

Fixed  in  earth  with  X-axis  aligned  with 
original  line  of  sight  and  positive  in  direction 
toward  target.  This  coordinate  set  is 
obtained  by  a  rotation  of  the  earth 
coordinates  about  the  Y-axis. 

Autopilot 

Fixed  in  the  missile  body  with  Y-  and  Z-axes 
at  45  degrees  to  the  planes  of  the  control 
surfaces,  and  with  the  X-axis  in  the  longitu¬ 
dinal  axis  of  the  missile,  positive  in  the 
direction  of  flight.  This  set  is  related  to  the 
inertial  system  by  three  Euler  rotations  in 
the  following  sequence: 

Earth- Yaw -Roll -Pitch -Autopilot 

The  positive  sense  of  these  rotations  is  the 
same  as  the  positive  sense  of  the  axes  about 
which  the  rotations  take  place. 

Missile 

Fixed  in  the  missile  body  with  the  Y-  and 

Z-axes  in  the  planes  of  the  control  surfaces. 

This  coordinate  set  is  obtained  by  a  rotation 
of  the  autopilot  axes  through  45  degrees 
about  the  positive  X-axis. 

Seeker 

Fixed  to  the  seeker  head  with  the  X-axis 
aligned  with  the  boresight.  This  set  is 
related  to  the  autopilot  system  by  two  gimbal 
rotations  in  the  following  sequence: 

Autopilot -Elevation -Azimuth -Seeker 

The  gimbal  rotations,  elevation,  and  azimuth 
are  taken  about  the  nominal  autopilot  Y-  and 
Z-axes,  respectively,  with  the  positive  sense 
of  rotation  being  the  same  as  that  of  the  axis 
about  which  it  takes  place. 
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TABLE  I.  SIMULATION  COORDINATE  SYSTEMS  (CONCLUDED) 


Coordinate  System 

Dose  ription 

Maneuver 

This  coordinate  system  is  related  to  the 
missile  coordinate  system  but  is  not  fixed  in 
the  missile  body.  The  X-axis  is  aligned  with 
that  of  the  missile  set,  but  the  Z-axis  is 
selected  so  that  the  missile  velocity  vector 
lies  in  the  XZ-plane.  The  direction  of  the 
lateral  component  of  missile  velocity  fixes 
the  positive  direction  of  the  Z-plane.  The 
angle  through  which  the  missile  axes  must  be 
rotated  about  the  negative  X-axis  to  coincide 
with  the  maneuver  axes  is  called  the  aero¬ 
dynamic  roll  angle,  0a.  When  no  lateral 
component  of  velocity  exists,  0a  is  taken  to  be 

45  degrees. 

NOTE: 

All  systems  are  in  right-handed  rectangular  cartesian 
coordinates. 

program  which  serves  primarily  to  direct  signal  flows  among  the  four  func¬ 
tional  subroutines  and  to  permit  the  input  and  output  of  data. 

(U)  A  FORTRAN  listing  of  the  call  program  appears  in  Table  II.  The 
flow  chart  and  block  diagram  for  this  program  appear  in  Figures  6  and  7, 
respectively.  Table  in  is  a  glossary  of  the  terms  used  in  this  program. 
This  includes  the  dimensions  and  coordinate  systems  referred  to  as  well  as 
the  subscripted  variable  or  constant  number  used  to  identify  the  term. 

2.  3.  3  Initial  Condition  Subroutine  (Setic) 

(U)  The  primary  purpose  of  this  subroutine  is  to  accept  the  data  which 
specifies  the  missile  conditions  at  the  time  of  launch  and  to  convert  this 
data  into  initial  conditions  useable  by  the  simulation.  Since  the  values  of 
certain  system  parameters  are  also  subject  to  change  over  a  series  of 
simulated  trajectories,  this  subroutine  also  provides  for  a  common  area  of 
data  input  shared  by  these  parameters  and  the  initial  conditions.  This 
common  area  is  the  T-array  provided  by  the  SADSAM  system.  The  inputs 
and  outputs  of  this  subroutine  are  shown  in  Tables  IV  and  V,  respectively. 
The  physical  relationships  between  the  various  input  quantities  are  indicated 
in  Figures  8  through  10.  The  definition  of  the  output  quantities  is  the  same 
as  that  of  Table  III. 
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TABLE  II.  MAIN  PROGRAM  FORTRAN  LISTING 


«  M1RTRAN  OFFK 

CCAU  KAVfcPICK  S I MULA1 1  ON 

C-iNMON  '/SSAH/  |TnV,  Nil,  TNEXl .  Vi  I  n ;  s  TPHx  .  SI  23  A* .  SIim;:?2 
1  .tFTA.  SETA.CNI),  SNH,  1  H  A  X  »  NZ.I  NVC»n  ) ,  M  TLE  r  250  )  .HH  10.  NT  I  TIXV  > 
2.  I  F (.! F N *  1  HI  OTN*  N|.  0E.N2,  I  FG2N 
C’iMmOM  /SiiAM/  PFAn.OH  T.A»Tnr.TI«F 
CoPMllN  7SSAK2/  V  (25(1),  T  (2 ‘tni.r,  (250  ) 

C 1 1 M H 0 N  /TKAKFk/  CtlUM.TR.Nl 
j  .crrx.r.iFY 


c 

c 

c 


^ I  VAI  ENCF 
1  i V(  1  )» ALT 
V  <V<  5  ) »  D  A 
3  (  V<  9  ) ,  V  Y  K 
*  <V(13),WZ 
(V(17)»AZCMn 
t  (V(?l), PITCH 
7  <V(?5),TEAP 
h  (  V  (  2  V  > ,  R  X 
•J  t  V  <  33  ) »  ALPHA 
Foil  1  VAI  ENCF 
J  ( V ( 37  ) , V  YE 
(  V  (  A  l  ) ,  A  H 
3  tV<45),DWY 
*■  (V(A9),ni)Yr 
<■  <  V  (  5  3  • .  H  X  P 
7  ( V ( 5  /  ) , DWVP 
»>  (  V  <  61  ) ,  A  ZP 

X. ( V(67),RXm 
i  <•!)  I  VAI  ENCF 
1 (V(66).TOTHIS) 

1  ( V ( 72  ) »  AXF 

2  (V(76),At)SV 

3  »  V(B0  >,  X  ISCHF 
A  ( V ( 85  ) »  DF 

E  JIM  VAI  ENCF 
)  I  V<U(jJ.Whb__ 
2  ( V ( 1 21 ) • ETAD 
.)  (V(125).XLPSD 
A  ( V ( 129 ) , YAWP 
5  (  V(133).IIFXSI) 
Fi,l  I  VAI  ENCF 


.  (  V  (  2  ) ,  I)  A  C 
,(V(  6), DP 
,(V(1H).V7M 
.  (V(H  >,AXM 
,  (  V(  10  )  ,  AYC«P 
,  <  V<  22 ) , RXS 
, ( V<  26  ) ,  TF  A  Y 
.  (  V  (  3  0  ) ,  R  Y 
, <  V(34  ) ,  AI.PMAP 


, <V(  3  ) « |>PC 
. (V(  /  ) »  P  Y 
, (V(1 l ),WX 
, ( VC  151, AYM 
. ( V( 1 9  ) , YAW 
. ( V(2J).PYS 

.  <  vc  ?  / )  ,srr,« 

, ( V  <  A  t  ) , FP1 
,  (  V(35),AI.Pf 


).(V(  A  )  ,  i) YC 
>.(V<  fll.VYP 
),(V(  I  ?  )  ,  V,  Y 
•  »(V(16),A7fl 
> , ( V( 20  ) , MOl L 
) ,  (V(  v<\  )  ,  R7S 
)  ,  (  V(  28  )  ,  SAGA 
) , ( V( i? ) , EP? 

A  Y ) , ( Y (  «6>. VXf 


,  <Y<38).  VZF 
» ( V <  A? ) , ACP 
. (Y(A6),UW7 
» ( V ( 50 ) , TSM I SS 
, ( V(5A ) , WYP 
,  (V(6B),IIWZP 
, ( V(62 ) , VXP 
, CV(68),HYH 


, ( V( JV  ) .  U 
.  < V(A3), ACY 
.  (  V<  4  7  ) ,  PDA'; 
,  <  V  <  S  1  ) ,  V  S  H  i 
, ( V(59),w;p 
, ( V ( 59  ) , A  XP 
. ( V(A3), VYP 
. I V(69),WZH 


)  .  (  V  (  4  0  )  ,  v  M 
)  .  (  V<  *4  )  .  |,WX 

>,(vu0>.onpc 
SS),(V»'.P),/SMlSS 
)  ,  (  V  (  9  A  )  ,  ||  W  X  P 
).(V(MI),AYP 
) , ( V ( A A ) , V  7P 
) 


(  V  (  7  0  )  ,  F  H  J 
, (V<73), AYE 
,  (V(77),II0HNT 
. ( V<81 ) , XLUS 
, ( V(B6),DFXS) 


(  V  (  71  )  .  l:MK  ). 

,(V(74),AZF  )  , ( V( 75 ) . ARWH  ). 

,(V(7rt),YERlBT)>(V(/9),!;|0HAF), 
. <V(H2>. YLOS  )  #  (  V  ( 1.3  ) ,  /  L  P5  ), 


, ( V( 1 1H ) , TEYH 
,  ( V ( 1 22 ) , M XP 
,  (  V(  126  ) ,  YLClSlI 
,  (  V  ( 13(1 ) ,  ROL I  P 
,  (V(13A),DH). 


,  (  V(119),  IF.'D  ) ,  (  V<  l?lt  )  .  1 1  ►  0 
,  (V(123),WY.i  ) .  (  V  ( 1  ?A  ) ,  W7D 

.  (  V(1  2  7  ) ,  Zl.USh  ) ,  (  V(  126  ) ,  ANT 
.  (  V( 131 ),P I  TP  ) .  ( V  < 1 32 ) . nFO 
(  V(1  I9),nn  ).(V(tJ7),EtAR) 


) . 
) . 
> , 
). 


1  ( (•  ( 109  > .  XK1  ),  (C(llQ).XK?). 

2("(J '( 1 1  j )  niu? )',  *  (T(114>,  if  ). 

EQUIVALENCE 

1 (  i :  ( 1 1 6  > , U 3 S  > ,  ( C ( 1 17  )  •  K2  T  ) 

2  <  0 ( 1 59 ) »  UE V ( 1 ) ) 


(I’.dlll.PMI,  (  C(  1 1  2  )  •  1  Alii  > 
(  C  ( 11*5  )  ,  T  C  ) 

(C(  lift  1 ,1'UPP)  ,  (C<  156  )  ,PC  ) 


THIS  MODEL  has  the  mTSSTlE  FLYING  at  a  A5  peg.  moli  angle 

S')P222  =  0.0  _ 

Rf  AD  =  1.0 
6  OulIN  1  =  0.0 
C  A I  I  L°AD 

T  m  X  =  fB _ _  __  __  _ 

“  trx=TC 

ntLllsIC-TH 

B1 7  =  0  , 


CALL-.  "1  0 


CALL" 08 U 
CA  L  L  i' n  90 
F  A  I  L  o  1. 0  « 
f  A  l.  L  0  1  1  H 
C  A  (  L  !'  1  2  n 
TAI  L"1  HI 
CAI  l  1  AO 
C  Al.l  MSO 
r  a  i  i  1 6  ii 
FAIL  1  7  (• 
r/ififtn 
r ai  i  n  on 

CAI  L  0?n  li 
CAI  I  '21 0 
CAI  I  "  22 II 
FALL" 2 30 
r  A  L  L  ?  4 II 
I  AI  L *  °5 n 
FALL "260 
CAI  I.'  271' 
F  A I  I.  •'  2H  ll 
TAI  L  ”  2  9 1* 
CALL  030  0 
CAI  L  3111 


r  A  L  L  0  3  2 II 
CALL " 330 
f  A  I.  L  "  3  A  0 
C  A 1.  L  o  350 
TALI  0360 
FALLP370 
CAI.LO3B0 
CALLH3B1 
C  A  L  L  *’  3  ft  2 
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TABLE  II.  MAIN  PROGRAM  FORTRAN  LISTING  (CONTINUED) 


ci  =  o. 


)U1 


GO  TO  330 
T1=TIHF*0ELT 
C 1  *  A  l  N  T  (  T 1  /  T  C ) 
B1=AINT((T1*DELT1)/TC) 


IF  (Cl.EO.ClZ)  GO  TO  31 
N 1  =  1 


330  CONTINUE 

ANGLE  RESTORATION  BIAS 
IT  (TINE. GT. C(*>>  CALL 


IF (FLAG. 10.1. LCALLPRF 


50  CAILNAERO 
120  1=2 


CALLnSiSO 

CALL'i*70 

Call«*ih 

CALL 0590 


ALLntl 
CAL  L 11600 
C  A  L  L  n  6  ?  (I 


TABLE  IL  MAIN  PROGRAM  FORTRAN  LISTING  (CONCLUDED) 


CALL  MPILOT 
130  I =  3 

- CALL  mFLTP - 

_  HO _ I  *4 _ _ _ 

CALL  R45R(WX, WY. WZ.MXP.MYP.WZP) 

_  CALL  EULANflIWXP.WYP.WZP, YAM,  ROLL. PITCH) _ _ 

CALL  R45A<DWX.DWV,OWZ,DWXP,DWYP.OW*#> 

_ CALL  R45R(AXH,AYH,AZ8,AXP,AYP,AZP) _ 

CALL  R45R(VXN, VYH.VZH. VXP,  VYP, VZP) 

_  CALL  EULTRNU.l,  VXE.VYE.VZE,  VXP.  VYP.  VZP.  YAW,  ROLL.  PITCH  > 
CALL  HSEEK 
150  1*5 

CALL  SfkTR  <0,-l,RSH.RYH»R7M7ftY5.HY57irZ5» SEtTA ,'SaOa > 

160  1*6 _ _ _ _ _ _  _  _ 

CALL  EULTrS  (O.iVrX.Ay.ALT.RXM.RVM.RZM.  Y  Ah,  r6u.  PITCH) 
170  1=7 

CALL  EULTRWIO.l,  AXE » A  YE . AZE.  A XP, A YP.  AZpTY A W, ROLL  .PITCH) 
RSO  =  V(?9)«*2*V<30  )«»2* V( 1 ) 

XLOSs ( V ( 30 ) • V ( 38  ) - V ( 1 ) • V ( 37 ) > /RSQ 

_ YLOS»(V<l)»V<36)-V<29)«V<38) )/RSO _ 

ZLOS*( V(29)*V(37  )-V(30 )«V(36) )/RS0 

XLOSD=XLOS«RTOn _ 

YLOSU=YLOS»RTOn 

ZLOSO=ZLOS»RTOD 

WXD=MX«RTOO 

_ wyd=wy*rt  on _ _ _ _ 

wzn=wz«RTon 

TE YD=TEAV«RTnD  _ 

TcPD*  TEAP*RTOD 

_ ilND  =  SEOA«RTOO _ _ _  _ 

ETAD*SAGA*RTOD 

_ YAWD  =  YAWRTOf) _ 

RnLLD=POLL*RTOP 

PI  TD«PI  TCH»RTOP _ _ _  _ _ _ 

l)LD*0E*RTon 

nEXSO  =  DEXS*RTOD _ _  _  _ 

"  ANT*S0RH  (AYM)**2*(AZM)««2) 

_ V<65  )=S1?345 _ _ 

C  I f  (T1HE.LT.THX)  00  TO  60 

c _ mxxtajLLTJL _  _  . 

C  60  CAI  L  TIES  T( T0X 1 

If  ( T  lhE.LT.JTCX)  10  10  7Q  _ 

TCX=TCX*TC 

70  CALL  TTEST(TCX) _ 

CALL  PPINT$(-V(2?) > 

_  CALL..  BEAU _ - 

IF  ( I  END  >  20.1000,1000 

10 DO  WRITE  (6.1001)  COUNT _ _ _ _ 

1  0  tl  l'  FORMAT  <  IH1 » 28H  TOTAL  NUPBFR  OP  I  TER  A  T  I  ONS*  .  1  PE15. 7  ) 

_ Qii  TO  6 _ 

END 


~  CA  L  O  6  3  0 

CALLo640 

CALL  0  650 
CALL0660 
CALL0670 
CALL  1/680 
CALLH690 
CALL0700 
CALLO710 

CALL072U 

7!TiX«Tll~ 

CALL6740 

CAI.LnY70 

CALLr.78U 

CALL|,790 

fALLonnn 


.CAlL'1810 

CALL6fl?X 

CALL'1821 

CALL 


CALL "840 
CALL!)  850 
CALLO860 
CALL ''920 
CALL0Q30 
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(U)  A  listing  of  this  subroutine  is  shown  in  Table  VI,  and  a  flow  chart 
appears  in  Figure  11. 

(U)  Several  points  must  be  clarified  regarding  two  of  the  inputs  to  this 
subroutine.  The  effective  tracker  time  constant  T(  13)  is  used  to  establish 
initial  tracking  error  angles  only;  it  is  not  used  thereafter.  The  assump¬ 
tion  is  made  that,  at  time  of  launch,  the  tracking  loop  has  achieved  steady 
state,  so  that  tracking  error  is  proportional  to  the  product  of  line-of-sight 
rate  and  tracker  time  constant.  This  initial  error  may  be  eliminated  by 
setting  this  input  to  zero.  In  this  case,  the  performance  of  the  simulation 
would  be  otherwise  unchanged. 


MISS 

INERTIAL  distance 


INERTIAL  AUTOPILOT 


AUTOPILOT  SEEKER 


Figure  4.  Piograms  Showing  Euler  Angle 
Relationships  Between  Coordinate  Sets 
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ya  AND  za  ALIGNED  WITH  AUTOPILOT  AXES 
y  AND  z  ALIGNED  WITH  MISSILE  BODY 


MISSILE  BOW  AND  AUTOPILOT  AXES 


(ALIGNED  WITH  X  AXIS) 


\ 

4 


MANEUVER  AXES 


Figure  5.  Simulation  Coordinate  System 
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15 

(The  reverse  of  this  page  is  blank) 


Figure  6.  Call  Program  Flow  Chart 


ALT 


Figure  7.  Call  Program 
Flow  Diagram 


(The  reverse  of  this  page  is  blank.  ) 


TABLE  III.  CALL  PROGRAM  GLOSSARY,  V  ARRAY 


Name 

Quantity 

Units 

Coordinate 

System 

V(l) 

hm,  Missile  altitude  above  ground 

ft 

Inertial 

V(2) 

5ac,  Aileron  deflection  command 

deg 

Missile 

V(3) 

5pC,  Pitch  deflection  command 

deg 

Missile 

V(4) 

6yc*  Yaw  deflection  command 

deg 

Missile 

V(5) 

5a,  Aileron  deflection 

deg 

Missile 

V(6) 

5p,  Pitch  deflection 

deg 

Missile 

V(7) 

5y,  Yaw  deflection 

deg 

Missile 

V(8) 

Vx,  Missile  velocity  X-axis 

ft/ sec 

Missile 

V(9) 

Vy,  Missile  velocity  Y-axis 

ft/sec 

Missile 

V(10) 

Vz,  Missile  velocity  Z-axis 

ft/sec 

Missile 

V(ll) 

u>x.  Angular  velocity 

rad/ sec 

Missile 

V(12) 

u>y,  Angular  velocity 

rad/ sec 

Missile 

V(  13) 

<*>z,  Angular  velocity 

rad/sec 

Missile 

V(  14) 

Ax,  Propulsion  and  aerodynamic 
acceleration 

g 

Missile 

V(1 5) 

Ay,  Propulsion  and  aerodynamic 
acceleration 

g 

Missile 

V(  16) 

Az,  Propulsion  and  aerodynamic 
accel  eration 

g 

Missile 

V(17) 

Azc,  Elevation  maneuver  command 

g 

Autopilot 

V(  1 8) 

A  ,  Azimuth  maneuver  command 

y  c 

g 

Autopilot 

V(  19) 

^  ,  Euler  yaw  angle 

rad 

V(20) 

<(>,  Euler  roll  angle 

rad 

V(2 1 ) 

0,  Euler  pitch  angle 

rad 

V(22) 

Rx,  Seeker  boresight  range 

ft 

Seeker 

V  ( 23 ) 

Ry,  Seeker  lateral  range 

ft 

Seeker 

V(24) 

Rz,  Seeker  normal  range 

ft 

Seeker 

V(25) 

«z.  Tracking  error  angle,  pitch 

rad 

Seeker 

V(26) 

€y,  Tracking  error  angle,  yaw 

rad 

Seeker 

V(27) 

v,  Seeker  elevation  gimbal  angle 

rad 

V(  28 ) 

Tl,  Seeker  azimuth  gimbal  angle 

rad 
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TABLE  III.  CALL  PROGRAM  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 


V(29) 

V(  30) 

V(31 ) 
V(32) 
V(33) 
V(34) 
V(35) 
V(36) 

V(37) 

V(38 

V(39: 

V  (40 
V(41 

V  (42 
V(43 

V  (44 
V(45 
V(46 
V(47 
V(48 

V  (49 
V(  50 
V(51 
V(52 
V(53 
V(54 
V(55 


Quantity 


R^,  Horizontal  longitudinal  range 
component 

Rj,  Horizontal  lateral  range 
component 

€gZ,  Gate  error  angle,  pitch 

«gy,  Gate  error  angle,  yaw 

a,  Total  miss  angle  of  attack 

a  ,  Missile  pitch  angle  of  attack 

ciy,  Missile  yaw  angle  of  attack 

Vj,  Horizontal  longitudinal  velocity 
component 

Vj,  Horizontal  lateral  velocity 
component 

V^,  Vertical  velocity  component 

q,  Dynamic  pressure 

Total  missile  velocity 

Missile  Mach  number 

aCp,  Acceleration  command  pitch 

aCy,  Acceleration  command  yaw 

0)  f  1 

x  Scalar  components  of  mis-  i 

<iy  sile  angular  acceleration 

.  in  missile  axes 

(a) 

Z 

5ac*  Aileron  command  rate 
5  ,  Elevator  command  rate 

5yC,  Rudder  command  rate 
Closest  approach  at  end  of  flight 
Range  component  in  Y  seeker  axis 
Range  component  in  Z  seeker  axis 
w  , 

,  I  Missile  body  rates  in 
y  j  autopilot  axes 


Units 

Coordinate 

System 

ft 

Inertial 

ft 

Inertial 

rad 

Seeker 

rad 

Seeker 

deg 

Missile 

deg 

Missile 

deg 

Missile 

ft/ sec 

Inertial 

ft/  sec 

Inertial 

ft/ sec 

Inertial 

lb/ft2 

ft/ sec 

g 

Autopilot 

g 

Autopilot 

rad/ sec2 

Missile 

deg/ sec 

Missile 

deg/sec 

Missile 

deg/ sec 

Missile 

ft 

ft 

Seeker 

ft 

Seeker 

rad/sec 

P  utopilot 
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TABLE  III.  CALL  PROGRAM  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

Coordinate 

System 

H 

0)< 

X 

0)' 

y 

0)' 

z 

Scalar  components  of  mis¬ 
sile  angular  acceleration 
in  autopilot  axes 

rad/sec^ 

Autopilot 

V(59) 

V(60) 

V(6 1 ) 

A' 

X 

Ay 

A1 

z 

Propulsive  and  aerodynamic 
acceleration  components  in 
autopilot  axes 

g 

Autopilot 

V(62) 

V(63) 

V(64) 

V' 

X 

V' 

y 

Vz 

Missile  velocity  components 
in  autopilot  axes 

ft/  sec 

Autopilot 

V(65) 

Special  test  variable  -  used  as 
system  diagnostic 

V(66) 

Total  miss  distance 

ft 

Miss  Distance 

V(67) 

x  component  of  range 

ft 

Autopilot 

V(68) 

y  component  of  range 

ft 

Autopilot 

V(69) 

z  component  of  range 

ft 

Autopilot 

V(70) 

y  component  of  miss 

ft 

Miss  Distance 

V  ( 7 1 ) 

z  component  of  miss 

ft 

Miss  Distance 

V(72) 

x  component  of  acceleration 

g 

Inertial 

V(73) 

y  component  of  acceleration 

g 

Inertial 

V(74) 

z  component  of  acceleration 

g 

Inertial 

V(75) 

y  component  of  acceleration  at 
blind  range 

g 

Miss  Distance 

V(76) 

z  component  of  acceleration  at 
blind  range 

g 

Miss  Distance 

V(77) 

Blind 

time  in  yaw  channel 

sec 

V(78) 

Blind 

time  in  pitch  channel 

sec 

V(79) 

Final  line  of  sight  angle  (vertical) 

rad 

Inertial 

V(80) 

Final  heading  angle  (horizontal) 

rad 

Inertial 

V(81) 

x  component,  LOS  rate 

rad/sec 

Inertial 

V(82) 

y  component,  LOS  rate 

rad/sec 

Inertial 

V(83) 

z  component,  LOS  rate 

rad/sec 

Inertial 

V(84) 

A,  Guidance  gain 
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TABLE  HI.  CALL  PROGRAM  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 


V(87) 

V(88) 

V(90) 

V(91) 

V(92) 

V(9j) 

V(94) 

V(95) 

V(96) 

V(97) 

V(98) 

V(99) 

V(100) 
V(101) 
V(102) 
V(103) 
V(104) 
V(105) 
V(106) 
V  ( 1  07 ) 
V(  1 08) 
V( 1 09) 
V(  1 10) 


DEXS 


FFE 

DFE 


FEXS 

DFEXS 


Quantity 


Total  yaw  precession  rate 

Total  pitch  precession  rate 

Yaw  gyro  inertial  angle 

Yaw  look  angle  (indicated) 

Forcing  function  cross- 
coupled  equation  1 

Derivative  forcing  function 
c ross -coupled  equation  1 

Forcing  function  cross- 
coupled  equation  2 

Derivative  forcing  function 
cross-coupled  equation  2 

Integral  forcing  function 
cross-coupled  equation  1 

Integral  forcing  function 
cross-coupled  equation  2 

Forcing  function  yaw  axis 

Derivative  forcing  function 
yaw  axis 

Forcing  function  pitch  axis 

Derivative  forcing  function 
pitch  axis 


Units 


Coordinate 

System 


NOT  USED 
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TABLE  HI.  CALL  PROGRAM  GLOSSARY,  V  ARRAY  (CONCLUDED) 


Name 


V(1 13) 

V(1 14) 

V(1 1 5) 

V(1 16) 

V(1 17) 
V(1 18) 
V(’  19) 


V(122) 
V(123) 
V(124) 
V( 125) 
V( 126) 
V(127) 
V(128) 
V(129) 
V(  130) 
V(131) 
V(  132) 
V(133) 


Quantity 


Suml  -  Tracker  sampler  bias 

TEAYD  Tracker  error  yaw  • 
RKAMG 

TEAPD  Tracker  error  pitch  • 
RKAMG 

TEAYS  -  Tracker  Z0H  output 
Signal,  Yaw 

TEAPS  -  Tracker  Z0H  output 
Signal,  Pitch 

VSYP  -  Tracker  output  signal 
Pitch 

VSPP  -  Tracker  output  signal  yaw 
TEYD  -  Tracking  error  -  yaw 
TEPD  -  Tracking  error  -  pitch 
UND  Seeker  elevation 
ETAD  Seeker  azimuth 
WXD 

WYD  Missile  angular  velocity 

WZD 

XLOSD 

YLOSD  LOS  Rate,  Inertial 


Units 


Coordinate 

System 


ZLOSD 

ANT  (New) 

i|i  Yaw  D 

4>Roll  D  Error  Angle 
0  Pitch 

DED  Total  precession  rate,  yaw 

DEXSD  Total  precession  rate, 
pitch 


deg/ sec 

deg/sec 

deg/sec 

deg 

deg 

deg 

deg 


deg/sec 


deg/ sec 
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TABLE  IV.  INPUT  TO  SUBROUTINE  SETIC 


Input 

Location 

Quantity 

Units 

T(  1 ) 

R,  Total  range  to  target 

ft 

Vq,  Launch  velocity 

ft/sec 

tr0,  Line  of  sight  angle 

deg 

T(4) 

YEV’ 

Heading  error,  vertical 

deg 

T(5) 

yEH’ 

Heading  error,  horizontal 

deg 

T(6) 

4*^,  Aircraft  roll  angle 

deg 

T(7) 

AaL* 

Aircraft  normal  acceleration 

g 

T(8) 

AaY ' 

Aircraft  lateral  acceleration 

g 

T(9) 

apQ,  Aircraft  angle  of  attack,  trim 

deg 

T(  1 0) 

Angle  of  attack,  gain 

deg/g 

T(  1 1 ) 

4>£,  Missile  mounting  angle,  roll 

deg 

1 

0£,  Missile  mounting  angle,  pitch 

deg 

t  ,  Effective  tracker  time  constant  used 
lo  calculate  initial  tracking  error  angle 

sec 

T(14) 

o 

X 

rad/sec 

T(  1  5) 

wy 

|  Missile  body  rates  in  autopilot  axes 

rad/sec 

T(16) 

rad/sec 

T(  1 7) 

RBH» 

Blind  range,  horizontal 

ft 

T(  1 8) 

rbv* 

Blind  range,  vertical 

ft 

T(  19) 

Steering  bias,  pitch 

g 

T(20) 

Steering  bias,  yaw 

g 

T(2 1 ) 

Roll  rate  bias 

rad/sec 
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TABLE  V.  SUBROUTINE  SETIC  OUTPUT 


Output 

Variable 

Quantity 

Units 

Coordinate 

System 

V(l) 

h,  Missile  altitude 

ft 

Inertial 

V(8) 

vx 

V(9) 

V 

Missile  velocity 

ft/ sec 

Missile 

y 

components 

V(10) 

V 

z 

V(ll) 

wx 

V(  1 2) 

Missile  angular  velocity 

rad/ sec 

Missile 

y 

components 

V(  1 3) 

uz 

V(19) 

* 

Missile  Euler  angles  in 

V(20) 

$ 

yaw,  roll  and  pitch 

rad 

V(2 1 ) 

e 

sequences 

V(22) 

^xs 

V(23) 

V 

Target  range  components 

ft 

Seeke  r 

V(24) 

V(27) 

v,  Seeker  elevation  gimbal  angle 

rad 

V(28) 

9,  Seeker  azimuth  gimbal  angle 

rad 

V(29) 

Ri,  Horizontal  range  component 

ft 

Inertial 

to  target 

V(67) 

Rxm 

V(68) 

Rym 

Target  range  components 

ft 

Autopilot 

V(69) 

Rzm 
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x,  y,  z  INERTIAL  COORDINATE  AXES 
V,,  V.,  Vk  AIRCRAFT  VELOCITY  COMPONENTS 

THE  INERTIAL  AXES  ARE  SELECTED  SO  THAT  THE  INITIAL 
RANGE  VECTOR  IS  CONTAINED  IN  THE  x-z  PLANE 

Figure  8.  Launch  Geometry  in  Inertial  Coordinates 
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u,  v,  w  AIRCRAFT  COORDINATE  AXES 
Vu#  Vy,  Vw  AIRCRAFT  VELOCITY  COMPONENTS 


ANGLE  OF  ATTACK  COMPONENTS  ARE  CALCULATED 

WITHIN  SETIC  AS  FOLLOWS: 

% 


Figure  9.  Launch  Geometry  in  Aircraft  Coordinates 
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INERTIAL  AIRCRAFT 

AXES  AXES 


AXES  AXES 

♦a  AND  0a  ARE  COMPUTED  WITHIN  SETIC. 

ALL  OTHER  ANGLES  ARE  PROVIDED  AS  INPUTS. 


Figure  10.  Euler  Angle  Relations  Between  Coordinate  Axis  Sets 
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TABLE  VL 


SUBROUTINE  SETIC  FORTRAN  LISTING 


t  INMTNAH  I'FCK 

c^f  1 1 c  si- 1  it.i  i  ml  cn  omoNS  i 

snpwoiiTiNt  seiif; 

COMMON  /SSAM2/  v  (250). T  i250),C  (250) 

F<l|i  |  VALENCE 

1  ( T  (  l). RANGE  >»(T(  2 )  ,  VEL  ).(T(  3), SIGMA  >.(!<  4),HFV  ), 

2  <T(  5  )  #  h£  H  )  <  (  I  (  6  I , ACROLL )  *  (  T  (  7), ACCEL  ).(!(  8  ) . ACLFL Y ) # 

3  t t <  q ) , alpo  ) » ( r < i n ) , oalda  i.min.Phii  ). < Mi2), ihltal). 

A  I  t  ( 13  ) ,  1 AUA  ).(T(14),HX  ).(T(15),MY  I. (1(16), HZ  > 

oata  kuw/57 .vi'tm/ 

0(0)31(17) 

C  (  3  )  =  T  ( 1 8  ) 

(5AMVs(SIGMA-llfeV)/»T00 

TAMi,V*TAN(GAMV) 

HMO<AII  =  MFH/«Tnn 
T  AMDs  T  AN(  ilEiiHAr)  > 

C  VI  MICITV  COMi'ONl  N t S  IN  EARTH  AXLS 

VI=VFL/Si)NT(l.n*TANGv'»»2*rAKH««?  > 

V  1  =  V  I  *  TANll 
VK  *  V  I  • 1 A MOV 

A  I  PI'AI’=(AI  P(i*UALnA«ACCFL)/filOU  * 

A I  ri'Aysl'ALPA#ArCtLY/PTOII 
T.M  APsTAMl  AI.PHAI’ ) 

TAMA YsTAIM  Al PHAY  ) 

C  VriOCITY  POMPOM  NTS  IN  a/C  AXES 

VH= VEl/SftMTI  t.n*TANAP»«2*IAMY«»2) 

V  V  *  V IJ  •  T  A  N  A  v 
VSsVIU  TANAP 
A'-MI  3ACR0I.L/R10M 
S'M.lsSiraAPHI  ) 

CPH  I  sCOS( APIH  ) 

C  FSTAHLIS'I  A/C  FILLS  ANGLES 
cnr'tsvv»SPni*vM»i:PMi 
CHM?*SORT(  VU»V|I«C<|N1»C0N1  ) 

TMF1  AA  =  AH5N(CON)/i:0M?)-AKSN(  VK/M0N2  ) 

C'IM.IsSIlpr  (  V  I  •  V  I  *  V  i»V  J) 

PS  I  A*AUSN(  VJ/00r.3)-Ai»SN(  (  VV  »CPtH  -  VW«  SP H  I  )/CuNJ) 

$PSI*SIN(PSIA> 

cpsi3cos(psia) 

STPt  =  SI'((  I  HE  T  A  A  ) 

CTl't  scns<  I  HE  r  A  A  )• 

T"f  f  I. IIs  THE  T  AL  /R  I  Of) 

ST.H  *  S  I  N  <  IHETLR) 

Cll'l  =  COS(  IHFTLR) 

pmi  RsphIl/ntud 
SMhi  SSIN(PHII.K) 

CPMI =CnS<PHlLK) 

A133-STHf.:*C1HL*CTrtE«SPMI»SPFl.*S)HI.-CTFE«CPH|*CPilL»SlHL 
A  ?3  =  CTHE»SPlM«Cprt|.*C!ME»CPKl»SlMI 

A?1=  (rPSI»STHE*SPHI-SPSI»C>H|  >«C8HI.  ♦  (  CPS  I  •  S  ( ML  »CPm  I  ♦  SI  S I  tSPHI  )• 
1  SPhL 

Cr>M3S0Rl(l.n-A?J»A2J) 

C  ESTABLISH  MISSILE  EULL-R  ANGLES 
P*i|=AHSN(  A23) 

7 nr  i  a  =  arsn(-ai3/com  i 
PSirRsARSN(-A21/C0NA ) 

SSM;*SIN(SIGMA/hfi)0) 

CS I  H*COS(  SI  GMA/H  TUI) ) 

C  RAMi F  COMPONENTS  IN  EARTH  AXES 

R|rlANGE»CSIO 


SFTIU010 

SLTI0P2U 

SET  10050 
Sf  TI0060 
SET I0H7C 
SFTIOOBI) 
SIT  I  0090 
SLTIOIOU 
SET  I  01 10 
SITI012U 

sftiouo 
c»  r 1 0140 
SFTI015U 
SET  I  0160 
SFTI017U 
St  TIOlflU 
SKTI019U 
ST  TI0200 
SFTI021U 
St T 10220 
Sf  TIU23U 
SFTI0240 
SET  I  0250 
SE  T I  0 26 0 
SETI0270 
SFT1U28U 
SI- T  I  02911 
SF  Tl 0300 
St  T I  0 3 1 0 
SFTI0320. 
SFTI033U 
SF  T I  034  U 
SF  ri035U 
St  TI036II 
SFTI037U 
St  TI0380 
SE  T  I  0390 
S  E  T  I  0  4  0  0 
SI  T  I  0410 
SLTM420 
SF  riC430 
SFTI044U 
SLTIU45U 
SF  T I046U 
St  T I  0  4  70 

se  r i o  4  a  n 

StT I 049U 
Sr  T  I  05 0  U 
SI  110510 
SE  r  10520 
SE  T  I  0 53 0 
SE  NU540 
SF  T I  0550 
SF  T I  0560 
SFTI057U 
SE  T I  0580 
SF II 0590 
SE  T I  06 0  0 
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TABLE  VL  SUBROUTINE  SETIC  FORTRAN  LISTING  (CONCLUDED) 


RKsRANGE*SSIG 

C  ESTABLISH  TRACKING  ERROR  VECTOR 

CONh*VK»CSIG-V|»SS|G 
EI*1AUA*CON6»SSlG 
FJ=  TAUA» VJ 
EK*-TAUA*C0N6»CS|B 
PI  SR  I -El 
PKsRK-EK 

C  ESTABLISH  SEEKER  OIMflAL  ANGLES 

CAI  L  EULTRNI  l.-l.EI  .EJ.PK.HX.RY.R/.ESIPR.PHI .  THEU) 

C  VFIOCITY  COMPONENTS  IN  MISSILE  AXES 

CAI  L  EULTRN\-1»-1»VI«VJ,VK«VX.VY,v2»PSIPR#PmI,TmETA) 

C  RANGE  COMPONENTS  IN  MISSILE  AXES 

CAiL  EUL'i<N<-l,-t.RI .U.U.RK.RXM.RYM.RZM.PSIPR.PMl.ThETA) 
C  VFIOCITY  COMPONENTS  IN  AOIOI'ILOI  AXES 

CALL  N45F(VX,VY.VZ.VXM,VYM.VZM 
UN=ATAN?(-RZ.RX) 

ETA  =  ATANi!<RY,SORriRX«MX*RZ»HZn 
C  RANGE  COMPONENTS  IN  StEKER  AXES 

CALI  SEKTR(l,l,MM*NTM,RZH.RXb.NYS.RZS.UN.EIA> 

V ( 8 ) - VXM 

V  1 9  )  *  V YM 
V(10 )sVZM 
V(19)=PSIPR 

V  (20 ) *PH I 
V(?l JsTHETA 
V<?7)*UN 
V(28)*ETA 

V ( 1  >  =  RK 
V (29 ) =R I 
V (67  > sRXM 
V(68)sRYh 
V ( 69 ) =HZM ' 

V ( 22  )  »RXS 
V (23 ) sRYS 
V( ? A ) sRZS 

CALI.  R45F(WX»NV»MZ»V(11)*V(12)#V(13)> 

RETURN 

END 


SFTI0610 
SET  I  0620 
SE  T I  0630 
SFTIU640 
SE  T I  0650 
Sb  T 1 0660 
SET  10  670 
SE  T  I  0680 
SE  T I  0690 
S  L  T  1  0  7  0  0 
SET  10710 
SET  I  0720 
SET  10730 
Si TID740 
S  F  T  I  0  7  5  0 
SFTI0760 
SFT  I  U 7 70 
SETI0780 
SE  T I  0790 
set  i  080.0 

SET  I  0810 
SF  T  I  0820 
SE  T I  0830 
SE  T  I  0840 
SE  T I  0850 
SET  I  0860 
SrTI0870 
SFTIOflOO 
SFTI0890 
SE  T  I  0 90  0 
SLTI0910 
SF  T  I  0920 
SET  I0R30 
SE  T  I  0940 
ST  f  10950 
SE  T I  0960 
SET  10970 
SE  T I  0980 
SE  T I  0990 


4 


t 

(U)  Also,  the  45-degree  rotation  of  the  missile  axes  relative  to  the 
mounting  hooks  is  not  contained  in  the  input  T(ll),  This  angle  will 

nominally  be  approximately  zero  or  ±90  degrees  depending  upon  whether 
the  missile  is  mounted  on  the  bottom  or  sides  of  the  pylon. 

2.  3.  4  Universal  Seeker  Subroutine 

(U)  Input  to  the  Universal  Seeker  subroutine  are  the  components  of  the 
missile  linear  velocity  and  acceleration  components  as  well  as  the  angular 
velocity  and  acceleration  components,  all  in  autopilot  axes.  This  subroutine 
performs  the  coordinate  transformations  of  these  quantities  to  seeker  axes 
and  performs  the  appropriate  integrations  to  yield  seeker  range  components 
from  which  tracking  error  angles  in  the  pitch  and  yaw  planes  are  computed. 
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DATA  *100/57.29578/ 
CI'HTUT) 

CO*-T<18> 

GAMV*<SIGMA-HEV)/VT00 

TANCV-TAWGAMV) 

HEHRAD-HEH/RT00 

TANH»TAIKH£HRA0) 

- V - 


VELOCITY  COMPONENTS  IN  EARTH  AXES 
V1=VEL/SQRTU  .0+TANCV**2+TANH**2> 
VJ»VI*TANH 
VK*VI*TANCV 

ALPHA  P=(ALPO+OALBA*ACCEL)/RTOD 
ALPHAY*DALDA*ACCELY/RTOD 
TANAP-TANIALPHAP) 
TANAY»TAPKALPHAY) 

*  - 

VELOCITY  COMPONENTS  IN  A/C  AXES 

VU’VEL/SQRTIl  .0+TANAP**2+TANAV*»2) 

W*VU*TANAY 

VW=VU*TANAP 

APHI«ACROLL/RTOO 

SPHI'SIMAPHII 

CPHI^COS(APHI)  _ 


ESTABLISH  A/C  EULER  ANCLES 

C0N1«W*SPHI+VW*CPHI 

CON2*SQRT(VU*VU+CON1*CON1) 

thetaa<aRsikconia:on2>-arsn(vka;on2i 

C0N3«SQRT(VI*VI+VJ*VJ> 

PSIAtARSN(VJ/C0N3>-ARSN((W*CPHI-VW*S  PHD/CONS) 

SPSI*SHKPSIA> 

CPSI-COS(PSIA) 

STHE»SIN(THETAA) 

CTHE«COS(THETAA) 

THETLR*THETAL/RTOC 

STHL*SIN(THEfLR) 

CTHL-COSITHETLR) 

PHILR*PHIL/RTOD 

SPHL*SIMPHILR) 

CPHL*COS<PHILR) 

A1J«STHE‘CTML+CTHE*SPHI*SPHL*STHL-CTHE*CPHI*CPHL*STHL 

A23<THE*SPHI*CPHL+CTHE*CPHI*SPHL 

A21a<CPSI*STHE*SPHI-SPSPCPHIPCPHL+<CPSI*STHE*CPHI+SPSI*SPHIP 

SPHL 

C0H4»SQRTU  ,0-A23*A23l 

- -  r-  - - 

ESTABLISH  MISSILE  EULER  ANCLES 
PHI=ARSIKA23> 

THETA=ARSN(-A13/C0N4) 

PSIPR=ARSN(-A21/C0N4) 

SSIG=SIN(SIGMA/RTOD) 

CSIG=COS(SIGMA/RTOO) 

- 1  _ 


Figure  1 1.  Subroutine  Setic  Flow  Chart 
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(U)  The  Seeker  subroutine  has  been  modified  for  the  CAS  system  simula¬ 
tion  to  accommodate  two  other  subroutines  which  simulate  tracker  and  gyro 
dynamics.  The  respective  subroutines  are  Track  and  Gyro;  their  addition 
facilitates  the  replacement  of  different  trackers  and  gyro  models  into  the 
simulation. 

(U)  The  Track  subroutine  accepts  error  signals  from  the  Universal 
Seeker  subroutine  and  simulates  the  tracker  dynamics.  Its  output  signals 
are  the  missile  acceleration  commands  to  the  autopilot  and  precession  rate 
signal  to  the  gyro  torquer. 

(U)  The  Gyro  subroutine  accepts  the  inputs  generated  by  the  Track  sub¬ 
routine  and  simulates  the  gyro  dynamics  including  drift.  Its  outputs  are 
the  gimbai  angles  and  gimbal  rates. 

(U)  A  FORTRAN  listing  of  the  Universal  Seeker  Track  and  Gyro  sub¬ 
routines  are  shown  in  Tables  VII,  VIII,  and  IX,  respectively.  The  subroutine 
block  diagram  and  flow  charts  appear  in  Figures  12,  13,  14  and  15,  respec¬ 
tively.  Tables  X,  XI,  and  XII  contain  a  glossary  of  terms. 

2.  3.  5  Aimpoint  Wander  Subroutine 

(U)  The  Aimpoint  Wander  subroutine,  GWAND,  is  called  out  within  the 
Seeker  subroutine  and  either  simulates  the  apparent  target  motion  caused  by 
the  wandering  of  the  seekei  airpoint  or  the  actual  motion  of  an  evasive  target. 
In  both  cases,  the  subroutine  input  is  the  boresight  range  Rxs. 

(U)  The  equations  implemented  in  each  of  these  options  are  shown  in 
Tables  XIII  and  XIV,  respectively.  If  neither  option  is  desired,  setting  the 
parameter  C(103)  to  zero  will  cause  the  entire  subroutine  to  be  bypassed. 

(U)  A  FORTRAN  listing  of  this  subroutine  appears  in  Table  XV  and  its 
flow  chart  is  shown  in  Figure  16. 

2.  3.  6  Angle  Restoration  Bias  (ARB)  Subroutine 

(U)  The  Angle  Restoration  Bias  subroutine  serves  to  implement  the  guid¬ 
ance  law  incorporated  in  the  CAS  missile.  It  operates  on  the  acceleration 
commands  from  the  Seeker  subroutine  to  provide  steering  commands  for  the 
Autopilot  subroutine.  A  FORTRAN  listing  of  the  subroutine  appears  in 
Table  XVI.  A  block  diagram  and  flow  chart  of  the  subroutine  appear  in 
Figures  17  and  18,  respectively. 

2.  3.  7  Blind  Range  Filter  (BRF)  Subroutine 

(U)  This  subroutine  simply  provides  the  filtering  for  commands  to  the 
autopilot  subroutine  when  blind  range  is  reached.  The  FORTRAN  listing 
for  this  subroutine  appears  in  Table  XVII.  The  subroutine  block  diagram 
and  flow  chart  appear  in  Figures  19  and  20. 
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TABLE  VII.  UNIVERSAL  SEEKER  FORTRAN  LISTING 


1  '  ffiRfRAN  LlTOU.DECK  "  '  . .  . 

CUSEfc  UNI VFRSAL  SEEKER 

- SHBWOTTBFTfsm - 

COMMON  /SSAMl/  RE*B^ DEL T.4IIT0T. TINE 
CiiNHON  /S5AE12/'  V  ( 290  ) .  f  I25#>.C  T?SoT 
CONN0N  /tRAKFR/  COUNT. TR.N1.OFFX.OFFY 
CUMHON  /TRAKZ  /  TEAPO .  TEA TO. DTFAP.OTEAY, WYSC.HZSC 
C  oM  H ON/OTR/ VST ,  V S P.WXS.WY S .  W Z_S 
"  F'W  I  VAI  ENCF*  -  .  -  - 

1  ( VI 22  > . RXS  ) > I  V ( 23  ) . RTS  >.(V(?4),RZS  ) . ( V I 25 ) . TEAP  >. 

2  (  VI26  ) .  TEXT  > .  <  V  <  1  7  )  .  AZCNII  > ,  (  VII  »  I ,  A  YCmD  >.(V(2  7). UN  ). 

3  (V(?8).ETA  I.  (V(Jl).EPl).  <V(3?),EP?> 

Poll  I  VAI  ENCF 

I  (V(SJ).NX  ).(V(*>4).WY  ).<V(Sl>>.h2  ),(V(*.A),PUX  ). 

t  IV(57>.niv*  I.IVIMI.OHZ  ).7v(RVf.AX  ).(V(A0).AV  >, 

3  I  V  (  A  1  ) »  AZ  )  >  (  V  (  62  )  .  V  X  )  >  (V(AJ).VY  ).(V(»,4),V/  ) 

F  III  |  VAI  ENCF 

|'<v<ioa).vs?>,  (vuot>),vsn,  ivmn>.wNi>> 

Roll  I  VAI  ENCF 

II  VC66  )  .  TOTNIS).  I  V  (  7  0  >,tM.I  ).(VI71),EHK  ) 

2.  (  VISU  ),  TSHSS).  C  Vt*>l  ),YN|SS  ) .  (  V  (  52 ) ,  ZH  I  SS  ) 

■  F  JUIVAI  ENCF 

i  n  il.!),  lAK'inr ),  i c 1 1 a > . r.K  ).  <i:i  15>. iai  ap>,  icu ai.mhfgal  > . 

2(1(17). Cl  >.(C(18).C2  ).  (  U  (  19  ) , C3  ), (1,(20). CA  ) , ( r ( 2 1 ) , CS  ), 

.'( I.  (  22  )  ,  C6  ).  <  C  (  2  3  ) .  (.  7  ),<i;c>4).CA  >.(.;(25).L9  ) .  ( U  ( 2  A  > .  C  Id  ) . 

Alt.  (27  )  ,dl  ).  (C(?«)»i;t2).  (C(29).C13).  (C(30  ).Ct  A  ),  (Clol  )  ,C11). 

Hi  (32  )  .Cl  A) ,  (U( J3).C1  7  ),  (C(.)A  >.C1  D.  ((.(  3*>  ),C1  » >»  (C(3A  )  ,C2  i ) , 

MUt  37  )  ,C21  1 ,  (Cl  JR  ) , C22  ) ,  (  C ( 39  ) . C23  I .  ((.(  AO  ) .  AN  I  ) .  ( 0  (  41  ) .  TG  > 

•■'.III  VAI  f  NCF 

1(.  (1119), XM).  (C(llO).  XY2I.  (C(ll  1  >.PM  I  .  (C(1  1?  I  .  1  AII1  >  , 

21.(113). Tai'2).  (C(HA).(H).  (ClitS). TO 

DATA  I<rcn/SY.29S77V5/ 

N4MIU  *■  I/iEAmS/NXS.  I  PXr..E*,J.I:M*  .  TOTHS,  T  |Mi-  ,  HFI  T 
ll  I  R)  7  |l.  f ...  'I  .  It  )  til  Til  N'l 

r  CtTJ)  IHHi|.«iM  C  ( 42 )  Al'F  RESFHVFi)  rui«  IMIS  SUHHOnTINE 
i:  <:(A2)  IS  IJI'ITT  COHNOl.  SET  Til  »I.O  TiJ  |N'I.PIH  i.rifi 

Tl  AP*-E  ZS/°YS 

ft  AY*i<vs/K»r; 

4  /  Sll  *  II  •  (1 
4  I  *|i  *11.0 

SIE'Sll.-SlMI  I  ( JT/NKifi) 

C'lSS  I  O*(;0S(  T  (31/RM  III 

ll  (C  <  1  (13  )  .1  0.1' .  0  )  (All  I'.wAiiiM  NXS,  Fk  Y,  FRZ  • 

C 

C  l.<ACK|E((  iX'ilN  ANliirS 

•ill  CmN  1  I  Nl  F 

II  (C(TOI).UI.lt.ll)  no  TO  H 
C/ll  OMANII  (9XS.FiM.FH7) 

i,  A I  l  i;  III.  IkMU.-l .  II.II.F.  PY.EH2.F;IXN,E«Y'',Ffl2N,Y#U.Rl.l  I.  ,P|  IC'D 
Oil  SIKTP  (0.1, Fifth. I.HYM.FH/  l.FXS.iYS.E  Zb.UN.F  I  A  ) 

•  MM-  <>S*fcXS 

i  ISS«NVS»FVS 
/I.  I  «.S  =  F  7S*I  ZS 

1C  s? 

-,1  F  .1  IISHES 
Y  .1  SS  =  F  YS 
7  ISSHZS 

'V  ll  AP  =  -/MISR/I  S''M 
I.  ATHE  ISj/I  S'in 

ii  (ahfiti  him.m  ,iARM;r)  nn  in  a? 


SFFK"0>0 
SEEK  i)  060 
SFFKOOOO 

SEEK0110 
SFFKOT  20 
SEEK "1 30 
SFEKOiau 


ST  F K o  1 5 n 
SFFKO)  All 
SFFKhT  7U 


Rf FK020D 
SI  FA022I1 
SF  F K 0 2 1  ll 
SI  FA  ?Yo 
Si  FKi  24(1 


SEFK1250 
RrrK’'2A0 
SEFK"?70 
S  F  F  "  0  2  A 
S(  FA  "  29i) 
SI  FK.  0 
SI  FK  "3?0 
SFFMIJ3II 
SI  F  A  n  3  4  ll 
SEFMUAll 
SE  FK  »  37  0 
SI  F K  0 3 R 0 
SI  F A  o :> g ii 
S  l  F  A  0  A  0  II 
SF  F  A  o  4 1  U 
SfFM'OPO 
S|  FA  431) 
SI FA0440 
SEFA''4R0 
SI  FA  4 AO 
SEFAi  m *7 1| 
SI  FK'  4Aii 
SEFA'  4 9 il 

si  rAiisnn 
sn  F.t5i 
si  fm*  ?o 

SI  FA'«MH 
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TABLE  VIL  UNIVERSAL  SEEKER  FORTRAN  LISTING  (CONCLUDED) 


It  < ARSI TEA Y ) .OF . tARRUF >  liO  T'l  62 
CALL  DjriTFAP.IlTEAP.tFAPO) 

call  <iif<  tfay,  utfay.'ti  ayo  >  ' 

It  (DHAP)  fH1.A2.fll 
Ml)  CAIL  SFTFSIITEAH.iIHAH.-TARNUE  t 

mi  io  H2 

»1  CAll  SHIEST  (TEAP#  II TFAP,  TARlfGEF 
»V  II  <  OTFAY )  B3.«6.b£ 

H3  CAI  L  bHTEST  (  TE  AY . IITEAY.  -  T ARNUP  ) 

Oil  in  f5 

'14  CAl  L  SHTESTlTEAY.IlTEAY.  TARNGFI 
*5 <  CuNI  I Nl't 
Cil  10  66 

ttAP*#.0 
!IIFAP  =  P.O 
If  A  POxI*  •  0 
Ft  AYxO.O 
l|IFAYr|..0 
It  AYOMI.O 

Ci/K  I  ROL  CO'IHANC 
.1*1  Ci.'MINIE 
CM  L  1  HACK 
CM  1.  HYPO 

CAIL  SFKTRI1  .  1  .-VX.-VY.-V/,  VRXS.  VHYS.  VR7S,l'tl.t  TA  1 

CM  L  VFCIV<VRXS,VRYS.VR/S,«XS,,»YS,R/S.wxS.HYS.W/S.nMX:;.liRfS.DH7M 

CM  I  SI  TFST(KXS.I)H«S.C(.I)  ) 

CAIL  SPTEbT(RXS.0HXS.C(9)  ) 

CM  L  SFTESTIRXS.  I1RXS.  IP.  ) 

If  Hvm.RT.O. •>>. AND. (FSIIrt.GT.il>. "It  00  10  14H 

ti  i  it  =  -(ES'H/ni'xs  i 

Il-Mixt.n 

YrtlSS=YMlSS*nTn»ncYS 

/MSS«7Hlss*nin*nn7s 

VI6I,  >xSf)R?IV(51 >»«7»VI*7>«»7> 

V<77>xr.(9)/VRXS 
Vl7i:>*L  Ut/VRXS 
VI79>x>TAM2(VC1  )  ,  V  I  29  )  ) 

V*m  txA  IAH?<  V<  Jil).V<7‘>)l 
111'  Call  IHTERIIlRXS.n^Yl.uRXSf.HXSit.RXS) 

CM  L  IME><(l'RYS,nnYl.llPY2.MVSi|,RYS1 
CAL  L  IMEiM|il»ZS,IIK/l.llF72.'»7S').R/M 
H|*VI2V>»C!'SSie*Vll  >  •  S  I H  S  I G 

R. ix  VI. II'  ) 

<»*.*- V<79>»S  I  NSIIi*  VII  >»COSSlfi 

vuvi  ja)«i;ossio.  v(.mi«sinslc 

V.lx  v  ( .17  ) 

V*x-V(.Ah)tSIIISIG*V(  JO  »«cnssiG 
f.".l*RJ-RI»V  J/VI 
F*>|r*HK-Kt»VK/VI 
TnTHISsS0RT(FMJ*,7*Ft1A»«?  I 
IMRXS.LE.l*..  >  WK  I  IF (  6  .  HtAMS  > 

1  •>!'  Ri  I  UR  if 
C'lP 


RIFK''*.5n 
St  E  KM  *  A  (J 
Sf  EKi>67FI 
St FK  *>71 
StFK  672 
SFFX  *>73 
SI  FK  *>74 
SHFK  *>75 
STFK  *>76 
SrEK  *>77 
SI  EK  *  711 
SF  FK  679 
SFFKflMA 
StEK"690 
SFFKAAMl 
St  E  K  i*  A 1  II 
Sf  FK"A2ll 
S  F  F  K  '■  A  3  H 
SI  FKI'AAII 
SFFK  i/  A|) 
SF  FK  aa  p 


SF  FK  I  70 II 
SI FK1 ??0 
SI  F K 1  ?  2 2 
SF  FK ' 724 
SI  F  K  ’  7  2  6 
SFFKI2AK 
SI  FK  1  24H 
SI  F K 1  7 A II 
SI  F  K 1  7  9  fl 
SFFKi.MHl 
SI  FK1  .M  h 
srr-K  1.179 
SlF K'13n 
srrK’.Mu 
SI  FK  1.169 
SI  F K  t .! A fi 
S  f  F  M  1 7  ll 
SI  FM.IAU 
St  FM190 
SFFK  1  40 1) 
St  FK1  41  li 
srrM4?» 
Sf  FK1  43,1 
SI  FK  1  4 4 U 
SF  FK'  4SII 
SF.FK  i  4 A II 
SF  EKl  47  II 

rrFKt4Ao 
SI FK I  400 


\ 
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TABLE  VIIL  TRACK  ROUTINE  FORTRAN  LISTING 


A2  HA V 1 1  TRACK! R  A -A  KITH  LAO  OHLV 

subroutine  track 

COMMON  /SSAH1/  RFAD.OELI.AUTOT. TIME 

COMmO*  /S$&m2/  V _ 1231.1  *  T  (251UB _ L21L1-  -  .  - 

COMMON  /SSAN/  IEN0. NO . TNEXT . VH I N , STMMX . $12345, SUH222 

1. CFTA.SMA.CNU.SNU.TMA)(.NZ.LNV<5I)iTIUE<25D>#DLITO,RI  I ILE<9> 

2.  IFOEN. iHrOEN.MFOEN?, IFG2N 
COMMON  /GVR/VSY,VSM,MK$,NVS,N7S 

COMMON  /TRAKZ  /  It APO . TE AYR, OTF Af . D1EA Y. NTSC , xZSC 

equivalence  _ 

1 <  V< 17 ) . aZCMD ) ,  ( V ( 1 A )  <  A  VCMD ) •  ( V ( ?5 ) . TE AM ) .  • V<  2*  >  #  TEA Y I , 

2<V<?2),RXS)>  (V(l 1 2  )  >  IEAYO ) ,  ( V( 113 ) . TE AMD ) ,  ( V( 1 16 > . VSYM) . 


(Cl 144) >  TS) . 
(CUlAIjTLBr). 


( C 1 146  )  •  OHFOLU )  ■ 
(CIHV)i  ROM ) •  . 


(C(lNl).TLQV).  ( C 1 1 52  ) ,  SMOT ) 


1 <  V ( 1 7 ) , AZCMD ) ,  ( V ( 1 A )  <  A  VCMD ) •  ( V ( 25 ) . TE AM ) .  I V( 26 ) »  TEA Y  ) , 

2<V(?2)*RXS)»  (V(1I2).  IEAYO).  <  V(  1 13  > .  TE  AMD  > ,  ( V<116  )  .  7SYMI . 

3( VI117). VSMM) 

EOU 1 VALENCE 

1 (C< 142 ) , SK  > .  ( C ( 1 43  ) . AK  T ) ,  (CI144I.TS).  C Cll 45 > . OHFQLU ) . 

2(C(  146 ) .  QKK  )  ,  <C(  147  I.RIiSl.  (CUIMjTLBDj  (CU4V  ) »  TLOM  ) .  . 
3(C<150), TLDY).  (C( 151 ). TLQY )<  ( C ( 1 52  ) . SMOT ) 

NAMELIST 

A/NAMB/ 

B TEA Y. TEA YD. TEAM. TEAMD .  TEA YS . TEAMS. SUM1 . VS YM, YSMM. AYCMT*. AZlhD. VSY. 
CVSM.TIME.U1  VaVO.OLVSMO.ILOM.TLRM.TLDY.TLQY. TEAMS . IEA Yfc . NKAHQ. SMXS 

DATA  RTOD/57. 2957795/  _ 

SMXS«RXS/SMOT 

CALL  rOENKRDUMG.SPXS.RKAMO.-l ) 
ir(KFAO.NE.O.H)  GO  TO  10 
GO  10  14 
SUM1*BIAS 
0  TO  16 

imiME.LT.SUMl  I  Of  TO  20 
TEAMM1E  AM»RK  ANG 
IEATK.TI  AfRKAMl, 

TEAYD*RTOO»TFAYK 

TEAPI)*RT00»TFAMK 

CALL  FRENI  (lOIIHY.  IE  A  YD.  TE  A  YS,  •  t  ) 

'  CAIL  FCENKTOUMM.  TEAMD.  TFAPS.-l  ) 

C0NH1. 

SUMl*SUMl«TS 
CALL  TTEST(SUHl) 

VSYM»AKT«  TFAYS 
VSMP«AKT»TFAMS 
CALL  DIFIVSYM.OVSYM, VOUMY ) 

call  nir< vsMM.nvspp. vuumm) 

CAI I  LAU( VSPM.OVSPM.OLVSMO.CLVSP.OICVSM, TLOP.tiDUMP) 

CALL  LAG(  VSYM.  IIVSYP.  OLVSYO.PL  VSY.  Dili  VSY.  TLQY.  GDUNV) 

AYCMD«QKK»OLVSY 

AZCMD*-0KK*GLVSM 

CAI  L  LiHlTIGLVSY.IILnvSY.OMEGLD.-OMECLD) 

CAI  I.  LIMI  I  <  01  VSP.  RlUVSP, OMFGLO •  -OMFGLD ) 

VSV«RLVSV«SK 

VSP«GLVSP*SK 

IFICONT.EO.l.  )  MR  I  I  E ( 6 . NAMR  ) 

CONI =U. 

RETURN 

END 


LLOV  20 
1  LOO  30 
LLtO  40 
LLOO  50 
LLOO  60 
LLOO  70 
LLOO  SO 


LLOO  90 
LLOO  100 
LLOO  110 

LLOO  120 
LLOO  130 
LLOO  140 
LLOO  160 
LLOO  160 
LLOO  170 
LLOV  160 
LLOO  19o 
LLOO  200 
I  LOO  210 
LLOO  2?U 

I  LOO  240 
LLOO  250 
LLOO  260 
LLOO  270 
LLOO  260 
LLOO  290 
LLOO  JOO 
LLOV  310 
LLOO  32U 
LLOO  330 
LLOO  340 
LLOO  35U 
LLOO  360 
LLOO  370 
LLOO  360 
LLOO  390 
LLOV  400 


LLOO  490 
LLOO  500 
LLOO  610 
LLOO  6/0 
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TABLE  IX.  GYRO  SUBROUTINE  FORTRAN  LISTING 


ccyro  gyro  sine  rail  _ onut  option 

SUBROUTINE  GYRO 

DIMENSION  F019(4).M2E<4).F04C(3) 

COMMON  7SSAH1/  READ . OEL T , AU TOT . T I MF 
COMMON  7SSAM27  V  (250),  T  (250), C  (251) 

COMMON  /TRAKER/  COUNT. TN.Nl 

A _  -GFOLiOFry 

COMMON 

1/OYR/VSV. VSP.CE.MYS.  WZS 

equivalence 

1  (  V( 27  ) ,  UN  l ,  (  V  (  2 0  > ,  E.T A  )  , 

EQUIVALENCE 

1  ( V ( 53 ) •  WA j_,  UI’DtJLUj  .  ULLiiijJZJ ,  <Yl>Q)»P NAJ' 

2(V(57),PUY>,  ( V ( 50  ) , 0U2  ) ,  (V(59),0A),  (V(60).AV). 

3 ( V(61 ) , AZ  ) 

equivalence 

l(V(B5).nE>.  (V(06).nEXS),  m07).E>.  (  V  (  MB  )  <EXS )  > 

2 ( V( 69  > , Cl ) •  ( V ( 9  0  ) ,  Q1 > ,  ( V (91 ) , D01 ) ,  (V(V2).U2). 

„i(V(93),n02),  (V(94  ),fll|L»j_  (lLS5J,Q2N)t  <V(V6),frt), 

4  (  V( 97  ) ,  |)EE  ) ,  (  V  (  98  ) , FEXS ) ,  ( V ( 99 ) . DFEXS ) 

EQUIVALENCE 

1 ( C( 1 16 ) , N3S ) .  ( C ( 1 17  ) , K2T ) ,  (C ( 118 >. DUMP > ,  (C(tl9>.HAIL> 

2 . ( L ( 136 ) . GNUT ) ,  ( C ( 1 37 ) . DTR I ,  ( C ( 138 ) .  DST ) ,  (  C ( 139 ) . DSU ) . 

3 ( C ( 140  )  > DAN ) .  ( C ( 141 ) , DDU ) ,  (C( 153 ) , CF1 ) .  (C( 154 ) , Cf ?  ) . 

4  ( C  ( 155 ) ,  Cp3 )  . 

REAL  IIP. I1RXE. I2E, l?S,  I2E>S. I3S. I3T. |4C. I4|i. I4DXC.  I  It 
REAL  MC.MP,LR.KO18.MI2.K12F,K04C,NPX.K3E.K3S.K2fc.K2S,K?EXS.KlE. 
1K1PxE.MJ.M23.M123.X2T 
REAL  KRR.KRT.KGT.ftr.K.M7JS.M123S 
NAME.LIS1 
a/nam  it. 

HUE) ,E.FHS.FFXS.ET.EXS.EXSS.EXST.QA2.0BP.GEP.MUXCD.KUR.Kl,T. 
CKI2T.KNH.KRI.M123S.M23S.SE.N|3.XE8M.X6EX8.X0Sfc,kR0.XKK. 

PXKRM.  XTfc.  XII 
E/NAHA/ 

lAl.A2.AGB,A0E.R.Bl.P2.Cl.CA.CV.CD.CE«Cr,rH.C0SR. CUSC  .COSD.COSfc. 
2C0UNT.CX1.CX2.CX3.  PIN, P2N. RB . 0*X, DC. ODB. DDE1 . DD»C « DL . DE 1 » 

30EfcXS.PFX.nFXS.ilFE.UFL XS . DOl , DP2 . DX . DMA. DNBXE. DHDXC.  E . El • EXS . 
4(R.)C.rUfl.FDC.r0F.rE.rFfc,FEXS,Cl,01N,ClZ.02.Q2N.S2Z.QH.UR2,OBL. 
5GBXE.UC.GE.0EL.CFXS.QS,HQDXC.HC.MR.REl.RF2.RXl.NX2.SRSF,Sfcn2.SFLR, 
6SECC. SfcCP. SFCE. S I NB.S INC, SINS. SINE. STAN,  TAN2. TANC. TANE, 

7IANY.TE.TEA.TEDU.TE)  . TEP.  I E S . TEU . TF XS . TEXS A , Tfc XSF , T!  XSP. fEXSS. 
BTEXSU,TXDU.Ul,U12.N2.MA,N8XE,ND,U0XCiNErUEBXE.HN.HT.  xFM. 
9XMK.XMP.XI-M.XIMC,  ES.XS 
A. TAND. TIME.DELT.CX4C 

MAVERICK  OYHO  FSTIMATLI)  AND  CALCULATED  PARAHfc (EPS  (7-15-69) 

DATA 

I I3S. IJT, I2E. IPS. I2EXS/4.27E  4.7. OAF  4.2.36E  4.P.98E  4.2.70E  4/. 
2I1H.IIF.I1BXE/1.46F  4.1.9JE  4.P.94E  4/. 

3I4C. I4P. I4DXC/1.75E  3.1.94F  3.2.8BF  37. 

4M3.MPJ.M12J/1.12E  3.2.53E  3,P.*»5E  37. 
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TABLE  IX.  GYRO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


5K3F  .XJS.«lF;KlRxr/A.42£-7.2.10r-7,0.52E-7,1.0iir-7/» 
6K2E.K3S.K2FXS/B.3AE-7.A.7BE-7.H.7BE-7/. 
7HCF.R2L.0CA.ePL/lA. . 4540 . , 2 A . , 4*4B ./ . 

8FS1H/33. 0,7.3,1.9,0.22/, 

9F12fc/22.5,2.1,fl.5,0.16/. 

AFMC/ll.a,1.3.A.4/. 

RY.K(i2.KOlB.K12E,K04C/1.2.l5d.,l.,9.,0./. 

CVAlB.yilE, VAAC/A.B.l.bx?,!/^,  .  .  _ 

iillS.UBXE.UE.UEXS.UD.UDXC.HPX/l.Al.l.AA.l.AA.l.OA.A.Sb  . 
ElR. 0L.0T. DUB. THE T/0. 215. 12. 1.21 .0372,1.90.0./, 
FO.A/9BO. .57.2957795/ 

IF  (READ.EO.A. )  00  10  3 

C 

CCCCC  INITIAL  COMPILATIONS 

if 

M23S«M23»«2 
Hl?3S*Ml23««2 
EXSS«|2F*S-I2S 
HPXCD*MDXC*I40 
ET* I2g* 14T  _  . 

BEE* I IBXE- I IE 
EEX$*I2E“ I2EXS 
SE*I2S-I2E 
E  XS1 *  1 3  T ♦ I ?t  XS 
ERStllR. I2S 
H  1 3*M3S* 1 35 
0HP*GPl/M123 
Cfcl *0PL/M23 

KRA*4.F-7#RCA/RPL»5.E-IU»RCA 

NR1  *1.4E-3/(RPL«RCA)0.9E-7/ SORT!  RCA  > 

K01*1.4E-3/<OPl«OCA)4l.VE-7/SOAT(OCA> 

X0R*9.E-7«GCA/0PL«9.E-l#»GCA 

XXP*KJE-X3S4XRR-AHI 

XKG*K2F*XGT-KGH/2. 

XKPN»XAR»N3««2 

XGSE*XPHH.H23S«(XAG-X2S> 

XGEkS:P?3S*(XK0-K?EXS) 

XTf *XKNH»M?3S»<N2fc*S-K2S) 

XEPH«M12JS*(K1F-K1HXF  ) 

0A?*GA»«2 

XII*0UB»N3S«»2/9BB. 

(U*GNU1 • XU 
WRITE  <6, NAM?) 

IF  (  PA  |l  ,  EO .  1 . 0  )  GO  10  29 
H*liN 
E1«ETA 
GO  10  26 
?5  CONI | HUE 

b*fia 

El  *0N 

26  CON  1 | NOE 
EXSxR 
fc*El 
ES*E 
<S*EXS 
KOI'NT«B 

KDUMP*ni)Hp..01 

c 

ccccc  missile  frahf  rates  and  accelerations 
c 


.50.2.0/. 


1 
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TABLE  IX.  GYRO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


3  CONTINUE 
K0UHT«H0UWT»1 
0S2>GB**2 
agb«abs(0b> 

ir  (A0B.lE.0RP)  00  ?0  21# 

OgisAOP-OBP 
00  10  ?20 
210  OBI *0 . 

220  CONTI  NOE 

)  B*f  0 1R  ( 1 )  *  I  01B(2)»ORL«FI1H(3)»SORT<OA?»OC««2>*I  •1BM)»UBL»»2 
FB*F8»CF1 

FC*I  04C<1  >«CF2MFU4C(2  )»AriS(0C  >  ♦FMC ( 3  >«$QR I  <  UA2»0B*  >  >*CFJ 
C 

C  HlSbJLE  BAIL  POSITION 

IFIKAIL.EU.1.0 )  00  10  40 
C  BOTTOM  RAIL 

U8-BY 
WC*WZ 
DUHsDNY 
OMCeOHZ 
OB*  AT 
0Cr  AZ 
b*Uh 

ei*tTA 
VfcsVSY 
YEXh*V*P 
00  10  42 
40  CONTINUE 
C  SIDE  RAIL 

-<B  =  wZ 
MC  =  HY 
OMH*0H? 

0MC«0riY 
«1H  =  AZ 
UC*AY 
N*F  TA 
bl*UN 
VE*VSP 
VE*b*VSY 
42  CUNUNIlE 

MEEXS*IiE*IIEXS 
S I  nb*  S I N( 0 1 
C0SH«C('S  ( U ) 

S  I  NT  *  S  I  N I  El  ) 

COSt«COS(El > 

SECF*l./COSE 
TANt«  IANIF1  ) 

IAN2«  TANE#»2 
SEt2*SFCE»«? 

STAN*SFCE»1ANE 
b  I  NC*COSB»S  I NE 
CI*aRS*(SINC> 

C0H*C0S(C1 ) 

SECL«1 ./COSC 
TANC*TAN(C1  ) 
coMi*ccs8»cnsE/cnsc 
SEr.n*i./cosn 

S  I  No*  S I NH/COSC 

iANn*siNn/coso 

NEcWA»S|NR*MC«COSrt 
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TABLE  DC.  GYRO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


MBXE«UA*rnsn*MC*SINU 

UD*-MA*!;iNC»MB»COSC_ 

MDXC*NA*C(lSC*HB*SINC 
»WBXE»n«lA»COSB-DMC»SING-llB«NE 
'  bunxC*liUA*COSC*DHB«SINC*UC«HD 

UEHXE«WE»MBXE 

DDXC«un«HIIXC 

M01lXC»HBX5B#||XSL 

•  C 

CCCCC  TE  COMPUTATIONS 
C 

HC=-SIM)»TANE«(1.*IR»SINC) 

MP*C0SD«  <COSB-LR«SINB»»2»SINE)/COSB 
X|H*MC«"P  _  .  . 

•:e>secf«mbxe»tane*dexs 

CF*2.*S|ND*SECC*SINR«TAN2 

CX?«S |NC*SIND**2«BE »cf*bexs . 

CA«FBS«SEC2»TANF 

CXl«TAPn*0HDXC«SECD**2*UDXC-TANC*SECC*0EXS**2 
CD«SECE*<  I1B*  (  TANE»DWUXt=  SEC2«BEIXE  MBEEtMEUXt ) 
CH«I2S*TANE«(SFCE*DMBXE-STAN*HFBXE> 
CX.1«TANC»SECC»DEXS*CF#DE 

CR*IANn«DM0XC-SErn««?»nDXC-SINC»ISIND«DE>«*2 

HE1*X|P»CH*H0XCD«HP*0bXC 

KE?*EXSS*CF»nEXS-UI3»0EXS*X|M«r.)il«HDXtD«MP*0Dxr 

c 

CCCCC  EE  MS  CPHPUT  A  T  I  OKS 
C 

X|HC«|4C»MC 

RX1*CI)»i:m»SE  •SECE*WHXE»DF»MI.'*DE4X|PC«C8*M0*CI)»MC»0liXC 
RX?rCH»CH  •X|HC*CXl*hDXCD*OflxC«MC 

C 

CCCCC  PRECESSION  TORQUE 
C 

IEXSP>( VF*S(  CE«VEXS*PC)*K2( 

IEP*-VPXS»K2T*MP 

C 

CCCCC  UNHALANCE  TORQUE 
C 

oE*iiA»siNH»nn«nosH 

UB XE ■ GA • COSH »0Ca SINN 

HS»UA*COSB«cnSE»OH*SINE-BC«S|NH«COSC 

SRSMSINB«SINE 

SECH>SINE*fOSB 

OEXS»-OA#SECB*OR«COSF*OC«SBSE 

UU*UUXr»TUA»SINC-OR»COSC>*UD«(OA»COSC«OH»5INC) 
XHP*MX*(l. -COSH«LR«S INHaSRSE ) 

IEXSU*(liS»UIIXE«SFCE-UbXS*IANE>*GF-IIE*SFCF*GHXt‘»UU*Ml 

J»XPt,*HF«MPX*(la*I.M*SFCBI«SINB*OS 

fEU«-l|S»OEXS»UFXS»OS*UU»HP-XNP#GEXS 

C 

CCCCC  FRICTION  TORQUE 
C 

OH*-NH*WA*COSR«TANE-WC*S|NB*TANL»l)EXS*SErE 
IIEU-MC»COSB-NA#SINB*QF 
BC*COSP«PEl-SINPiSINF»PN 
f DR*S|fiN(1..0B» 
ir<PR.FQ.I.)  FOH*0 . 

F 0rrS|GN(1..0C) 

ifcdc.eq.o.)  f nn*u • 
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GYRO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


T\BLE  IX. 


F0E«SIGN(1..0E1> 
ir(DE;.to.n. >  fqj-.u. 

AOI <ARS(Gfc ) 

IF  (ARE.Lfc.GEP)  GO  TO  230 
OEI «A0F*GEP 
00  10  240 
230  OEL*«. 

,  240  CONUHl'E 

FE1<F12F<11<CF2 

F E2< (FI 2E( 2 ><0EL<F12E( 3 ><SOHT(GB2*OBXF**?MF1?F(4> <GEL<<21<CFJ 

FE«F£1<FE2 

XFM<FC<F0C<V94C<DC 

I EXSF  *- ( F8  <FnH*V01B<DB)<SECE-XFN<MC 
IE>*-(Ffc  -•£DE*012E»BEl)-J(FJ)tHP 
C 

CCCCC  SPHNQ  I0HQUF 
C 

TANY* I A  N ( Y<F1 ) 

XMX<KU4C<C0SR<TANY 

.  I 1*01  B<SELE<KU2I<  I  ANCUlT^lllKtfU: 

TES*- <K12F<K02<C0SH I  < T AN Y- XMK<MP 

C 

CCCCC  ANISOEl  ASTICITY  TORQUES 
C 

1EXSA<GS<CF<XGSI  -GF* RE  XS<IANE<XGEXS<GRXE •6E*SfcCt*XEHH 
1EAs-uEXS<GS<XTL 
C 

CCCCC  DYNAMIC  UNRAIANTR  TORQUE 
C 

NT=N3S<1 IMF*THF|/H 

r  xuii> xn« s  i  n <  u 1 1 

iEiiu«xMCo$<wn 

c 

CCCCC  riNAL  COMPIIT  A  1 1  ('NS 

c 

Tfc*G<(  U  P<OFR<TFF<RST< It S< RSU<TEU<DAM<TE A<0DIKTE0U> 

TFXS»«<( 1EXSP<UF R<T EXSF<OST<TEXSS< RSU<TEXSU<OAN< TtXSA«nnM<TXIlU) 

if(gnut.eo.i.o)  go  in  2»o 

C 

cccc  nvkn  without  nuiaiiun 

CX4c<l4C<nw II XC*SIN0/CQSn>|4C<Dr. XC<SECI<SFCO<H IiXCD<0I  xc 
rXG|a(nil<CH)<CX4C<M(. 

IEG I »CX4C<MP 
IE»S*IFXS<Txni 
TE*  IE*TFG| 

0EMEXS/NI3 

•>E»S«-TF/MI3 

CAIL  |NTFR(DF.0fcMl.DEM2.FU0,E> 

CAIL  |NTER(0EXS,nXNl.|IXH2.FX0.EXS) 

•  OAI  l  IMFHIDR.rpMl.URN2, ROO, R) 

CALL  |NTER(nFJ.neiMl.RFlH2.E10.Et ) 

GO  10  26« 

2*U  CON  I  I Ndf 
C 

CCCC  GYRO  WITH  NUTATION 
A 1  =  E  T  ♦  X  |  H<0I1SD 
H1*-X|M.S|N|I<TANF 
A2<XlMr«C0Sn 

R2M-XST-X  IMC<SIN0<TANE<  I  |«<SEC?<  1 2S<TAN2 
IIX*A1<R2-A?<RI 
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TABLE  IX.  GYRO  SUBROUTINE  FORTRAN  LISTING  (CONCLUDED) 


N1<H-A?*<EXSS*CE>U|3>XIH*CX3) 

i*J  >  »<CA.nE«t£XS.TAN£.DE»SE»TAIIEtPE.IIMC«CxaH/PX _ _ 

M2*CA1*<CA»DEXS»SE»CE«NI3-X|MC»CX2fEEXS*DIXSaIANE)*A?»|(|M*CX2)/DX 
01«lt2*<TE-RFL>-|l«UtJt$-RXUj/DX 
02»<A1«<TEXS-RX?)-A2#<TE-RE2>>/DX 
CALL  DIFI01. 001. DUMA) 

CAl L  0IFI02. 002. DUMB) 

CM  l  blL'UliJLl  JJPldlUiAKLOliLimmiUL 
CAl L  GRAlE(1.02.D02.G2Z.Q2N,D2N.0UHX2> 

CAl  1  OIFIDA.pjlB.OUNR) 

CALL  OIKUEl.ODEt.DUHU) 

CALL  QRATEIl. OB.  DUB  >  A2>  B »  DAX. DOM7  > 

CALL  OA ATE ( 1. DEI. U0fcl.E12.El.DEX. CURB) 

Jil?*  ll«Hl  .. 

HN^SL  'T1N12) 

rrL>(01»Ul«022)/H12«ES 
CALL  ninrrE.DFE.DliHC  ) 

FEXS*<P?-W?»012)/W12»XS 
CALL  DIHHXS.OFEXS.OUMU) 

.  CAU  I  DSFC(FFE.nFE.E2.E.Dfc.  It.  U^iJUUl.iBllJUrJ.  UMX3J 
CALL  LOSECIFEXS.  DFEXS.EXS2.FXS .DEXS>B..l..k' DUNX5. DUHXa ) 

26V  CONTINUE 

IF(PAIL.Eu.l.l)  00  TO  2711 
MYS*DEX$ 

U2S>0E 
UH*H 
ETA*F1 
00  TO  ?H n 
2 /II  CONTINUE 
UYSzUE 
N2s«nExs 

UN‘H 

ETA*« 

?tt(l  continue 

K  X*MOD 1  FOUNT. AOUNP) 

IF  ((KX.tO.Q  )  .OR.  (COUNT. LE. 5.  )  )  NR  I  TE  ( 6.  NAHA  ) 

RETURN 

END 
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(The  reverse  of  this  page  is  blank.) 


UNIVERSAL  SEEKER 
SUBROUTINE 


TRUCK 

SUBROUTINE 


RXS  RYS  RZS 


AYCMO 


12.  Universal  Seeker/ 
b routine  Block  Diagram 
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rse  of  this  page  is  blank.  ) 
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Figure  14.  Universal  Subroutine  TRACK  Flow  Chart 
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CONTINUE 


Figure  15.  Gyro  Subroutine 
Flow  Chart  (CAS  6  DOF 
Simulation ) 
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(Ihe  reverse  of  tnis  page  is  blank.  ) 


FRICTION  TORQUt 


Figure  15.  Gyro  Subroutine 
Flow  Chart  (CAS  6  DOF 
Simulation)  (Continued) 
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(The  reverse  of  this  page  is  blank.  ) 


AGE-ABS(GE) 


GEL=0.  GEL-AGE-GEP 


Figure  15.  Gyro  Subroutine 
Flow  Chart  (CAS  6  DOF 
Simulation)  (Continued) 
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(The  reverse  of  this  page  is  blank.  ) 


return 


TABLE  X.  SEEKER  (MSEEK)  SUBROUTINE 


Coordinate 

Name 

Quality 

Units 

System 

V(  1 7) 

AZCMD 

Azc>  Elevation  maneuver 

g's 

Autopilot 

command 

V(  18) 

AYCMD 

Avc ,  Azimuth  maneuver 

g's 

Autopilot 

command 

V(22) 

RXS 

Rx,  Seeker  boresight  range 

ft 

Seeker 

V(23) 

RYX 

Ry,  Seeker  lateral  range 

ft 

Seeker 

V(24) 

RZS 

Rz,  Seeker  normal  range 

ft 

Seeker 

V(25) 

TEAP 

«p,  Tracking  error  angle, 

rad 

Seeker 

pitch 

V(26) 

TEAY 

€y,  Tracking  error  angle, 

rad 

Seeker 

yaw 

V(27) 

UN 

v,  Seeker  elevation  gimbal 

rad 

angle 

V(28) 

ETA 

rj,  Seeker  azimuth  gimbal 

rad 

angle 

V(53) 

WX 

u>x 

V(  54) 

WY 

(  Missile  body  rates  in 
y  autopilot  axes 

rad/sec 

Autopilot 

V(55) 

WZ 

V(56) 

D  WX 

Ux  Scalar  components  of 

V(57) 

DWY 

.  missile  angular 

Y  acceleration  in  auto- 

,,  2 
rad/ sec 

Autopilot 

V(58) 

D  WZ 

uz  pilot  axes 

V(59) 

AX 

Ax  Propulsive  and  aero- 

V(60) 

AY 

dynamic  acceleration 

Y  components  autopilot 

g’s 

Autopilot 

V(61) 

AZ 

Az  axes 

V(62) 

VX 

V ' 

x  Missile  velocity 

V(63) 

VY 

Vy  components  in  auto- 

ft/  sec 

Autopilot 

V(64) 

VZ 

.  pilot  axes 
v  z 

V(66) 

TOTMISS 

Total  Miss  Distance 

ft 

Miss  Distance 

V(70) 

EM  J 

Y  Component  of  Miss 

ft 

Miss  Distance 

V(7 1 ) 

EMK 

Z  Components  of  Miss 

ft 

Miss  Distance 
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TABLE  X.  SEEKER  (MSEEK)  SUBROUTINE  (CONTINUED) 


Name 

Quantity 

Units 

C(13)  TARNGE 

~~Y~  »  Half  the  seeker  field  of  view 

rad 

C(  14)  GK 

K  ,  Guidance  gain 
,g 

t^,  Tracker  time  constant 

g's/deg/ sec 

C{  1  5)  TAUAP 

sec 

C ( 1 6)  OMEGAL 

Precession  rate  limit 

rad/sec 

C ( 1 7)  Cl 

Seeker  drift  term 

rad/sec 

C(  1 8)  C2 

Seeker  drift  term 

rad/sec 

C(  19)  C3 

C(20)  C4 

Seeker  drift  term 

Seeker  drift  term 

rad/sec/g 

C(21 )  C5 

Seeker  drift  term 

1  /  sec 

C(22)  C6 

Seeker  drift  term 

1  /  sec 

C(23)  Cl 

Seeker  drift  term 

sec 

C(24)  C8 

Seeker  drift  term 

sec 

C(25)  C9 

Seeker  frift  term 

rad/sec/g 

C(26)  CIO 

Seeker  drift  term 

rad/sec/g 

C(27 )  Cll 

Seeker  drift  term 

1 /sec /g 

C (28 )  C 12 

Seeker  drift  term 

o 

rad/sec/g 

C(29)  C 13 

Seeker  drift  term 

2 

rad/ sec/gc 

C(30)  C  14 

Seeker  drift  term 

2 

rad/ sec/g 

C(3 1 )  C 1 5 

Seeker  drift  term 

rad/sec 

C(32)  C16 

Seeker  drift  term 

rad/ sec 

C(33)  C 17 

C(34)  C 18 

Seeker  drift  term 

Seeker  drift  term 

rad/sec/g 

C(35)  C19 

Seeker  drift  term 

1  /sec 

C(36)  C20 

Seeker  drift  term 

sec 

C(37)  C21 

Seeker  drift  term 

rad/sec/g 

C(38)  C22 

Seeker  drift  term 

rad/sec/g 

C(39)  C23 

Seeker  drift  term 

rad/sec/g^ 

C(40)  AK1 

Kj,  Tracking  loop  velocity  gain 

1  /  sec 
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TABLE  X.  SEEKER  (MSEEK)  SUBROUTINE  (CONCLUDED) 


Name 

C(41 )  TG 
C(42) 


Quantity 

Units 

Gimbal  preload 

g 

Drift  control,  set  to  1 .0  to  include 

drift 

Not  used 

TABLE  XI.  TRACKER  GLOSSARY  OF  TERMS 


Name 

Quantity 

Units 

Coordinate 

System 

19 

V  Array 

A  ,  Elevation  maneuver  command 
zc 

g 

Autopilot 

BBS 

A  ,  Azimuth  maneuver  command 
yc 

g 

Autopilot 

V(22) 

Rx>  Seeker  boresight  range 

ft 

Seeker 

V(25) 

Tracking  error  angle,  pitch 

rad 

Sepker 

V  (26) 

€y,  Tracking  error  angle,  yaw 

rad 

Seeker 

V(1 12) 

TEAYD  Tracker  error  yaw 

deg 

V ( 1 13) 

RKAMG 

TEAPD  Tracker  error  pitch 

deg 

V(1 16) 

RKAMG 

VSYP  -  Tracker  output  signal 

deg/sec 

1  ' 

V(  1 17) 

Pitch 

VSPP  -  Tracker  output  signal 

deg/ sec 

C(  142) 

Yaw 

C  Array 

SK  Torque r  gain  coefficient 

V/deg / sec 

C(  143) 

AKT  -  Tracker  gain  constant 

lsec 

C(  144) 

TS  -  Sampling  period 

sec 

C(145) 

OMEGLD  -  Precession  rate  limit 

deg/sec 

C(146) 

GKK  -  Guidance  gain 

g / deg/ sec 

C(147) 

BIAS  -  Sampling  rate  offset  bias 

sec 

C(  148) 

TLDP  -  Tracker  filter  lead 

sec 

C(149) 

time  constant  pitch 

TLGP  -  Tracker  filter  lag  time 

sec 

C(150) 

constant  pitch 

TLOY  -  Tracker  f  ilter  lead  time 

sec 

C  ( 1 5 1 ) 

constant- Yaw 

TLGY  -  Tracker  filter  lag  time 

sec 

C(152) 

constant- Yaw 

SPOT  -  Tracker  spot  size 

ft 
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TABLE  XII.  GYRO  GLOSSARY  OF  TERMS 


V(85) 

V(86) 

V(87) 

V(88) 

V(90) 

V(91) 

V(92) 


V{96) 

V(97) 


Quantity 


V  Array 

v.  Seeker  elevation  gimbal  angle 
rl,  Seeker  azimuth  gimbal  angle 
u> ' 

X 

,  Missile  body  rates  in 

wy  autopilot  axes 


Scalar  components  of 
missile  angular  accelera¬ 
tion  in  autopilot  axes 


Units 


Coordinate 

System 


rad/sec  Autopilot 


rad/ sec  Autopilot 


x  Propulsive  and  aerodynamic 

A'  acceleration  components  in  g 

autopilot  axes 
z 

DE  Total  yaw  precession  rate 

DEXS  Total  pitch  precession  rate 

E  Yaw  gyro  inertial  angle 

Cl  Yaw  look  angle  (indicated) 

Gl  Forcing  function  cross- 

coupled  equation  1 

DG1  Derivative  forcing  function 

cross -coupled  equation  1 

G2  Forcing  function  cross - 

coupled  equation  2 

DG2  Derivative  forcing  function 

cross-coupled  equation  2 

GIN  Integral  forcing  function 

cross-coupled  equation  1 

G2N  Integral  forcing  function 

cross-coupled  equation  2 

FFE  Forcing  function  yaw  axis 

DFE  Derivative  forcing  function 

yaw  axis 


Autopilot 
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TABLE  XII.  GYRO  GLOSSARY  OF  TERMS  (CONCLUDED) 


Name 

Quant  ity 

Units 

Coordinate 

System 

V  Array  (Continued) 

FEXS  Forcing  function  pitch  axis 

V(99) 

DFEXS  Derivative  forcing  function 
pitch  axis 

C  Array 

C ( I 16) 

3S  -  Gryo  motor  Speed 

rad/sec 

C ( 1 17) 

K2T  -  Precession  torque 
coeffic  ient 

gcm/V 

C ( 1 18) 

Dump  program  control  logic 

B  =  0 

C(1 19) 

-  Rail  control  logic 

S  =  1.  0 

C(136) 

GNUT  -  Program  logic  control  - 
W/O-O,  W  =  1.  0 

C(137) 

DFR  Coulomb  friction  drift 
f  actor 

Dim 

C(  138) 

DST  Spring  torque  drift  factor 

Dim 

C(  1 39) 

DSU  Unbalance  drift  factor 

Dim 

C(140) 

DAN  Anisoelastic  drift  factor 

Dim 

C  ( 14 1 ) 

DDU  Dynamic  unbalance  factor 

Dim 

C(  1 53) 

CF1  Friction  factor  coefficient 

D 

C(154) 

CF2  Friction  f  actor  coefficient 

D 

C  ( 1 55) 

CF2  Friction  factor  coefficient 

D 
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TABLE  XIII.  SUBROUTINE  GWAND  USED  TO  SIMULATE 

AIMPOINT  WANDER 


(This  option  is  exercised  when  C(106)  ^  0) 


For 


C(  103)  <  RXS  S  8.35*  C(103), 


ERY 


_ H*RXS _ 

C( 1 03)*  C(  1 05)*  C( 106) 


ERZ 


-  Z  *'•'  RXS 

C(  103)*  C ( 1 0 5 ) *  C(  1 06) 


Apparent  target  motion 
is  y  and  z  earth  axes 


Othe  rwise 

ERY  =  ERZ  =  0. 

Where 


H  =  fj(a) 
Z  =  f2(a) 


functions  f]  and  are  described  by 
function  generators  1A  and  IB 


and 


a  =  C(  104)* 


V 


-0.563  +  %  2.45  -  2.42 


/,  C(1 03  )\ 
\  "  RXS  / 


or 


a  =  C(104),  whichever  is  smaller 
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TABLE  XIV.  SUBROUTINE  GWAND  USED  TO  SIMULATE 

TARGET  MOTION 

(This  option  is  exercised  when  C(106)  =  0) 

y-p  is  target  displacement  in  the  positive  earth  fixed  y 
direction. 

yT  =  V £  [t  -  r  ( 1  -  e_t/T)] 
yT  =  Vf  (1  -  e_t/T) 


Where 

a  =  C(103)*  32.2,  initial  target  acceleration 
V£  =  C(104),  final  target  velocity 
t  =  Vf/a 

t  =  time  measured  from  the  point  when  boresight  range  equals 
C(  105) 


nnonoonnon 


TABLE  XV.  SUBROUTINE  GWAND  FORTRAN  LISTING 


f  FORTRAN  DECK 

rw*HD  AIM  POINT  UANOEH 

SUPROOTINE  GWANlMRXS.ERY.ERZ) 

COMMON  /SSAhl/  KFAO.nfcLT.AUIOT.TjHE 
COMMON  /SSAI'2/  V  <25 0>,T  <2501. C  (250  ) 

FOUIVAlfcNCE  (C(lnA).Rr). <C<1<H).A).<C(105).PLOTK),  ( C( 106 

II  C ( 1 0  3  > IS  SET  TO  0,  THIS  suhrouiine  will  oe  bypasseu 
if  c ( i o a > is  non  o.aim  point  hander  kill  be  simulated 
if  c < t o 6 ) is  set  tu  o,  tablet  notion  will  re  simulated  where 

1)  C(10J)  =  |N|  T A L  TARGET  ACCIL.  IN  OS 

2)  C(104)*FINAl.  TARGE  I  VELOCITY  IN  FPS 

3)  C< 105)=StE*ER  RANGE  AT  START  OF  1ARGET  MOTION 
TARGET  MOTION  OUEYS  THE  FOLLOwlNQ  EGUAT I  OKS 

DDYs  A*EXP(-T/TAIJ» 

DYs  T»A«(l-l XPI-T/TAU) ) 

Y  s  T  •  A  <T-TAU*<l-EXP<-T/rAU>n 
IF  (RC-aO. £0.0.0  ICO  TO  50 
GX  =  l'I.OTK*PHOTUK»RF 

ERYsO.O 
ER2*0 • 0 
SNA-0 . 0 

CALL  FGENt  < IA.SNA.H.-1> 

CAll  FGfc’Nl  (  I  FT.  SNA  •  Z  .  -1 ) 

IF  <C(l(i3).EO.O.O)  GO  TO  100 
RTST*fl.35»RF 


WAND0010 

KAND0020 


)  . PHOTOK ) WANOt  050 
WAND0060* 
WAND0070 
WAND00P0 
WANDQ090 
WAND0100 
KAN00110 
KAN00120 
RANU0130 
WANU0140 
WAND0150 
WAN00160 
WANU0170 
KANU0180 
RANU0J90 
WANU0200 
K  A  N  li  0  2 1 
UAND022 
WANU0230 
KANU0240 


AC=32.2*C(103) 

T  A C  =  C  < lo4  ) 

T  Alls  T  AC/ AC 
GO  TO  100 

50  IF  <C(103).EQ.0.0)  CO  TO  100 
IF  (C(106).EO.O.O )  GO  TO  200 

ir  (pxs.ot.htst)  00  to  100 

IF  (RXS.LT. RF)  GO  TO  100 
RMO*<  \  .li-RF/RXS) 

SMA=A»<-. 563 ♦ SQM T<2. 45-2. 42«RH0)) 

,  IF  ( SNA . GT . A  )  ShA  =  A 
CALL  FGENK  IA.SHA.H.-1) 

CAI  L  FGT:N1<  IB. SNA. 2. -1) 

ERY=H«RXS/GX 
FHZ=-Z*RXS/OX  , 

100  RETURN 

200  IF  (RXS.LT. C(105)l  00  TO  250 
TSTART«TIHE 

CALI  II  IF  (RXS.pRXS.RXSO) 

CALL  SPTESK-HXS.-DRXS. -0(1051) 

GO  TO  100 

250  TTT=TIHE-TSTART 

ERY=TAC*<TTT-TaU«(1.0-EXP(-TTT/TAU)») 

GO  TO  101 

END 


RAND0250 

WANDU260 

WAND0270 

RANU0260 

WAND0290 

KAND03II0 

WANU0310 

RANU032U 

KAN0U330 

KANU034  0 

WANOU  150 

WAND0J6 

K  A  N  DU  3  7 

WAN00300 

WAND0390 

LAN00400 

WAN 00410 

WAND0420 

WANU0430 

WAND0440 

WANDO450 

WAND0460 

WANDO  4  70 

WAND04B0 

WANUU490 
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Figure  16.  Subroutine  GWAND  Flow  Chart 


TABLE  XVI.  SUBROUTINE  ARB  FORTRAN  LISTING 


t  CORTHAN  DECK 

CAPM  CIIIOANCh  I  AW  DniMMCAIION 


SMIlROlJT  IME  AKU 

ARK 

0020 

cupimii  / bSAP?/ 

v  i?bo),r  c/soi.i;  (2*«> 

COHfON  / SS AM  1 / 

PFAn.DHI.AUIOI.HHI 

EUlllVAl  I  NCE 

AH') 

0060 

1  (V(  n.Al.T 

>.(V(2  I.OAC 

».iv(  Ji.nnc 

>.(V(  4 ) . 0 VC 

>. 

AIM) 

00/0 

2  (vi  s  > » o  a 

).(V(  6). IIP 

l.(V(  7  > . D  V 

),(V(  Bl.VXh 

Aim 

0080 

.1  (V(  9J.VVM 

>.  (  V(  1  ll  )  ,  VZM 

i.ivm  >,wx 

>.( Y( 12). MV 

>. 

AMI 

U 119 1) 

4  (  V  ( 1 .1 )  >  W  2 

).(V(14).AkM 

1 . 1  V  ( 1 S  ) .  A  V  M 

).(V( 16), A2M 

>. 

AIM 

0100 

*>  (V(|/>,A/CMD 

).(  VdAl.AVCMn 

>.(  vim.VAu 

)  .  (  V(  20  )  .  HOI.I. 

), 

ADI) 

0110 

A  (V(2t ).P| fCH 

> » ( VI 22  I . KXS 

) »  (  V I  2  J  > »  K  Y  S 

).(V(24).I<ZS 

>. 

ARM 

0120 

7  <  V(2s>.  IMP 

).(  V1PAI.TF.AV 

I ,  (  V  (  ? / ) , Sh  0 A 

) .  (  V(  20  )  .  SAI.A 

>. 

ADM 

01.10 

H  (V(29>.RX 

i.ivijm.NY 

).(  VI.H  ),0I  AP 

).(VI.I2).UEAY 

). 

ARM 

0  MO 

0  (  V(  .1.1 ) .  ALPHA 

)  >  (  V(  .14  1 .  ALPHA P  )  >(  V(  .ID  >.  ALPHA  V 

) .  (  V  ( .16  ) .  V  XT 

) 

ARII 

0160 

fdiii  vai.fngf 

AIM 

0  160 

1  (V(J7>.VVE 

i.ivum.vif 

I.IVMU.O 

)  .  (  V ( 4  0 )  .  VM 

>. 

ARM 

0  1/0 

2  (  VI41  > »  AM 

> . ( V ( 4? I , ACP 

I.IVM. M.ACY 

)  .  (  V(  44  ),  |IUX 

>. 

AIM 

ill  HO 

.i  (VMS i.nwr 

>.(V(1A).I)U7 

>.(  V( 4/ > .  IIUAC 

1 ,  (  v  (  4  8  )  ,  ll  ni'C 

). 

AM) 

0190 

s  (VMsi.nnvi; 

>.  (  v(*»9  >.  f  SHI  SS  i .  (  VtSI  >.  VSMISS 

>.(V(S2>. /SMISS). 

ARM 

0200 

ft  (VCj.U.MAP 

>.  (  VISA  I  ,  W  v  |> 

l.l  V(SS),W/I> 

>  .  I  VISA  > .  IIHAP 

). 

AIM 

0210 

/  (VCi/l.lHIVP 

i.(  visai.iH/p 

).  (  VlMH.AXP 

1  ,  I  VI  AH  1  ,  AVH 

). 

ARM 

0220 

H  (  V ( A l  )  .  A  /P 

>. ( V(h2), VXD 

)  »  (  V  1  A  J  )  »  V  Y  P 

),(V(A4>, V/P 

> 

ARM 

0  2.10 

E  HIM  VII  MCE 

1  ( 0  ( 1 119  ) .  I'S  1  PM  I .  ( 0  ( 1 1  0  ) .  (IP  ) .  ( C  ( 1 1 1 ) .  AK  S  1  OP  > 

ir (UFAP.f D.e.  > 

in  To  in 

P'illfl-O. 

in  r wimp-si  »;a»s7 . 

P'»*>/P*C(7> 

If  (IWHIiD.ni.U 

.in  Chucks# . I 

ADD 

02S0 

cam  mi  (ppmip.iiL.iirn 

AIM 

1126  0 

CM  1  1  Ali  (  III  1  . 

HI  .IIIASI  .HI  AS!  N 

.  Ill<  |Af.  .CIS), DUMDUM) 

AH8 

0270 

pl'KII'sA/CMD/lMH) 

C  A I  I  INK  R(  t'PSI  H,  lltIHl.  DIIM2.PSIP.PSIPN) 

p‘inin»si,sip-psip* 

ll  ii’SMi  P.u  .u. )  PSIOl  P«». 

ACI’s  A/CMD-MP*  Ak!>  I  OP*Pt»  I IICP 
AiIY-irCMD*  H2H  ) 

ah  acpmcp-iiiascmma)*  m»> 

If  (PI  All.  1:11. 0  •  )  IIU  10  80 
AOPM. 

aoy=o. 

•»1  PFTIIRM 

END 
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TEST 

READ 
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Figure  18.  Subroutine  ARB  Flow  Chart 


TABLE  XVIL  SUBROUTINE  BRF  FORTRAN  LISTING 


t 

CbrfF 


5 


10 

20 


40 

50 


FORTRAN  DECK 


BLIND  RANGE  FILTER 

HRF 

0010 

SUBROUTINE  BRF 

BRF 

.0020 

COMMON  /SSAM2/  V  (250), T  (250), C  (250) 

COMMON  /SSAM1/  READ, OfcLT.AUTOT, TIME 

EOUI VALENCE 

RRF 

0070 

1  (V(  1 )  , ALT  ) , ( V ( 2  ) ,  DAC 

),(V<  3 ) , DPC 

>,(V(  4 ) , DYC  ), 

BRF 

00B0 

2  (V(  5  ) , DA  ).(»(  4). DP 

).(V(  7 ) , D Y 

),(V(  8) . VXM  >, 

HRF 

0090 

3  (V(  9 ) »  VVN  ) , ( V ( 1 0 ) , V/M 

). ( V( 11 >.MX 

) ,  (  V ( 12 ) , M Y  >, 

RRF 

0100 

4  ( V ( 1 3 ) , MZ  ),(V(14),aXm 

>,(  VU5),AYM 

),(V(16).AZM  ), 

RRF 

0110 

5  ( V ( 1 7 ) , AZCM0  > , ( V ( 1 0  > . A  YCHb 

), ( V(1V), YAU 

) ,  (  V(  20 ) , ROLL  ), 

HRF 

0120 

6  ( V ( 21 ) , PITCH  ), ( V( 22 ) , RXS 

>. ( V(2J),RYS 

) ,  (  V  (  24  > . RZS  ), 

HPF 

0130 

7  ( V(25 ) , TEAP  ) , ( V ( 26 ) , TF A V 

).<V(?7),SE0A 

>,(  V(28),SAGA  )* 

UNF 

0140 

6  ( V ( 29 ) , RX  )»(V(3Q),RY 

). ( V(3l ),GEAP 

) ,  (  V(  32 ) , OEA  V  ). 

OKF 

0150 

9  (  V ( 33  > ,  ALPHA  > ,  (  V  ( 34  ) .  AI.PN  AP 

), ( V(3b), ALPHAY 

) ,  (  V( JA ) , VXE  ) 

RRF 

0160 

EQUIVALENCE 

DRF 

0170 

1  (V(37),VVE  ).<V(3«).VZE 

> ,  ( V ( 39  > ,  0 

>.(V(40),VM  ), 

ORF 

0180 

2  ( V(41 ), AM  ) , ( V ( 4 2 ) , ACP 

), ( V(4J),ACY 

),(V(44),DMX  >, 

OP.F 

0190 

3  ( V ( 4  S ) , DM Y  ),(V(44),D<Z 

) .  (  V  (  47  ) ,  DliAC 

).(V(4B),nOPC  ># 

RRF 

0200 

5  (V(49),nOYC  ).  (  V(50  ),  ISrllSS 

) . ( V ( 51 ) , YSM 1 SS 

), ( V(b2),ZSM|SS), 

BRF 

0210 

6  ( V(53 ) , MXP  )«(V(54).uYp 

) , ( V(55), H2P 

),(V(56).0MXP  ). 

RRF 

0220 

7  ( V(57 ) , DMVP  ),(V(58),04ZP 

)»<V(59),AXP 

>,(V(60),AVP  >, 

HRF 

0230 

8  ( V(6t ) , A2P  ), ( V( 62  ) , VXP 

), ( V ( 63  > , VYP 

) , ( V ( A4 ) , VZP  > 

BRF 

024U 

IF  (READ. EO. 0.0)  GO  TO  5 

ORF 

0250 

CUSSIO  =  COS(T  (3 ) *3 . 1  415927/1 8 0.  0  ) 

S 1  NS  1 U  =  S 1  NIT  ( 3 )«3. 1415927/ t B0. 0) 

IF  (RXS.LT  .C(9) )  .'.0  YO  10 

RRF 

0280 

CALL  DIF(V(43),I)SI0Y,DIIM?) 

HRF 

0290 

CALL  LARI  V(43).DSIOY.  XX  ill.  XX YN . UXX Y , C(5 > ,  DUM6) 

BRF 

0300 

V ( 75 )  =  V( 73  ) 

BRF 

0310 

GO  TO  20 

PHF 

0320 

AC Y=XXYN 

RRF 

0330 

IF  (RXS.Lt.C(3) )  GO  TO  40 

RRF 

0340 

CALL  OlF(V(42),DSIGP.nUMl  ) 

DRF 

0350 

CALL  LAO( Vf 42), PS  IBP. XXPQ. XXPN, 0XXP,C(5 ) , DUM5 ) 

RRF 

0360 

V(76)*-V(72)»S|NSIG‘ V(74)«CUSSIG 

RRF 

0370 

GO  TO  50 

BRF 

0380 

ACPsXXPN 

BRF 

0390 

RETURN 

BRF 

0400 

end 

RRF 

0410 
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Figure  19.  Subroutine  BRF  Block  Diagram 


Figure  20.  Subroutine  BRF  Flow  Chart 
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2.  3.  8  Autopilot  Subroutine  (MPILOT) 

(U)  This  subroutine  simulates  the  behavior  of  the  missile  autopilot.  It 
accepts  inputs  from  the  ARB  and  BRF  subroutines  as  well  as  linear  and 
angular  acceleration  components  of  the  missile.  Its  outputs  are  the  command 
control  deflections  in  yaw,  pitch  and  roll.  A  FORTRAN  listing  of  this  sub¬ 
routine  appears  in  Table  XVIII.  The  corresponding  flow  chart  and  block 
diagram  appear  in  Figure  21  and  22,  respectively.  A  glossary  of  terms 
used  in  this  subroutine  appears  in  Table  XIX. 

2.  3.  9  Flipper  Subroutine 

(U)  The  Flipper  subroutine  simply  accepts  the  outputs  of  the  Autopilot 
subroutine  and  processes  them  to  obtain  the  actual  control  surface  deflections 
in  the  three  control  axes.  The  listing,  block  diagram,  and  flow  chart  are 
shown  in  Table  XX  and  Figures  23  and  24,  respectively.  A  glossary  of  terms 
appears  in  Table  XXI. 

2.  3.  10  Aerodynamic  Subroutine 

(U)  The  Aerodynamic  (Aero)  subroutine  is  the  most  complex  of  the  entire 
program.  The  complete  set  of  aerodynamic  equations  for  forces  and  moments 
on  the  missile  (as  available  from  wind  tunnel  tests  or  actual  flight  tests)  is 
programmed.  Table  XXII  contains  the  subroutine  listing,  and  Figure  25  is 
the  block  diagram  of  this  portion  of  the  simulation.  A  flow  chart  is  shown 
in  Figure  26,  and  a  glossary  of  terms  is  contained  in  Table  XXIII. 

(U)  The  aerodynamic  data  is  readily  available  only  in  maneuver  axes; 
therefore,  the  force  and  moment  coefficients  are  generated  in  these  axes. 
Because  of  the  complexity  of  the  final  equations  of  motion,  the  aerodynamic 
coefficients  are  generated  in  a  series  of  steps  which  are  labeled  intermediate 
expressions,  secondary  expressions,  and  primary  expressions;  these  steps 
are  shown  in  sections  (A),  (B),  and  (C),  respectively,  of  Table  XXIV.  It 
is  desirable  to  integrate  the  equations  of  motion  in  missile  coordinates;  so 
the  coefficient  generates  in  maneuver  axes  are  transformed  to  missile  axes, 
as  shown  by  the  equations  in  Table  XXV.  The  final  six  equations  of  motion 
in  missile  axes  are  shown  in  Table  XXVI. 

(U)  This  subroutine  contains  numerous  parameters  which  are  input  as 
function  generators  and  are  used  in  the  generation  of  the  aerodynamic  coeffici¬ 
ents.  This  is  achieved  through  curve  fitting  techniques  applied  to  raw 
aerodynamic  data.  A  separate  list  of  the  parameters  contained  in  the  function 
generators  is  supplied  in  Table  XXVII. 

(U)  Among  the  options  available  in  this  subroutine  are  the  selection  of 
different  values  of  missile  mass  and  moments  of  inertia  during  several 
stages  of  thrusting.  Three  points  in  time  are  chosen  (corresponding  to  the 
thrust  interval)  and  corresponding  values  of  mass  and  moments  of  inertia  are 
also  selected.  During  the  portion  of  the  simulation  contained  within  this 
initial  interval,  the  mass  and  moments  of  inertia  are  varied  linearly  between 
the  selected  values. 
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TABLE  XVIII.  AUTOPILOT  SUBROUTINE  FORTRAN  LISTING 


f 

f ONTXAN  DECK 

cr  ii( 

MAVERICK  AUlOPILOr 

P 1  LOO  010 

SMRROiil  IKE  *iPIL0T 

PIL0002U 

COMMON  /SSAMl/  MEAD, DELT.AUIOI. TIME 

COMMON  /SSAH2/  V  <  250  » »  T  <  250  » .  C  ( 250  > 

EQUIVALENCE 

P 1  LOO  050 

1  <  V < 1 2  > • WY  ) , ( VI 13), W2 

),  (Vtll ).UX 

) . ( V<  20 ) , PH  1 

), 

PILO0060 

2  <VI4/)»ACP  ) ,  ( V( 1 4 ) , AX 

) , ( V ( 43 ) , AC  Y 

>, ( V(2  ) , DAC 

>. 

P 1  LOO  070 

3  ( V  <  3  ) , DPC  ) , ( V ( 4  ).DVC 

) 

piLOoneo 

F"UI VAI  ENCE 

PIL00090 

l  (  V  I  6  0  ) ,  A  Y  ) ,  (  v  (  6 1  ) ,  A  Z 

),(V<47),DDAC 

) , ( VI 46 ) , DDPC 

). 

PILOOIOO 

2  (  V  (  44  )  ,  DWX  )«(V(45),UrlY 

),( V(46),DNZ 

) , ( V( 49 ) , DD YC 

) 

PIL0U110 

Emil  VAI  ENCE 

PILO0120 

1  (  C ( 43  > » CSN  ) ,  ( C ( 44 ) ,  AK 

) .  ( C ( 45 ) ,  0J 

).(C(46),PHIJ 

>, 

PILOOIOO 

?<r.(47),ACCLIM),(Ct4H),nirLlM 

» .  (f.(4R  ),  1  YALD 

>» (C(5U  )»  TYALO 

), 

r  |  LOO  140 

3(C<51).1YHLH  ),  IC(52).TYI)LC 

1.  (C(53),IiEA0l 

> , (C( 54 ) , OMOL INI, 

PIL00150 

4 ( C ( 55  ) »  TPCLC  ),<Ct5fil.RSR 

) . <  C  <  5 7  ), TAUACC),  (C(5ft).  TAURO 

), 

DL00160 

5(0(59). TiPHIJ  ) ,  (  C  l  Ml  )  >  P  A  L  T 

) .  (C(M  ),  1R1 

),  (C(62), TR2 

), 

PILO0170 

6  (  C  ( (>3  ) ,  liALl  ),<CI«4),DA12 

).  (C<65 ),kGL 

) 

PILOOIOO 

c 

r  l  4 ;u  THROUGH  C(72)  RESERVED 

EOH  THIS  SURROMT INE 

P IL00190 

II  (NEAII.EU.O.O)  00  TO  50 

PILU0200 

CDNsinn. 0/3. 141^927 

PILO0210 

PK=CON«HJ 

PILO0220 

P'i  I  K=CON»PH  1  1 

PILO02J0 

ni’HIKsOPHI  J»CON 

P  1  LOO  24  0 

Y 1 N1 =0 . 0 

PIL00250 

PIM1-U.U 

PIL00260 

PMINsPHI 

PIL00270 

PMlOaPHl 

PILO02A0 

PPM  1  20 . 0 

P 1  LOO  290 

ISW  =  0 

PILU0300 

DAL  =  flALi 

PIL0031U 

PH|(lsT<6)/C  TRCL«»CON) 

PIL00320 

PIIN3TC6)/I IRCLCaCON) 

P ILO0330 

Y  INSOsli. 

PIL0034U 

V<XL  =  0.P 

PILO0350 

wxLiisn.u 

PILOU3AO 

TMETAX=0.0 

PIL00370 

TETAXMsO.O 

P  IL00300 

DAClsn.O 

PILO0390 

A/N=0.0 

PIL00400 

AZ1*0.0 

P  1  LOO  41 0 

P IN?=0.0 

PIL00420 

P!N20=0.0 

PIL00430 

UYL*0.0 

PIL00440 

UYLOsfl.O 

P  IL00450 

PIN3=0.0 

PIL00460 

P 1 N3O=0 • 0 

P  1  LOO  4  7 II 

A  YNsO.D 

PIL00400 

AYlsQ.fi 

PILO0490 

Y 1 N?  s  o . 0 

PIL00500 

YIN?0:D.O 

PILO051I) 

WZLsO.O 

P1LO0520 

hZLOM.I 

PIL00530 

Y  1  N3C  0 . 0 

P IL00540 

A  Z0  =  0  . 

PIL00550 

nilPINsO. 

PIL00560 

AYO  =  0. 

PIL00570 

DUYINaO. 

P  IL00500 

Y IN30*0. 

P IL00590 
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TABLE  XVIII.  AUTOPILOT  SUBROUTINE  FORTRAN  LISTING 

(CONCLUDED) 


DM1 =0. 

0N2*0 . 
pm3=o. 

DM4=0. 

0M5=0. 

0M6*0. 

0M7=0. 

OHB=0. 

pmpso. 

OHIO*#. 

0*11*0. 

DM12*g. 

GO  TO  65 

c 

C  PITCH  CONTROL 

50  (FI ABS(UX).UE.RSM)  ISM*1 

IF < <AX. UT.CSW). OR. (TIME. LT. DEAD! ) )  GO  TO  65 
60  CALL  l)|F  (AZ.OAZZ.AZO) 

CAl  L  LAG IAZ. RAZZ, aZ1.AZN.0A2, I  AO ACC. DPI) 

CALL  0  I r  < ACP.DACP.ACPO) 

CALI  LlMlTIACP.UACP.  ACCLIM.-ACCL If > 

PIN1=AK*<ACP-AZ1 I 

CALL  OIF  (PIN1.UPIN.0UPIN) 

CALL  LDLAC<P|U1,0?IN  .  P  I  n?0  ,  P  I  N2  ,  OP  I  N2  ,  T  Y  ALI) .  I T  ALG .  DM2  ) 

CAl  L  LIMI T(P|N2,0P|N?.niFLlh.-niFL IP) 

CALL  LAGlHY.PHY.HYLO.  MYL.  UMYL.  TAUNG. DP3) 

CALL  LIM!T<NYL,MMYL.rfCL.-WGl> 

CALL  L  OL  AG  (  U YL0*BK.  il  I  YL*fl*  .  P  I  N3(J  .  PI  N  J.  OP  I  N3 ,  T  Y  OLD,  T  Y  OLG  .  DMA  ) 
C 

C  YAH  CONTROL 

CALL  Olf'l  AY.PAYY.AYO) 

CAl  L  LAGlAY.nAYY.AYl.AYN.il  AY. lAUACC.DPS) 

CALL  OIF (ACY.OACY.ACYU) 

CALL  LIMI  MACY.liACY.ACCL  Irt.-ACCLIf  ) 

YIN1sAK«(AY1-ACY) 

CAl  L  DIM  YINl.pVfH.OdVIN) 

CALL  LPLAGIYIMl.nYII.  >  Y I N20,  Y  IN2.  ii  Y  IN?.  TYALIi.  TYaLG,UM6) 

CALL  LIMT(Y|N?.nYI|i?.niFLlP.-niFL  IP) 

CAl  L  LAO(UZ.nA2.WZLO.kZL.lih2L.TAUHR.DP7> 

CAl  L  UNI  UwZI.liuZL.wGl.-kGL  ) 

CAl  L  LDLAG(HZLO»Rk.l'W2L»BK.YIN30,Y  I  A3 , 0  Y  I  N3,  T  YOLO.  TYULG.  DMA  ) 
DPC=PIN20*PlN3n«YlN?0 
DnPC=0PIN2*nPIN3*UYlN2 
0YC*YIN?O*YIN30-PIN2O 
DDYC*0YIN2*DYI)F3-UPIN2 
C 

C  ROLL  CONTROL 

65  IFMTINt-.LT. PEART). AND. (ISM.Eu.ll)>  GO  TO  i 0 1) 

70  CALL  LAG(UX.nhX.HXLn.MXL.I)UXL.TAUI<r.,D)'9> 

CAl  L  LIMI  TIMXl.l'MXL.GKOL  IM.-GHOI.  l«  ) 

CALL  LARI  0.0.  O.U.PMIO.PHIN.OPUI. ThCLG.DHlO) 

WXJ  =MXL0*PHI0«T(21 ) 

DMXl *PfcXL*DPMI 

90  CAl  L  CRATEll.MXl. 0WX1.THE TAX. IE1AXN.PTAX, DM11) 

RIN1 =OPHlK«HXi*PHlK»TMETAX 
0R|N1*DPHIK*0MX1 .PHI  MPT  AX 

CAl  L  LDLAG(R|«t .DRINl.nAC.DACl ,  DP  A  C  .  IF  1 ,  TKZ.  DM12  > 

IF  (TIME. OT. PALI)  0AL-DAI.2 
CAl  L  LIMIT!  liACl. 01/ AC, UAL. -DAL* 

100  CALL  OIFIAX.PAXX.UAXXO) 

CAl  L  SPTESTI-AX.-OAXX.-GSU) 

CAl  L  TTFSTI OEAPT ) 

CAl  L  TTESTI UALT > 

RETURN 

END 


PILOOY.OO 

PILO0610 

PILO0620 

PILO0630 

PILO0640 

PILO0650 

PILO0660 

PIL00670 

PILO0680 

P I  LOO 690 

P I  LOO  70  0 

PILOU71U 

PILOD72U 

PILOfl/30 

PIL00740 

PI  LOO  750 

PILO0760 

PIL00770 

P I  LOO  780 

PILO  785 

PILO0790 

PIL0080U 

P I  LOO  81 0 

PILO0H20 

P I  LOO  83  0 

PILO0840 

PILO0H50 

PILO0860 

P  I  L  0  0  8  7  0 

PIL00880 

PILO0890 

PIL00900 

PILO  905 

PIL00910 

P  I  L  0  0  9  2  0 

PILO0930 

PILO0940 

PILO0950 

PILO0960 

PILO0970 

P  I L  0  0  9  8  0 

PILO0990 

P  I LO10  0U 

PILO1010 

PILO1020 

PILO103U 

PIL01040 

P1L01050 

PILO1U60 

PILO107U 

PILD1080 

P  ILO1090 

PILOUQO 

PILO1110 

PILO1120 

PILO1130 

PILO1140 

PILO1150 

PILO1160 

PIL0117U 

PILO1180 

PILO1190 

PILO1P00 

P  ILO1210 

PILO1220 
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MAVERICK  AUTOPILOT 
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Figure  21.  Autopilot  Subroutine  Flow  Chart 


±OIFLIM 


Figure  22.  Autopilot  Subroutine  Block  Diagram 


TABLE  XIX.  AUTOPILOT  (MPILOT)  SUBROUTINE 


Name 


V(2) 


DAC 


6 


Quantity 


Units 


Coordinate 

System 


ac,  Aileron  deflection 
command 


deg 


Missile 


V(3) 

DPC 

V(4) 

DYC 

V(ll) 

WX 

V(  1 2) 

WY 

V(13) 

WZ 

V(  14) 

AX 

V(20) 

PHI 

V  (42) 

ACP 

V  (43 ) 

ACY 

V(44) 

D  WX 

V(45) 

DWY 

V(46) 

D  WZ 

V(47) 

DDAC 

V(48) 

DDPC 

V(49) 

DDYC 

o 

> 

AY 

V(61) 

AZ 

T(2 1 ) 

C(43) 

GSW 

C(44) 

AK 

6pCi  Pitch  deflection 
command 

6yC ,  Yaw  deflection 
command 

<*>x,  Angular  velocity 

u> y,  Anguiar  velocity 

uz,  Angular  velocity 

Ax,  Propulsive  and 
aerodynamic  acceleration 

<f>,  Euler  roll  angle 

ACp,  Acceleration  com¬ 
mand  pitch 


ACy,  Acceleration  com¬ 
mand  yaw 


Scalar,  components 
of  missile  angular 
acceleration  in 
autopilot  axes 


6ac,  Aileron  command 
rate 


6pC,  Elevator  command 
rate 


,  Rudder  command  rate 

a  missile  lateral 
acceleration 
z  components 
Roll  rate  bias 


Autopilot  lateral  channel 
activation  switch  level 

Ka,  Lateral  channel  gain 


deg 

deg 

rad/ sec 
rad/ sec 
rad/sec 
g's 

rad 

g's 

g's 


rad/ sec^ 

deg/ sec 

deg/sec 

deg/ sec 

g's 

g's 

rad/ sec 
g's 

deg/g 


Missile 

Missile 

Missile 

Missile 

Missile 

Missile 


Autopilot 

Autopilot 


Missile 

Missile 

Missile 

Missile 

Missile 

Missile 

Missile 


TABLE  XIX.  AUTOPILOT  (MPILOT)  SUBROUTINE  (CONCLUDED) 


Name 

Quantity 

Units 

C(45)  BJ 

K^,  Damping  gain 

deg/deg/ sec 

BK 

K^,  Damping  gain 

deg/rad/sec 

C(46)  PHI  J 

Roll  channel  gain 

deg/deg 

PHIK 

Roll  channel  gain 

deg/rad 

C(47)  ACCLIM 

Acceleration  limit,  lateral 
channels 

g's 

C(48)  DIFLIM 

Command  limit,  lateral 
channels 

deg 

C(49)  TYALD 

T  ,  Lead  time  constant 
'sn 

sec 

C  ( 5  0 )  TYALG 

Tg(j,  Lag  time  constant 

sec 

C(5 1 )  TYDLD 

t\  ,  Lead  time  constant 

sec 

C(52)  TYDLG 

TYni  r  *  time  constant 

sec 

C(53)  DEADT 

Autopilot  activation  delay 

sec 

C(54)  GROLIM 

Roll  rate  signal  limit 

rad/ sec 

C(55)  TRCLG 

Tra 

sec 

C(56)  RSW 

Roll  rate  switch  level 

rad/sec 

C(  57)  TAUACC 

Tacc,  Lateral  channel  time 
constant 

sec 

C(  58)  TAURG 

trG,  Lateral  channel  time 
constant 

sec 

C(59)  DPHIJ 

0^.,  Roll  rate  gain 

deg/deg/ sec 

DPHIK 

0^,  Roll  rate  gain 

deg/ rad/sec 

C(60)  DALT 

Roll  channel  limit  change 
time 

sec 

C(6l)  TR1 

r  , ,  Lead  time  constant 
r  1 

sec 

C(62)  TR2 

T  £,  Lag  time  constant 

sec 

C ( 6 3 )  DALI 

Roll  command  limit 

deg 

C(64)  DAL2 

Roll  command  limit 

deg 

C(65)  WGL 

Lateral  channel  rate  limit 

_ 

rad/ sec 
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TABLE  XX.  FLIPPER  SUBROUTINE  FORTRAN  LISTING 


%  F OUTRAN  HECK 

CM  I  f  rUPPFR  WITH  THRESHOLD  PROVISION 

SIIHNOIIT  IMt*  HFL  I  P 

COHMOS  / SSAM1  /  REAO.'lELT.AUTOT.TIPE 
COHMON  / SSAM?/  V  (250). T  (250), C  (250) 

EOHI VALENCE 

)  ( V  <  2 ) • D AC  >.  (v(47),D'lAC  >.(V<  3>,0PC  ) ,  ( V ( 48 ) , DDPC  >» 

2  (V(  A ) ,  D  YC  ). (V(49),ni>YC  ) ,  (  V  (  5). DA  ).(V(  6). DP  )• 

■s  <v(  7).r»Y  ) 

F')i>|  VAII  NCE 

1  (r(97).UA|N  ).(C(98>.TAU  ) , ( C < 99 ) ,  VL I H  ) , ( C ( 1 0 0  ) , PLI H  ), 

2  ( r ( 1  n 1 ) . T HUES ) 

ir  (RrAii.EO.n.n >  on  to  10 

n  i = o .  o 

n  2  =  ii .  o 

03-11 . 0 

IM  =  U.0 

UMIsO.O 

1)112  =  0.0 

nn.iso .  0 

OilAsO.il 
I  M  = 1  .fl/OAIN 

in  II  i  THAfS.EO.ll. 0 )  GO  ID  15 
r » =  < n vr-HAC-ru  )*ga I n 
nri  =  ( 'MlYC-HDAC-l'D)  )»GA  IN 
r  2  r  ( It  AC  ♦  !ir*C- nss  )»GAI  M 
III  2  =  (  mi  AC  *  DP  PC  - 1*1)2  )  •OA  IN 
r.isi  PYr*DAc-ru)*GAiN 
nr.3=(  IIIIYC  »n!JAC-DD3)»GAIH 
r')*(iipoii«c-n4)*GAiN 
HI  As(  nni'C-niUC-DDA  )»OAIN 
r.o  id  ?u 
15  risHYC-DAC 
nn  slUlYC-HDAC 
rpniAonrc 
012  =  III)  AC  *  DDPC 
F  3  s  li  y  c  ♦  II A  C 
niMsnnvc*  ddac 
fasopc-uac 
nr«i  =  iinpc-DDAi: 

CAI  1  L  A  l!  (  F 1  .liFl.'l'in.Ol.nm.TAl.OUl) 

CAM  VLIHITiniO.Ol.lllil.VUH.-VUH) 

F  Al  I  Lin  T(H1.  DDE.  PI.IM.-PLIR) 

CM  l  LAC(F2,IiF2, D20, H2.HII2.TA1. 002) 

CALI  VI  |l‘IT(l)2n.n2.Pb2*  VLIP.-VLIM) 

CAM.  Lin  T(li?.ni'2.Pl .  Ih.-PLIP) 

CAM  LAG  (  r  j,  OF 3  .  H3 (i .  D3.  Uli3,  T A 1  >  DU3  ) 

CAI  I  VLIM|T(!i30»D.»,Dr3.  VI  IH.-VLIN) 

CAM  Lin  f(l'3»nl‘3.PL|H,-PLIP) 

CAM  l  AG(F4,br4,.)40,ll4.PD4.TAJ,DU4> 

CAM  VI  I  M  I  TIDAO.  DA.riPA.  VLIH.-VLIM) 
cai  i  Lin  r ( im. dm. pl i h. -pl no 
C'1  10  30 

2  n  CAM  LAGCFl.DI  1.710. 71.  DZ.TAU.Dl'l) 

CAI  I  DHAHIK/III.DZ.Z/I.TMPES.-IMKES) 

CAM  LF'tl  MZZl.t  Zl.  VLIM.-VLIH) 

cmi  oRATLd.zzi.nzi.Dio.ni.ooi.nu?) 

CAI  L  LIHl  f  (Dl.DDl.F'UH.-PLIH) 

CAI  I  L AG (F2.  nr 2.  Z20.Z2.DZ.T All. DD3> 

CAM  DBANIM  Z20, DZ.Z/2. T HRfrS, - IHPES  I , 


FLIP0010 

rL!P0020 


rLIP0050 

FLIPOOISO 

FLIP0070 

FI  IPOOfll) 

FL  IP0II90 

FLIP0100 

FLIP011U 

FI IP0120 

FL  IP0130 

ri.  IP014D 

FL  I PO 150 

FL  1 PU 1 60 

FL  IP017U 

FL  IP0180 

K  IPD190 

F  L  I  P  0  2  0  0 

FLIP0210 

FL  IP0220 

FL1PU230 

FL  IP0240 

FLIP0250 

FL  IP0260 

FLIP027U 

EL1P02H0 

FI.  IPU290 

FLIP0300 

FI  IPOJIO 

rLIP0320 

FLIP0330 

FLIP0340 

FLIP0350 

FLIP0360 

FLIP03/0 

FLIP0380 

FLIP0390 

FL  IP0400 

FLIP0410 

FLIP0420 

ri IP0430 

FL I PO  4  4  0 

F  L l PO  450 

FLIP0460 

FLIP0470 

f LIP0480 

FLIP0490 

FLIP0500 

I LIP0510 

FI IP0520 

mposso 

FL IP0540 

FL  I  P055II 

TLIP0560 

FLIP0570 

FLIP058O 

FLIP059# 


76 


T>  B LE  XX.  FLIPPER  SUBROUTINE  FORTRAN  LISTING  (CONCLUDED) 


CAI.I  UMl  T<7Z2.1'2.  VUN.-Vim 
CAM.  GP  A  TF  d.  Z  7?.  JZ.  020.02.0  02.004) 
C Al  t  LIMIT  (D2.li02.FUH.-PI  IN) 
call  LAiura.nr 3. 7.10, z;», »/, tau. 0115 i 

CALI  ON  AND  (LTD.  I'?.  Z7  3.  THRES.  -  IHfiFd 
TAIL  LIHl K2Z3.L7.Vl IN.-VLlH) 

TAIL  GRATE d .  Z7J.  UZ,  DJD.  U3.  D0J.IIU6) 
CAU  LIPlKD3.Dli3.Pl  IM.-PLIH) 

CALI  LAC(F4.0F4.74O,Z4,l)Z,TAU.0U7) 
CA(  I  DS  AND  (  2  4  U.  ilZ  >  ZZ4 .  TMRES.-THRES) 
CAIL  I.IMI  KZZ4.IIZ.  VI  IM.-VLIM 
TAU  GR A  TF.d, 274,  Al.  04  0.04.  D04. 0118) 
CALI  LIMITCD4.0U4.PI.IM.-PLI)') 

30  DA  =  .?!>»(D2-D4-01*U3) 

OI>s.5»(D2*04) 

11 1  * , 5»  ( I'l  ♦  D3  ) 

RFTURN 

END 


FLIP0600 
F  L  IP0C.1O 
fl  IP0620 
r L  I P0630 
FLIP0640 
FI IP0650 
FI  I P0660 
F 1. 1  P  0  67  0 
FL IP068U 
FLIP0690 
FL I P0700 
F  1. 1  P  l)  7 1 0 
FLIPO/20 
Fl IP0730 
Fl.  IPO  74  0 
FLIP075O 
FLIP0760 
FLIP0770 
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FLIPPER  WITH  THRESHOLD  PROVISION 
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Figure  24.  Flipper  Subroutine  Flow  Chart 


TABLE  XXI.  FLIPPER  (MFLIP)  SUBROUTINE 


Name 

Quantity 

Units 

Coordinate 

System 

V(2) 

DAC 

6ac>  Aileron  deflection 
command 

deg 

Missile 

V(3) 

DPC 

6pC ,  Pitch  deflection 
command 

deg 

Missile 

V(4) 

DYC 

6yC,  Yaw  deflection 
command 

deg 

Missile 

V(5) 

DA 

6a,  Aileron  deflection 

deg 

Missile 

V(6) 

DP 

6p,  Pitch  deflection 

deg 

Missile 

V(7) 

DY 

6y,  Yaw  deflection 

deg 

Missile 

V(47) 

DDAC 

6ac »  Aileron  command 
rate 

deg/ sec 

Missile 

V(48) 

DDPC 

6pc»  Elevator  command 
rate 

deg/sec 

Missile 

V(49) 

DDYC 

6yC»  Rudder  command 
rate 

deg/sec 

Missile 

C(97) 

GAIN 

Kfc,  Servo  velocity  gain 

1  /sec 

C(98) 

TAU 

t§,  Control  surface 
time  constant 

sec 

C(99) 

VLIM 

&L»  Control  surface 
velocity  limit 

deg/sec 

C(100) 

PLIM 

&L,  Control  surface 
angle  limit 

deg 

C(101) 

THRES 

6th  Control  surface 
rate  threshold 

deg/ sec 

D 1 

6j ,  Control  surface 

No .  1  deflection 

deg 

D2 

^2»  Control  surface 

No .  2  deflection 

deg 

D3 

63,  Control  surface 

No .  3  deflection 

deg 

D4 

64,  Control  surface 

No .  4  deflection 

deg 

80 


o  o 


TABLE  XXII.  AERO  SUBROUTINE  FORTRAN  LISTING 


%  FORTRAN  DECK 

OAERft  MAVIRIC*  AERO 

SoRkOilT  INE  ('AERO 

COMMON  /SSAM1/  REAn.DFIT.AU10l.TtN 
COMMON  /SSAM2/  V  (250), T  <25!!),C  (25U) 

equivalence 

1  <  V  t  L  > . RLLAL f ) »  <  V  <  51, DA  ).(V(  6). OP  ).<V(  7  ) . DV  ). 

2  <  V  t  81.VXH  ).(V(  Vl.VYh  ).  (Vtlll).VZM  >,(V(111,WX  ), 

3  ( V(17  ),bV  >. ( V< 1 3  I, WZ  ) > ( V ( 1 A  ) . AX  ) , ( V< 15). AY  ), 

A  (V(IA).AZ  ).(  Vt.%3).  ALPHA  ).(  V  (  34  ).  AlPrtAP ),(  V  (  35 ),  ALPHA  V  > 

EOUI VALENCE 

2  ( V  ( 4*i )  >  DWY  ).  (  0(46), DWZ  ).<V<39>,0  >,(V(«n),VM  ). 

3  ( Vt  41  ) . AH  ). ( V(44 >, UOX  ). (V(«4>,CAPLAM> 

foui valence 

1(0(731. S  ).(C(74),n  ). (C( 75  l.FSL  1 . ( C < 76 ) . TUOOST >  . 

2(0(77  ),  TSuST  1.  (C(7M,Fl.Til  > , ( C ( 79  1 , XBAR  >,(C(Hu),CLP  1. 

3(0(81), AE  >. (C<8?).AJX0  >.  (C(<<3),AJX1  1 ,  <  C  (  84  ) »  A  JXT  ), 

4(0(85). AJYO  1.  <C(f>M,AJU  ) . ( C ( U7 > , A JY T  ) ,  (  C  (  88  >  ,  AMASSU  )  » 

5(0(84 >. AMASS1 ).(C»9« )» AMASS! > 
o, (C<91 ), TSEOArt). (C(92), TGI  AL T) 

DATA  RD2/114.59155/ 

UAT»  KTOD/57. 29578/ 

tr  (RFAO.EO.n.n >  on  ro  5 
E  LIW  =  1. 4142  L  3/ <.»?./«57. 24578) 

sn  *  s*d 

(5=32.2 
0MV2=n/2.0 
RXB=RH2»(XHAR/D>**2 
AF(IS  =  AF/S 

C*)NALF  =  ?.0»XflAR»»2/(  3.(l»U»*2  ) 

1101  Ills  II /EL  Til 

PAJX1  =  (  AJX0-AJX1  1/TK.lOST 
nAjYi  =  (AJrn-AJyi  >/TbO«ST 
noASS=( AMaSSD-AMASSI )/THnoST 
T'm*TSUS!-liniOST 
PAJX7=(A.)X1-A  JXI  )/TI'ir 
l)AJY?=(  AJY1-A.IY1  )/T|ltr 

n  ha  ss?  =  (  amass  j-a<4  ass  n/rmr 

5  VYZ2  =  VVH»*2»Vi'H»*2 
VV7=SORT(VY22) 

VS0=VY77» VXM»»2  ' 

VM=S0HI( VSO) 

T  I  ME* T I M* TSLPAR 
A I  T  =  RFL  AL  I  ♦ TO  I Al  T 
C 

0  MACH  NUMBER 

CALL  rcrNK ll.Al T.SVEL.-l 1 
AM=VH/SVEL 

AEHO'iYNAMTC  ROIL  ANGLE.  I’M  I A 
IF  (VY7.NE.U.U)  G'l  TO  10 
CPH|A=. 70711 
SPIIIA*.  70711 
r.n  to  ?u 

10  C»,MIA*VZH/VYZ 
SPH I  A* VYM/VYZ 
C 

C  AKOLF  (IE  ATTACH,  alpha 

20  Al PHAsATAN(VY//VXM)»RTOO 
ALPHA2  =  ALI’MA»*2 


AFR00010 
ALR00020 
AERO  30 

AER00050 

AERU0060 

AFK00070 

ALR00080 

AERUUII90 

AERU0100 

AEROOUO 

AER00120 

AEH00130 

AL  ROO 14 0 

AIR 00150 

ALRU0160 

AL-R001  70 

At  KUO  180 

ACRUU190 

ALR00200 

AFRUU210 

Al.  R00220 

AtR00230 

ALRU0240 

AFRU0250 

AFRO 0260 

AERO02/U 

AFROnZflO 

AFRU0290 

AL  RU03P0 

Al R00310 

A  I  RUG  32 0 

ALRUU330 

AL  RUG 3 40 

AFR00350 

At  RU0360 

AtR0037U 

AFROn.180 

AFROU39D 

At HU04U0 

AFRU0410 

AFRO0420 

AFRUU425 

AFRU0430 

AIRU0440 

ALRUU45 

AFRU0460 

Al RU0470 

Al R00480 

At  R0049U 

AFROUSO  0 

AI-R00510 

ALRO0520 

AER00530 

ALRUA540 

AFRU0550 

AFR00560 

Al RU0570 

AFR00580 
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TABLE  XXII.  AERO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


AI.Pma3  =  ALPMA2*ALPHA 
ALPMAP*  ALPMA*CPli  I A 
ALP('AY*aLPHA«SPH|  A 
CALI  n!FULPHAP,DALF«».DALFPl» 

CALL  DIFULPIWY.DALFY.OALFYl) 

C 

C  rilPPEK  UEFLlCTIUNS  in  maneuver  axes 
0T=0P»CPHIA-DY«SPH1A 
PR  =  l,Y*C»’HIA*OP*SPM|  A 
c 

C  At  XII.  I  ARY  FUNCTIONS 

COS4PH*l .0-b.Uti:PNlA**2*SPMI  A**2 
SIN4PM=b.l8*(ApS(CPHIA)-ARS(SPHlA>)*CPHIA*SPH|A 
l)ll2V  =  nOV2/VH 
C 

C  MANEUVER  A  XFS  AERODYNAMIC  COEFFICIENTS 

r 

c  intermediate  expressions 
c  c  sur»  small  mo 

N  =  0 

CALI  FCFNK  12. AM. SMI, N) 

CALL  ECENK  I3.AM.SM2.N) 

CAI  L  FOCNl(  I4.AM,SH3,N) 

CA1L  FOENK  19, AM, SNA, N> 

CSMO=SMl*ALPHA*(SM2*$h3*COS4PH)*ALPhA2  ♦  SNA* ALPHAS 
C 

C  DELTA  C  SUM  SMALL  M 

CALL  FCEMI  I6.AM.SM5.N) 

C  A I  L  FCENK  I7.AM.SM6.N, 

CALL  FGENK IA.AN.SN7.N) 

DCSM=  (SM5+ISM6  ♦  SM7«COSAPM)»ALPMA2)»OT 
C 

c  c  sun  shall  no 

CALL  FOENK  I9.AM.SN1.N) 

CALL  rCEHK  I1M.AN.SN2.N) 

CSN«=(SNl*ALPMA2  ♦  Sn2*ALPHaJ)*SIN4PH 
C 

0  DELTA  C  SUU  SMALL  N 

CALL  FGEMK  IU.AM.SN3.N) 

CALL  FGENK  I12,A,N»SN4.N> 

CALL  FGHN1L llJ, AM. SN9.N) 

CALI  FCENK  m,AM.SN6,N> 

CALL  FGENK  119. AM. SN7.NJ 
CALI  FOENK  lift.  AM, SMB. N) 

OCSNs ( SMS  ♦  (SNA  ♦  SN9*COS4PH)*ALPHA2)*DR 
1  ♦  ( ( SN6  ♦  SN7«COSAPH)*ALPHA  ♦  SN8* ALPHA?* DT >*DA 
C 

C  C  SUB  NO 

CALL  ECENK  117. AM. CNt.N) 

CALL  FGENK  118. AM. CN2.N) 

CALL  FOFNl( IIV.Ah.CNS.N) 

CALL  rnr.NK  120.  AM.CNA.N) 

CN0=CN1*ALPHA  ♦  (CN?*CN3*COS4PH)*ALPHA2  *  CR4*ALPHA3 
C 

C  DELTA  C  $UH  N 

DCN=OOLTB*DCSM 

C 

C  C  SUB  VO 

CALL  FOENK  121,  AH.  Yl.N) 

CALL  rCfcNIC I22.AN.V2.N) 


AERU0590 

AFRU0600 

ALROOMO 

AER0U620 

AER00630 

AEROOAAO 

AERO069O 

AFRO0860 

AERU067Q 

AERO0A8O 

AER00690 

Af  R00700 

Af  R00710 

ACR00720 

AER00730 

AER007AU 

AF  RO0750 

AER00760 

AERO0770 

AERU0/80 

AERO 079 

AEROOBO 

AkROOBl 

AERU0B2 

AFR09830 

AEROOBAO 

AERO0H50 

AERO0S6 

AFR0087 

ALRU0B8 

AERO0B90 

AERU0900. 

AERU0910 

AERU092 

AERO093 

AEROU9AO 

AER00990 

AERO0960 

AER0097 

AERU098 

AERO099 

AFRO100 

AER0101 

AERO1020 

AER01030 

AEROIOAO 

AER01050 

ALR01060 

A FRO 10 70 

AL  RU1UB0 

AERO1U90 

AfcROllOO 

ALRO1110 

AER01120 

ALROUSO 

AERUUAO 

Al R01150 

AER0116I 

AFRO 1170 

ALROUSO 
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TABLE  XXII.  AERO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


CYOM  YMALPFIA2  ♦  Y2*ALPFFA3MSINAPH 

OFLIA  C  SUR  V 

QCY=n0LT8<nCSN 

c  sijr  ro 

CAM  FGbNl< I2J.A1.CDU.N) 

CALL  EGENM  I2A.AH.C2.NF 
CAI  L  FGFNM  I25.AM.C3.N) 
cail  roFNic i2;,AH.cnrt.H) 

CAM  rr.l  Nl(  I2«,  AM.CIT2.N) 

I r  (  (TML.GT.TSUST  ) . 4  Ml  .  ( RE A D . t U , U 
CALL  EGFN1M26,AM.CDH,NF 
GO  YO  AO 
3n  COH=0.0 


An  CI»r  =  ALY*(CDU  ♦  ALMC0F2F 

ci=cdh  ♦  cnf  ♦  cdiuaeos 
CCn  =  Ol*<:2«ALPHA*CJ»ALPHA2 

PFLfA  C  SIJU  C 

CAM  FGFMM30.AH.CA.NF 
CAI  I.  FfiENU  I31.AM.C5.NF 


ncc=ci»< ALPMAP»UP-ALPMAY«0Y>»C5»DI*«2 


SLR  LO 

CAM  FGFNM  132,  A».SLl.N> 

CALL  FGFNM  I  3  J .  AM .  SL  2  .  N  ) 

CALI  FGI.NM  I3A.AH.SL3.NF 
CALL  FOF.M  (  I35.AH.SLA.NF 
CAM  FGFMt  I3o,AH.$L5,NF 
CALL  FGI.NM  137,  AM. SLA,  N) 

CALI  FGFNM  I3«, AM, SL7.N) 

CAM  FGl  NM  139,  AM. SL8.NF 
CALL  FGFNM  I A  U  ,  AM.SL9.N) 

CAM.  FGfeNMMU,  AM, SLtll, NF 

C  L  0  *  <  S  L 1  ♦SL2*ALPmA  ♦  SL3«Al.PHA2  ♦  SLA»ALPHA3)»ALPHA2»SINAPH 


r* F I  IA  C  Sim  L 

IJCI  =(  SL'jF  (SI.6*SL7»C0SAPH)«ALPHA 
1*A|.  F’HA?)*bR 
IF  <  (  1  IME.GT.TSOSIF. ANfl. (READ. EQ.O. OFF  GO  To  50 
IF  (TIPF.r.T.YHOOSI)  GO  TO  A5 


AMASS  =  AMA5S0-0MASS*T I  ME 
AJY  =  AjYli*llAJYl»l  I  ME 
AJX  =  A.|X0-l)AJ<l*TlME 
GO  TO  A A 

A 5  T  H I F  =  T I  ME- THOUST 

AMASSsAMASS1-UMASS2»TDIF 
A  JY  =  AJYl-DA.IY2*TniF 
AJX  =  A.lXl-DAJX2»TniF 
A*  AJXYM  AJT-AJXF/AJY 
C  J.O  -  PEL T A  L/LTrt 

CALI  FGfcNM  IA1.1  IMF, PLOD. -1> 

CALL  FGFNM  142, 1ME.OPDLOL.-M 
CAI  L  FGFNM  I A  A ,  T  |  ME ,  THRUST. -1> 
CAM  FGLNM  IA5.ALT.PRfcS.-l) 

TOM=< 1HHUST  ♦  AF* <PSl -PRES F  F/ AMASS 
GO  TO  60 
50  TOMsO.O 

An  CALI  FGFNM  IA3, ALT. MHO. -1) 


AI.RO1190 
AFRU1200 
AFR012T  0 
AFR01220 
AERO1230 
ALR012AO 
AFRU1250 
ALRU126 
ALR0127 
AFR0120 
AFR0129 
AFR013QO 
AFH01310 
AFRU1320 
AF.RU1330 
AFRU13A0 
AFR01350 
ALRO1360 
AER01370 
AFR013AO 
AFRU1390 
AFROIAOO 
AEROLA10 
AfcROl A20 
AFPOl A30 
AL-R01AAU 
AL-R01A50 
At  Rill  A60 
AF RO1A70 
AFKO1A80 
AFR01A90 
AFR01500- 
AFR01510 
At  ROl 52  0 
ALRU153 
AFR015AO 
AFRO1550 
ALR0156# 

)«PAMSLA»ALPllAMSLV*SLin«COSAPH)  A I  R0157  0 

At RU15G0 
AfcR01590 
AF.R01AP0 
At  ROlAlO 
AF  R01620 
AFROJ630 
AFROIMO 
AF.  RU1A50 
AFR01660 
AFPO1670 
AF  RO 1 6  8  0 
At RU1690 
At  R01710 
AFRU1700 
ALK01720 
AFRU1 730 
ALRU17A0 
AER01750 
At  ROl 766 
AERU1770 
AF.R01780 


U ) ) GO  TO  30 
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TABLE  XXII.  AERO  SUBROUTINE  FORTRAN  LISTING  (CONTINUED) 


c 

r  pki*«apy  f xpressiohs 

rstii-M=csHo-r)Loo«cwo*nKniOL*iicsH 
rsi  CN  =  CSNO*liLOn»C»0*0H0Ll)L*I.CSN  .  • 

CStll<N:rN0*0L‘N 

rsimr=rYO-ncy 
csi'nc*cco*occ 
rsi'hi.scio  ♦  not. 

CNnALP  =  CM*  ALPHA  «<Cf,2*CN3»CT»$4PH*ALPMA»CN4> 
rhPrRXft»CNOALP 
rM0ALr  =  !)fl2V»CiJMALf»CM0ALP 
C 

C  ' 

0  AFNOBYNAMIC  COEFFICIENTS  IN  MISSILE  AXES 
CS|i|!HHsCSUBH»CPH|A*CSI  CN»SPHIA 
1  ♦  CH'JAI.  F»0ALPP 
CSI .C*i|3i;SLCN«CP»<l  A-CS0BH«SPNIA 
i  -  cpn ai  f  •  OAi.f  y 
CliliPN«  =  CSURN«CPH|A  -C5U6Y»SPHIA 
CSHhY«  =  CSilBY«CPh|  A  <  CHIIHN*  SPH  I A 
C  DYNAMIC  PRESSURE#  0 
«=PHO*VSO/2.0 
C 

c:  ACCELERATIONS 

CONI  =  o«S/AMASS 
COM2  =  0*SD/AJX 
COM3  =  0»SD/AJY 
C0Y=-n02V*CLP 
coy^s-iio2v«i;ho 
TAItxsl.n/(CuX«CON2) 

TAIIZ  sl.0/(i:nYZ«CON3) 

TS=WX»AJXY«TAUZ 
AXMsCOM»CSOBC  ♦  TOM 
aYP  =  C')M*CSO0YB 
A/M  =  C0H1  *CSIIBNB 

CALL  CUL THN ( 0. -1,0.0. H. 0 . 0. G X# C Y , CZ. YAH# ROLL*  P I f CH ) 
CALL  R45F(0X»*JV»fiZ«t»XM#6VM#(»ZH) 

C 

TX=csnnL/cox 

CALL  II I E  (  TX.  DTX#/IX1> 

ty=csobmb/coyz*ts«hz 

CALL  i)IF<TY.I)fY.DYl> 

TZ=CSLCNB/C0YZ-1S*MY 
CAIL  nir(TZ.OTZ.DZl) 

IF (RFAO  .EO.  U.OJOO  TO  70 

TX  =  WX 

TY  =  MV 

TZ  *  HZ 

RTX  =  0.0 

0TY  *  0.0 

nT7  *  0.0 

70  0VX=AXM ’ OXM*HZ«VYM-MY«VZN 
nVY=AYM*CYM*MX« VZM-MZ»VXM 
nvZsAZM*GZM*HY«VXM-MX*VYM 
CALL  INrER<l>VX#XXl,XX2,VXO.VXM> 

CALL  INTERIPVT. V Y l » Y Y2 , V YO. V YM ) 

CALI  |NTER(tiyZ,ZZl#ZZ2.VZO.VZM) 

CALL  LAO  (TX.IITX.HXO.k'X.OHX.TAUX.OUPX) 

CALI  LAO  (TY.UTY.*iro.HY,nMY.TAU7,DUHYJ 
CALL  LAO  (TZ,UTZ,«ZO.HZ#DnZ,TAUZ.DUMZ> 


AERO 1790 

At R01B0  0 

AER01B10 

AFR01820 

AER01B30 

AER01840 

AFRU1850 

ALRU1860 

AERO1B70 

AFR01880 

AER01090 

AER01900 

AERO1910 

ARR01920 

ALR 01930 

ACRU1940 

A  L  R  0 1 9 1>  0 

AER01960 

AERO1970 

AER01980 

AFHU1990 

AEKU2U00 

ALRU201 0 

AERU2020 

AFRU2P30 

A»  R02M0 

AERO2050 

AFRU2060 

ALRU2070 

AERU20B0 

AER02U90 

ALRU21 0 0 • 

AER02110 

AFRU2120 

AERO2130 

AER02140 

ALRU2150 

AERU2160 

ACRO2170 

AERO21B0 

AER02190 

AFR0220  0 

AER02210 

AER02220 

AERU2230 

AFR02230 

A t  RU224  0 

ARR02250 

AER02270 

AER02280 

AERO2290 

AER02300 

AER02310 

At R02320 

AER02330 

AER02340 

ACRU235I 

At R02360 

AEH02370 

AER02380 
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TABLE  XXII.  AERO  SUBROUTINE  FORTRAN  LISTING  (CONCLUDED) 


A*=AXH/G 
A YsAVH/fl 
AZsAZrt/0 

ZAP1*C(44)»C<14>»1. 41421 J 
2AP2sC<44)*V(40>»n.UB 
t AP3=2. • AMASS/ (RM0»V( 40 )«S) 


AFR02390 

AFR0240U 

ALRU2410 

AFR02420 

AIR024JO 

AFRU244U 


ZAP4*  (Cm*(CN2*CN4»C0S4Pli)«  ALPHA  »CM*ALPHf2)/<  SHI  ♦  <  SM2*SM3»CUS4PH  )  AI  R0245U 
1  •ALPHA*SH4*ALPHA2)-D0LT8  AFRO2460 

ZAP5sSM5»(SM6*SP;»C0S4PII)»ALPHA2  AFRO2470 

CAPIAM  *  ZAP1/(ZaP2*C(45)-ZAP3/<ZAP4»ZAP5«5?.29SZ»7))  AFRO2460 

RFTURN  ALRO2490 

END  AFR02i>00 
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SIN4PH 


C0S4PH 


VM 


ALPHAP 


ALPHAY 


ALPHA 


AM 


ALT 


Q 


AM  TIME  DA  DP  r 


I _ 


Figure  L 
B1 


(The  revers 


AM 


F  UNCTION 
GE  Nt  HAi  ORS| 
I?  TO  140 
K40 


TIME 


DA  DP  DV 


ALT 


TIME 


FUNCTION^ 
jENERATORT 
141.  14? 


l0S4Ph 


ALPHA 


alphap 


AlPHAV 


GENERATE  AERO 
COEFFICIENTS 
IN 

MANEUVER 

AXES 


IV  X 


w  z 


V  M 


dalfp 


dalf  V 


z 

Z 

> 

<J 

ai 

CO 

D 

D 

in 

1/5 

V) 

O 

'  -  ' 

u 

u 

t  1 

RESOLVE 

THROUGH 


u 

GENERATE 

MASS. 

PRESS 

T/m 

CSUBC 


CSUBVB 


CSUBNB 


±±±± 


VXM 

VYM 

VZM 

WX 

WY 

WZ 


GENERATE 

LINEAR 

ACCELERATIONS 

IN 

MISSILE  AXES 


DVX 


VXM 


DVY 


VYM 


DVZ 


r 


GENERATE 

ANGULAR 

ACCELER 

IN 

MISSILE 

AXES 


TIME 


GENERATE 

MOMENTS 

OF 

INERTIA 


_ 

i 

i 

i 

qsd2 

*?J  F  H) 

X 

i 

i 

uni2 

JTV  Cm., 

L- 

1 

i 

c 

rrv  ^ 

• 

VZM 


WX 


WY 


WZ 


NOTE  IT  IS  COMPUTATIONALLY  EXPEDIENT  TO 
SOLVE  FOR  THE  MISSILE  ANGULAR  RATES 
AS  OUTPUTS  OF  FIRST  ORDER  LAGS  RATHER 
THAN  BY  DIRECT  INTEGRATION  FOR  THIS 
REASON  THE  FEEDBACK  DAMPING  IS 
INTRODUCEDSEPARATELY  FROM  OTHER 
ELEMENTS  OF  THE  ANGULAR  ACCELERATION 
EXPRESSIONS,  ANDTHE  EQUIVALENT  LAGS 
ARE  PROGRAMMED  IN  THE  SIMULATION 


Figure  25.  Aero  Subroutine 
Block  Diagram 
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(The  reverse  of  this  page  is  blank.  ) 


I 


vwvttrCt  *$*0 


>  T0OQC'T 


£ 


AMASS  AMASSO-OMASS*TIME 
AJY  AJYO-DAJYl'TIME 
AJX  AJX0-DAJX1  'TIME 


TDIF  TtME-TBOOST 
AMASS  AMASS  1-DMASS2"  TDIF 
AJY  AJY1-DAJY2*TDIF 
AJX  AJX1-DAJX2*TDIF 


AJXY  (AJY-AJXl/AJY 


CALL  FG£N1(I41,TIME,DL0D,-1) 

CALL  FCEN1(I42,TIME,OMDLOL,-1) 
CALL  FGEN1II44, TIME, THRUST, -1) 
CALL  FGENKI45, ALT, PRES, -1) 

TOM  <THRUST*AE'(PSL-PRES»/AMASS 


CALL  FGEN1(I43,ALT,RH0,-1' 


CSUBM  CSMO-DLOD»CNO‘OMOLOL‘DCSM 

CSLCN  CSNO»DLOD*CYO«OMDLOL*DCSN 

CSUBN  CNO-DCN 

CSUBY  CYO-DCY 

CSUBC  CCO*DCC 

CSUBL  CLO»OCL 

CN0ALP=CN1*ALPHA*(CN2*CN3*C0S4PH‘ALPHA‘CN4' 

cmq-rxb*cnoalp 

CMDALF  -  002V*C0NALF'CN0ALP 


CSUBMB  CSUBM'CPHIA  -CSLCN-SPHIA 
•CMDALF'DALFP 

CSLCNB-CSLCM-CPHIA-CSUBM'SPMIA 

-CMOALF'DALFY 

CSUBNB  CSUBN‘CPHIA-CSUBY*$PH!A 
CSUBYB  =  CSL'BY*CPHIA  -CSUBN’SPHIA 


Q  RHO‘VSQ/2.0 


CONI  Q‘S  AMASS 
CON2  Q‘S0  AJX 
CON3  Q*SD  AJY 
COX  -D02V"CLP 
COYZ  -D02V-CMQ 


TAUX  1  0/(C0X’C0N2) 

TAUZ  i  u/TCOYt ‘COM3! 

TS  WX'AJXY’TAUZ 
AXM  CONI  ‘CSUBC  •  TOM 
AYM  CONl’CSUBYB 
AZM  CONl'CSUBNB 

CALL  EULTRN(0,-1  O.O.O.O.G.GX  GY  GZ  YAA  ROLL, PITCH!] 
CALL  R4SF(GX,GY,GZ  GXM.GYM  GZM 


TX  CSUBL  COX 
|  CALL  DIF(TX.DTX.DXl) 

TY  CSUBMB, 'COYZ  *  TS*AZ 
|  CALL  DIFfTY  OTY.DYl) 

TZ  CSLCNB'COYZ-TS'TVY 
|  CALL  OIFITZ  DTZ.DZ1) 


DVX  A  X  U  •  G  X  M  •  VVZ  ‘  VY  M  -  AY  •  VZ  M 

DVY  -AY  M  •  CY  M  •  AX  •  VZ  M- VVZ  ‘  VX  M 

DVZ=AZM‘GZM*AY'VXM-AX*VYM 

CALL  INTER  DVX, XXI, XX2  VXC.VXM 

CALL  INTER'0VY,YY1,YY2.VY0,VYM 

CALL  INTER  DVZ,ZZ1,ZZ2.VZ0,VZM 

CALL  LAC:TX,DTX,AXO,WX  DAX.TAUX  DUMX) 

CALL  LACITY  DTV.AYO.AY  DAY ,  TAUZ ,  OUMY) 

CALL  LAG!TZ,DTZ,AZ0,AZ,DAZ,TAUZ-DUMZ; 

AX-AXM/G 

AY  AY  M/C 

AZ  AZM  G 

ZAP1  X!44)*C'14)*1. 414213 

ZAP2-C(44)‘V'40)‘FLUB 

ZAP3  -•2,‘AMASS/TRH0‘V(40)*S) 

ZAP4  ICNl-tCN2-CN3‘COS4PH) 
■ALPHA"CN4’ALPHA2)/ 

(SMI  -ISM2  •  SM3‘C0S4PH' 
■ALPHA-SM4-ALPHA2  -OOLTB 
ZAP5  =  SM5-!SMb*SM7‘COS4PH  -ALPHA 2 
CAPLAM  ZAPl/(ZAP2-C(45)-ZAPy' 

(ZAP4-ZAP5‘S7.29S787'I  _ 


f  return) 


Figure 


26.  Aero  Subroutine 
Flow  Chart 
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TABLE  XXIII.  AERODYNAMICS  (MAERO)  SUBROUTINE 


Coordinate 

Name 

Quantity 

Units 

System 

V(l) 

RELALT 

hm » 

Missile  altitude  above 

ft 

Inertial 

ground 

V(5) 

DA 

6a » 

Aileron  deflection 

deg 

Missile 

V(6) 

DP 

6p » 

Pitch  deflection 

deg 

Mis  sile 

V(7) 

DY 

6y  > 

Yaw  deflection 

deg 

Mis  sile 

V(8) 

VXM 

V*. 

Velocity  X-axis  missile 

ft/ sec 

Missile 

V(9) 

VYM 

v 

Velocity  Y-axis  missile 

ft/ sec 

Mis  sile 

V(10) 

VZM 

Vz’ 

Velocity  Z-axis  missile 

ft/ sec 

Missile 

V(ll) 

wx 

wx 

Components  of  angular 

V(12) 

WY 

CO 

y 

velocity  in  missile 

rad/ sec 

Mis  sile 

V(  13) 

WZ 

wz 

axes 

V(14) 

AX 

Ax 

Propulsive  and 

V(  15) 

AY 

Ay 

ae  rodynamic 

g 

Missile 

V(16) 

AZ 

y 

Az 

acceleration 

V(33) 

ALPHA 

a,  Total  missile  angle  of 

deg 

Missile 

attack 

V(34) 

ALP HAP 

aP« 

Missile  pitch  angle  of 

deg 

Missile 

attack 

V(35) 

ALPHAY 

Qv» 

Missile  yaw  angle  of 

deg 

Mis  sile 

attack 

V(39) 

Q 

q,  Dynamic  pressure 

lb/ft2 

V(40) 

VM 

Total  missile  velocity. 

ft/ sec 

V(41) 

AM 

Missile  Mach  number 

V(44) 

D  WX 

(0 

Scalar  components  of 

V(45) 

DWY 

U) 

y 

missile  angular 
acceleration  in 

rad/ sec^ 

Missile 

V(46) 

D  WZ 

“z 

autopilot  axes 

C(73) 

S 

S,  Missile  ref.  area 

ft2 

C(74) 

D 

d,  Missile  ref.  diameter 

ft 

C(75) 

PSL 

Psl»  Sea  level  pressure 

lb/ft2 

C  ( 76 ) 

TBOOST 

V 

Booster  burn  time 

sec 

C(77) 

TSUST 

ts’ 

Sustainer  burn  time 

sec 

C{78) 

ELTB 

1  TB 

,  Tail  length  (burnout) 

ft 
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TABLE  XXIII.  AERODYNAMICS  (MAERO)  SUBROUTINE  (CONCLUDED) 


Name 

Quantity 

Units 

Coordinate 

System 

C.  ( 7  9 )  XBAR 

X,  eg  to  control  surface 
trailing  edge  distance 

ft 

C(  80)  CLP 

Ci  ,  Roll  damping 
r  coefficient 

1/rad 

C(81 )  AE 

Ae,  Nozzle  exit  area 

ft2 

C ( 8 2 )  A JXO 

Jxoi  Launch  roll  inertia 

slug-ft2 

C ( 83 )  A JX 1 

Jxj ,  End-of-boost  roll 
inertia 

slug-ft2 

C(84)  A  JXT 

Jx{.,  End-of- sustain  roll 
inertia 

slug-ft^ 

C  ( 8  5 )  A  JYO 

Jyo,  Launch  lateral  inertia 

slug-ft^ 

C(86)  A JY 1 

Jyl,  End-of-boost  lateral 
inertia 

slug-ft^ 

C(87)  AJYT 

Jyt,  End-of- sustain  lateral 
inertia 

slug-ft^ 

C(88)  AMASSO 

Mqi  Launch  mass 

slugs 

C(89)  AMASS1 

Mj,  End-of-boost  mass 

slugs 

C(90)  AMASST 

Mt,  End-of- sustain  mass 

slugs 

C(91)  TSEPAR 

Burn  time  prior  to  launch 

sec 

C(92)  TGALT 

Target  altitude 

ft 

Inertial 
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TABLE  XXIV.  AERODYNAMICS  SUBROUTINE 


Intermediate  Expressions  (A) 

C1 

~  CDco  +  ^CDF1  f  CDF2  h)  h  +  1  CDB 

A 

— 
s  1 

6t 

=  6p  cos  -  6y  sin 

6r 

=  6y  cos  ^  +  6p  sin  *a 

COS 

4*a 

2  2 

=  cos  (4^  )  =  1-8  sin  0  cos  <f> 

E  EL  cL 

CO 

H* 

a 

4* 

=  5 . 1 8  (  | c o s  ^a  |  -  |sin^a|)cos^>a 

sin  6 
ra 

Secondary  Expressions  (B) 

Cmo 

= 

2  3 

+  (m^  +  cos  4^a)  a  +  m4  a 

CNo 

= 

NjOr 

+  (N2  +  N3  cos  40a)  a  +  N4  a3 

AC 

m 

= 

(m5 

2 

+  (m^  +  rri^  cos  4<#>a)  a  )  &T 

C 

no 

= 

(nj  c 

2  3 

+  n_  «  )  sin  4 <f> 

L  a 

C 

yo 

= 

(y !  o 

2  3 

+  y,  a  )  sin  4# 

L  E 

AC 

n 

= 

2 

(n3  +  (n4+n5  cos  4^&)  a  )  6R  +  ((nfe  +  n7  cos 

4*a>  «+n8o2  6T)6a 

C 

CO 

= 

C1  +  CZa  +  C3  ° 

AC 

c 

= 

(C4or 

4  Cg  bj,)  bj, 

C, 

to 

= 

(l1  + 

t^a  +  l3  »2  +  £4  a3)  a2  sin  4^& 

ACj 

= 

(!5  +  (*6  +  l7  cos  4^)  a)  6a  +  (ig o  +(l9+  l]()cos  40a>o2)  6R 

d 

2V 

cm  • 

s 

Tf  (N1  +  "(N2  4  N3  cos  4*a  +  “N4»  (fjf] 
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TABLE  XXIV.  AERODYNAMICS  SUBROUTINE  (CONCLUDED) 


Primary  Expressions  (C) 


C  =  C  +  AC 

C  CO  c 


Ct  =  C4  q  +  AC| 


TABLE  XXV.  AERODYNAMIC  COEFFICIENTS  IN  MISSILE 

BODY  AXES 


mb 

Cnb 

Cyb 


'Nb 


C  cos  d  +  C  sin  # 
m  a.  n  a 


C  cos  4  -  C  sin  6 
n  a  m  a 


C  cos  d  +  C  sin 
y  a  n  a 


C.T  cos  6  -  C  sin  6 
N  a  y  a 
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TABLE  XXVI.  EQUATIONS  OF  MOTION  IN  MISSILE  BODY  AXES 


TABLE  XXVII.  FUNCTION  GENERATORS  IN  AERODYNAMIC 

SUBROUTINE 


FGEN  No. 


Input 


Output 


Symbol 


Units 


I  1 


ALT 


SVEL 


Sonic  velocity 


I  8 


I  10 
I  11 
I  12 
I  13 
I  14 
I  15 
I  16 
I  17 
I  18 
I  19 


I  22 
I  23 
I  24 
I  25 
I  27 
I  28 
I  26 
I  30 
I  31 


CDW 


CDF1 

CDF2 

CDB 


'DF2 


ft/ sec 
1/deg 
1  /deg2 
1  /deg2 
1  /deg3 
1  /deg 
1  /deg3 
1  /deg3 
1  /deg2 
1  /deg2 
1  /deg 
1  /deg3 
1  /deg3 
1/deg2 
1  /deg2 
1  /deg4 
1  /deg 
1  /deg2 
1  /deg2 
1  /deg3 
1  /deg2 
1  /deg3 


1  /deg 
1  /deg' 
1/ft 


1  /  deg 
1  /deg' 
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TABLE  XXVII.  FUNCTION  GENERA  TORS  IN  AERODYNAMICS 
SUBROUTINE  (CONCLUDED) 


FGEN  No. 

Input 

Output 

Symbol 

Units 

I 

32 

AM 

SL1 

*1 

1 /deg2 

I 

I 

33 

34 

AM 

AM 

SL2 

SL3 

lZ 

i3 

1 /deg3 

1 /deg4 

I 

35 

AM 

SL4 

*4 

1 /deg5 

I 

36 

AM 

SL5 

'5 

1  /deg 

I 

37 

AM 

SL6 

4 

1 /deg2 

I 

38 

AM 

SL7 

h 

1 /deg2 

I 

39 

AM 

SL8 

<8 

1 /deg2 

I 

40 

AM 

SL9 

i9 

1 /deg3 

K 

40 

AM 

SL1 0 

*10 

1  /deg 

I 

41 

TIME 

DL0D 

Ai/d 

— 

I 

42 

TIME 

0MDL0L 

1  -  Af/|  TB 

— 

I 

44 

TIME 

THRUST 

Rocket  motor  thrust 

lb 

I 

45 

ALT 

PRES 

Ambient  pressure 

lb/ft2 

I 

43 

ALT 

RH0 

Air  density 

slug/ft3 

2.3.11  Program  Glossary 

(U)  A  master  program  glossary  defining  all  the  elements  appearing  in 
both  the  V  and  C  arrays  is  contained  in  Tables  XXVIII  and  XXIX,  respectively 

2.4  TYPE  OF  SIMULATION  FACILITIES  TO  BE  USED 

(U)  The  only  equipment  that  is  required  is  the  GE-635  or  any  other  digital 
computer  capable  of  compiling  FORTRAN  IV  source  decks.  No  training  equip 
ment  or  mockups  are  required. 


2.5  INSTRUMENTATION 
(U)  None  is  required. 

2.6  DATA  REDUCTION  AND  ANALYSIS  TECHNIQUES 
(U)  No  special  data  reduction  techniques  are  required. 
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TABLE  XXVIII.  MASTER  GLOSSARY,  V  ARRAY 


Name 

Quantity 

Units 

Coordinate 

System 

V(l) 

hm,  Missile  altitude  above  ground 

ft 

Inertial 

V(2) 

6  .  Aileron  deflection  command 

ac 

deg 

Missile 

V(3) 

6  ,  Pitch  deflection  command 

pc 

6  ,  Yaw  deflection  command 

yc 

deg 

Missile 

V(4) 

deg 

Missile 

V(5) 

* 

6  ,  Aileron  deflection 
a 

deg 

Missile 

V(6) 

6  ,  Pitch  deflection 

D 

deg 

Missile 

V  (7 ) 

6  .  Yaw  deflection 

V 

deg 

Missile 

V(8) 

V  ,  Missile  velocity  X-axis 

ft /sec 

Missile 

V(9) 

V  ,  Missile  velocity  Y-axis 

ft/sec 

Mis  sile 

V(10) 

V  ,  Missile  velocity  Z-axis 

ft/ sec 

Missile 

V(ll) 

u>x<  Angular  velocity 

rad/sec 

Missile 

V(12) 

w  ,  Angular  velocity 

rad/sec 

Missile 

V(13) 

u>z,  Angular  velocity 

rad/sec 

Missile 

V  (14) 

Ax,  Propulsion  and  aerodynamic 
acceleration 

g 

Missile 

V  ( 15) 

Ay,  Propulsion  and  aerodynamic 
acceleration 

g 

Missile 

V(16) 

Az,  Propulsion  and  aerodynamic 
acceleration 

g 

Missile 

V(17) 

A  ,  Elevation  maneuver  command 
zc 

g 

Autopilot 

V(I8) 

A  ,  Azimuth  maneuver  command 
yc 

g 

Autopilot 

V(19) 

4*,  Euler  yaw  angle 

rad 

V(20) 

<#>,  Euler  roll  angle 

rad 

V(2 1 ) 

0,  Euler  pitch  angle 

rad 

V(22) 

R^,  Seeker  boresight  range 

ft 

Seeker 

V(23) 

R  ,  Seeker  lateral  range 

ft 

Seeker 

V(24) 

Rz,  Seeker  normal  range 

ft 

Seeker 

V(25) 

€  ,  Tracking  error  angle,  pitch 

rad 

Seeker 

V(26) 

e  ,  Tracking  error  angle,  yaw 

rad 

Seeker 

V(27) 

v9  Seeker  elevation  gimble  angle 

rad 

V(28) 

t].  Seeker  azimuth  gimble  angle 

rad 
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TABLE  XXVU1.  MASTER  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

Coordinate 

System 

V(29) 

Ri,  Horizontal  longitudinal  range 
component 

ft 

Inertial 

V(30) 

Rj,  Horizontal  lateral  range 
component 

ft 

Inertial 

V(3 1) 

£  ,  Gate  error  angle,  pitch 

rad 

Seeker 

V(32) 

c  #  Gate  error  angle,  yaw 

rad 

Seeker 

V(33) 

a,  Total  missile  angle  of  attack 

deg 

Missile 

V(34) 

a  ,  Missile  pitch  angle  of  attack 

deg 

Missile 

V(35) 

a  ,  Missile  yaw  angle  of  attack 

deg 

Missile 

V(36) 

Vi,  Horizontal  longitudinal  velocity 
component 

ft/ sec 

Inertial 

V(37) 

Vj,  Horizontal  lateral  velocity 
component 

ft/sec 

Inertial 

V(38) 

Vfc,  Vertical  velocity  component 

ft/ sec 

Inertial 

V(39) 

q.  Dynamic  pressure 

lb/ft2 

V(40) 

Total  missile  velocity 

ft/sec 

V  (4 1 ) 

Missile  Mach  number 

V(42) 

a^.  Acceleration  command  pitch 

g 

Autopilot 

V(43) 

a  ,  Acceleration  command  yaw 

g 

Autopilot 

V(44) 

V(45) 

V(46) 

Scalar  components  of  mis- 

^y  sile  angular  acceleration 

in  missile  axes 
z 

rad/sec^ 

Missile 

V(47) 

6  ,  Aileron  command  rate 

ac 

deg/sec 

Missile 

V(48) 

h  ,  Elevator  command  rate 
pc 

6  ,  Rudder  command  rate 

°yc 

Closest  approach  at  end  of  flight 

deg/sec 

Missile 

V(49) 

deg/sec 

Missile 

V(50) 

ft 

V(51) 

Range  component  in  Y  seeker  axis 

ft 

Seeker 

V(52) 

Range  component  in  Z  seeker  axis 

ft 

Seeker 

V(53) 

V(54) 

V(55) 

wx 

Missile  body  rates  in 
wy  autopilot  axes 

co- 

z 

rad/sec 

Autopilot 

TABLE  XXVIII.  MASTER  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

Coordinate 

System 

V(56) 

V(  57 ) 

V(  58 ) 

d)' 

X 

(I)1 

y 

(j 1 

7. 

Scalar  components  of  missile 
>  angular  acceleration  in 

autopilot  axes 

rad/sec^ 

Autopilot 

V(59) 

V(60) 

V(6 1 ) 

>  >  > 

*  -*<  -  x  - 
_ 

Propulsive  and  aerodynamic 
acceleration  components  in 
autopilot  axes 

g 

Autopilot 

V(62) 

V(63) 

V(64) 

Vxl 

V' 

y 

V' 

z 

1 

Missile  velocity  components 
in  autopilot  axes 

ft/  sec 

Autopilot 

V(65) 

Special  test  variable  -  used  as 
system  diagnostic 

V(66) 

Total  miss  distance 

ft 

Miss  Distance 

V(67) 

x  component  of  range 

Autopilot 

V(68) 

y  component  of  range 

ft 

Autopilot 

V(69) 

z  component  of  range 

ft 

Autopilot 

V(70) 

y  component  of  miss 

ft 

Miss  Distance 

V(7 1 ) 

z  component  of  miss 

ft 

Miss  Distance 

V(72) 

x  component  of  acceleration 

g 

Inertial 

V(73) 

y  component  of  acceleration 

g 

Inertial 

V(74) 

z  component  of  acceleration 

g 

Inertial 

V(75) 

y  component  of  acceleration  at 
blind  range 

g 

Miss  Distance 

V(76) 

z  component  of  acceleration  at 
blind  range 

g 

Miss  Distance 

V(77) 

Blind  time  in  yaw  channel 

sec 

V(78) 

Blind  time  in  pitch  channel 

sec 

V(79) 

Final  line  of  sight  angle  (vertical) 

rad 

Inertial 

V(80) 

Final  heading  angle  (horizontal) 

rad 

Inertial 

V(8 1 ) 

x  component,  LOS  rate 

rad/sec 

Inertial 

V(82) 

y  component,  LOS  rate 

rad/sec 

Inertial 

V(83) 

z  component,  LOS  rate 

rad/ sec 

Inertial 

V(84) 

A  ,  Guidance  gain 

J 
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TABLE  XXVIII.  MASTER  GLOSSARY,  V  ARRAY  (CONTINUED) 


Name 


V(87) 

V(88) 

V(90) 

V(91) 

V(92) 

V(93) 

V(94) 

V(95) 

V(96) 

V(97) 

V(98) 

V(99) 

V(100) 
V(I01) 
V(102) 
V(103) 
V(I04) 
V(105) 
V(  1 06) 
V(  1 07) 
V ( 108) 
V(I09) 
V(110) 


DEXS 


DG1 


DG2 


DIN 


G2N 


FFE 

DFE 


FEXS 

DFEXS 


Quantity 


Total  yaw  precession  rate 

Total  pitch  precession  rate 

Yaw  gyro  inertial  angle 

Yaw  look  angle  (indicated) 

Forcing  function  cross- 
coupled  equation  1 

Derivative  forcing  function 
cross -coupled  equation  1 

Forcing  function  cross- 
coupled  equation  2 

Derivative  forcing  function 
cross -coupled  equation  2 

Integral  forcing  function 
cross -coupled  equation  1 

Integral  forcing  function 
cross -coupled  equation  2 

Forcing  function  yaw  axis 

Derivative  forcing  function 
yaw  axis 

Forcing  function  pitch  axis 

Derivative  forcing  function 
pitch  axis 


Not  Used 
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TABLE  XXVIII.  MASTER  GLOSSARY,  V  ARRAY  (CONCLUDED) 


Name 

Quantity 

Units 

V(lll) 

Sum  1  -  Tracker  sampler  bias 

sec 

V  ( 1 12) 

TEAYD  Tracker  error  yaw  • 
RKAMG 

deg 

V(1 13) 

TEAPD  Tracker  error  pitch  • 
RKAMG 

deg 

V( 1 14) 

TEAYS  -  Tracker  Z0H  output 
signal,  yaw 

deg 

V( 1 15) 

TEAPS  -  Tracker  Z0H  output 
signal,  pitch 

deg 

V(1 16) 

VSYP  -  Tracker  output  signal 
pitch 

deg/sec 

V( 1 17) 

VSPP  -  Tracker  output  signal  yaw 

deg/ sec 

V(  1  18) 

TEYD  -  Tracking  error  -  yaw 

deg 

V(  1 19) 

TEPD  -  Tracking  error  -  pitch 

deg 

V(  1 20) 

UND  Seeker  elevation 

deg  1 

V(  1  21 ) 

ETAD  Seeker  azimuth 

deg 

V( 1 22) 

WXD 

V( 1 23) 

WYD  Missile  angular  velocity 

deg 

V(  1 24 ) 

WZD 

V(  1 25) 

XLOSD 

V(  1  26) 

YLOSD  LOS  Rate,  inertial 

deg/sec 

V(  1 27) 

ZLOSD 

V(  1 28) 

ANT  (New) 

V(129) 

^  Yaw  D 

V(  1 30) 

Roll  D  Error  angle 

deg 

V(  1 3 1 ) 

9  Pitch 

V(  1 32) 

DED  Total  precession  rate,  yaw 

V(  1 33) 

DEXSD  Total  precession  rate, 
pitch 

deg/sec 
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TABLE  XXIX.  MASTER  GLOSSARY,  C  ARRAY* 


Name 

Quantity 

Units 

C(l) 

Not  used 

Autopilot  activation  time 

sec 

Blind 

range,  pitch 

ft 

C(4) 

Gate  error  angle /^-  V  (18) 

mills 

g/sec 

C(5) 

Blind 

range  filter  time  constant 

sec 

C(6) 

Kp,  Angle  restoration  gain 

g/deg 

C  ( 7 ) 

Restoration  angle  bias 

deg 

C(8) 

To,  Angle  restoration  filter  time  constant 

sec 

C(9) 

Blind 

range ,  yaw 

ft 

C(10) 

1 

C(  1 1) 

S  Not  used 

C(  12) 

J 

I 

C(13) 

Half  field  of  view,  seeker 

rad 

C(  14) 

Guidance  gain 

g/deg/sec 

C(  15) 

T 1  ,  Tracker  time  constant 

sec 

C(16) 

Precession  rate  limit 

rad/ sec 

C(  17) 

C 1  "I 

rad / sec 

C(  18) 

C2 

rad/sec 

C(  1 9 ) 

C3 

rad/ sec/g 

C(20) 

=4 

C(21) 

C5 

1  /sec 

C(22) 

C6 

1  /  sec 

C(23) 

C7 

)  Seeker  drift  terms 

sec 

C(24) 

S 

sec 

C(25) 

C9 

rad/ sec/ g 

C(26) 

C 10 

rad/sec /g 

C(27) 

C11 

1  /sec/g 

C(28) 

C 1 2 

2 

rad / sec/g 

C  ( 29 ) 

C13  - 

2 

rad/sec/g 

*See  Appendix  V  -  Volume  III  -  System  Analysis  Document  for  typical  values 
of  the  C  array. 
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TABLE  XXIX.  MASTER  GLOSSARY,  C  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

COO) 

C14^ 

rad/sec/g^ 

C(3 1 ) 

°,s 

rad/sec 

C(32) 

C16 

rad/sec 

C(33) 

C17 

rad/ sec/g 

C(34) 

C18 

k  Seeker  drift  terms 

C(35) 

C19 

/ 

1  /sec 

C(36) 

C20 

sec 

C(37) 

C21 

rad/sec/g 

C(38) 

C22 

rad/sec/g  1 

2 

C(39) 

C23  > 

rad/sec/g 

C(40) 

Kj,  Trackin  Toop  velocity  gain 

1  /sec 

C(4 1 ) 

TG,  Gimbal  load 

g 

C(42) 

Drift  control  vset  to  1  .0  to  include  drift) 

C(43) 

Autopilot  lateral  channel  activation  switch 
level 

g 

C(44) 

K  ,  Autopilot  lateral  channel  gain 

deg/g 

C(45) 

K^,  Autopilot  damping  gain 

deg/deg/sec 

C(46) 

0^,  Roll  channel  gain 

deg/deg 

C(47) 

Acceleration  lir  lateral  channels 

g 

C(48) 

Command  limit,  lateral  channels 

deg 

C(49) 

T  ,  Lead  time  constant 
sn 

sec 

C(50) 

Xg{j,  Lag  time  constant 

sec 

C(51 ) 

Tt_  »  Lead  time  constant 
bn 

sec 

C(  52) 

tYDLG’  time  constant 

sec 

C(53) 

Roll  channel  deadtime 

sec 

C(54) 

Roll  rate  signal  limit 

rad/sec 

C(55) 

Tra,  Time  constant  in  autopilot 

sec 

C(  56) 

Roll  rate  switch  level 

rad/ sec 

C(57) 

Tacc  * 

Lateral  channel  time  constant 

sec 

C(  58) 

tRG’ 

Lateral  channel  time  constant 

sec 
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TABLE  XXIX.  MASTER  GLOSSARY,  C  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

■■ 

Roll  rate  gain 

deg/deg/sec 

•  H 

Roll  channel  limit  change  time 

sec 

C(6l) 

t  , ,  Lead  time  constant 
r  1 

sec 

C(62) 

rr2*  Lag  time  constant 

sec 

C(63) 

Roll  command  limit 

deg 

C(64) 

Roll  command  limit 

deg 

C(65) 

Lateral  channel  rate  limit 

rad/sec 

C(66) 

C(67) 

C(68) 

C(69) 

^  Not  used 

C(70) 

C  (7 1 ) 

C(72) 

J 

C(73) 

S,  Reference  area 

ft2 

C(74) 

d,  Reference  diameter 

ft 

C(75) 

Sea  level  pressure 

lb /ft2 

C(76) 

Booster  burn  time 

sec 

C(77) 

Sustainer  burn  time 

sec 

C(78) 

Tail  length,  burnout 

1  D 

ft 

C(79) 

x,  eg 

to  flipper  distance 

ft 

C(80) 

C*  ,  Roll  damping  coefficient 

Xp 

1  /  rad 
,  2 

C(81 ) 

A  ,  Nozzle  exit  area 
e 

ft 

C(82) 

J  ,  Launch  roll  inertia 

slug-ft2 

xo 

2 

C(83) 

J  ,  ,  End-of-boost  roll  inertia 
xl 

slug-ft 
.  2 

C(84) 

J  End-of- sustain  roll  inertia 

xt 

slug-ft 

C(85) 

J  ,  Launch  lateral  inertia 

slug-ft£ 

yo 

2 

C(86) 

J  ,,  End-of-boost  lateral  inertia 

yi 

slug-ft 
.  2 

C(87) 

J  End-of-sustain  lateral  inertia 

yt 

slug-ft 

C(88) 

M  ,  Launch  mass 
o 

slugs 
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TABLE  XXIX.  MASTER  GLOSSARY,  C  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

m 

M.,  End-of-boost  mass 

i 

slugs 

ESI 

Mt,  End-of-sustain  mass 

slugs 

C(9 1 ) 

Burn  time  prior  to  launch 

sec 

C(92) 

Target  altitude 

ft 

C(93) 

C(94) 

C(95) 

Not  used 

C(96) 

C(97) 

Control  surface  servo  gain 

1  /sec 

C(98) 

Control  surface  time  constant 

sec 

C(99) 

Control  surface  velocity  limit 

deg/ sec 

C(100) 

Control  surface  angle  limit 

deg 

C(101) 

Control  surface  rate  threshold 

deg/sec 

EBBI 

Not  used 

BBSS 

Dual  purpose  input: 

Min.  range  of  aimpoint  wander  when 

ft 

C ( 1 06 )  i  0.  Initial  target  accelera¬ 
tion  when  C(  106)  =  0.  Subroutine 

g 

bypassed  when  C ( 103)  =  0. 

C ( 1 04 ) 

Dual  purpose  input: 

A. 

used  in  aimpoint  wander  when 

C(106)  i  0.  Final  target  velocity 
when  C(  106)  -  0 . 

ft/ sec 

C(  105) 

Dual  purpose  input: 

PL0TK,  used  in  aimpoint  wander 
when  C(106)  i  0.  Boresight  range  at 
which  target  motion  starts  when 

C(  106)  =  0. 

ft 

C( 106) 

PH0T0K,  used  in  aimpoint  wander 

C(107) 

C  ( 1 08) 

Not  used 

C(109) 

LOS  memory  threshold  in  guidance  law 

deg 

C(110) 

G-Bias 

g 

C(lll) 

LOS  memory  gain  in  guidance  law 

g/deg 
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TABLE  XXIX.  MASTER  GLOSSARY ,  C  ARRAY  (CONTINUED) 


Name 

Quantity 

Units 

C  ( 1 12) 

C  ( 1 13) 

C(  1 14) 

TB  Sampling  time  (function  of  distribution 
to  target  (0 . 0800) 

C  ( 1 15) 

TC  Sampling  time  control  (0.1000) 

C( 1 16) 

3S  -  Gyro  rotor  speed 

rad/ sec 

C  ( 1 17) 

K2T  -  Precession  torque  coefficient 

gem/ V 

C  ( 1 18) 

Dump  program  control  logic 

B  =  0 

C ( 1 19) 

-  Rail  control  logic 

S  =  1 .0 

C( 1 20) 

C(  1 2 1 ) 

C(  1 22) 

C(  1 23) 

C(  1 24) 

C(  1 25) 

C(  1 26) 

C(  1 27) 

C(  1 28) 

C(  1 29) 

C(  1 30) 

C(  1 3  1 ) 

C( 1 32) 

Not  used 

C(  1 33) 

PSIPGJ  LOS  Memory  threshold  guidance  law 

deg 

C(  1 34) 

BP  -  G-Bias 

g 

C(  1 35) 

AKSIGP  LOS  Memory  gain  in  guidance  law 

g/deg 

C(136) 

GNUT  -  Program  logic  control  -  W/0-0, 

W  =  1 . 0 

C(  1 37) 

DFR  Coulomb  friction  drift  factor 

C(  138) 

DST  Spring  torque  drift  factor 

C(  139) 

DSU  Unbalance  drift  factor 

C(140) 

DAN  Anisoelastic  drift  factor 
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TABLE  XXIX.  MASTER  GLOSSARY,  C  ARRAY  (CONCLUDED) 


Name 

Quantity 

Units 

C(  141 ) 

DDU  Dynamic  unbalance  factor 

HEX 

SK  Torquer  gain  coefficient 

V/deg/sec 

mm 

AKT  -  Tracker  gain  constant 

1  sec 

C  ( 1 44 ) 

TS  -  Sampling  period 

sec 

C(  145) 

OMEGLD  -  Precession  rate  limit 

deg/ sec 

C(  146) 

GKK  -  Guidance  gain 

g/deg/sec 

C(  147) 

BIAS  -  Sampling  rate  offset  bias 

sec 

C(  148) 

TLDP  -  Tracker  filter  lead  time  constant  - 
Pitch 

sec 

C  ( 1 49 ) 

TLGP  -  Tracker  filter  lag  time  constant  - 
Pitch 

sec 

C(  1 50) 

TLOY  -  Tracker  filter  lead  time  constant  - 
Yaw 

sec 

C(  1 51 ) 

TLGY  -  Tracker  filter  lag  time  constant  - 
Yaw 

sec 

C(  1 52) 

SPOT  -  Tracker  spot  size 

Ft 

C(  1 53) 

CF1  Friction  factor  coefficient 

C( 1 54) 

CF2  Friction  factor  coefficient 

C(  1 55) 

C(  1 56) 

C(  1 57) 

CF3  Friction  factor  coefficient 

C( 1 58) 

C(1 59) 
C(l60) 

Not  used 
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INITIAL  DISTRIBUTION 


HQ  USAF  (RDP\)  1 

HQ  USAF  (RDQRM)  1 

HQ  USAF  (XOOW)  2 

A?  SC  (DLP)  1 

AFSC  (DLTW)  1 

AFSC  (DLXP)  1 

AFSC  CSDWM)  1 

TAC  (DRPW)  1 

’  AFAL  (TE)  1 

AFAL  (NV)  1 

ASD  (SD)  2 

*  ASD  (ENYY)  2 

FTD  (PDJR)  5 

AFFDL  (FYS)  1 

RADC  (IRAP)  1 

AUL  (AUL/LSE-70-239)  1 

SAND I A  CORP  1 

ARMY  AIR  MOBILITY  DIV  2 

DEP  CH  STAFF  USMC  CODE  AAW-IC  1 

CH  NAV  OPP.S  (OP-982E21)  1 

OODR5E,  OAD/COMBAT  SUPPORT  1 

WPNS  SYS  EVAL  GP  1 

ARPA  (TIO)  1 

NAV  AIR  SYS  COMD,  AIR-350B  1 

DASA  INFO  S  ANAL  CTR  GE  CO-TEMPO  1 

RAND  CORP  1 

NWC  (CODE  456)  2 

NWC  (CODE  753)  1 

REDSTONE  SCI  INFO  CTR  1 

USA  ENGR  R5D  LABS  1 

BATTELLE  MEMORIAL  INST  1 

INST  FOR  DEF  ANALYSES  1 

DDC  2 

DL  1 

DLB  1 

DLW  1 

DLU  1 

DLX  1 

DLOSL  2 

AFSC  (DLW)  1 

t  ASD  (SD-65)  1 

DLWS  1 

PACAF  (IGY)  1 

•  ASD  (ENYS)  1 
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